@@ -9,9 +9,9 @@ import Poseidon.Core.ColumnTypesUtils (ListColumn (..))
99import Poseidon.Core.Janno (JannoRow (.. ))
1010import Poseidon.Core.Utils (LogA , PoseidonException (.. ),
1111 PoseidonIO ,
12- envInputPlinkMode ,
13- logInfo , logWarning ,
14- logWithEnv , padLeft )
12+ envInputPlinkMode , logInfo ,
13+ logWarning , logWithEnv ,
14+ padLeft )
1515
1616import Control.Monad (forM , forM_ , unless , when )
1717import Control.Monad.Catch (MonadThrow , throwM )
@@ -325,9 +325,9 @@ getConsensusSnpEntry strandCheck snpEntries = do
325325 (i: _) -> return i
326326 _ -> return $ head uniqueIds
327327 genPos <- case uniqueGenPos of
328- [p] -> return p
328+ [p] -> return p
329329 [0.0 , p] -> return p -- 0.0 is considered "no data" in genetic position column
330- _ -> return $ maximum uniqueGenPos -- multiple non-zero genetic positions. Choosing the largest one.
330+ _ -> return $ maximum uniqueGenPos -- multiple non-zero genetic positions. Choosing the largest one.
331331 (ref, alt) <- if strandCheck then getConsensusAllelesStrandCheck snpEntries else return (getConsensusAlleles snpEntries)
332332 return (EigenstratSnpEntry chrom pos genPos id_ ref alt)
333333
@@ -354,7 +354,7 @@ checkAlleleFlipNeeded
354354 else if (refA', altA') == (consAltA', consRefA') then Right True -- alleles flipped
355355 else if strandCheck && (revComp refA', revComp altA') == (consRefA', consAltA') then
356356 -- simple concordance on reverse strand
357- Right False
357+ Right False
358358 else if strandCheck && (revComp refA', revComp altA') == (consAltA', consRefA') then
359359 -- alleles flipped on reverse strand
360360 Right True
@@ -366,7 +366,7 @@ checkAlleleFlipNeeded
366366 if refA' == consRefA' || (strandCheck && revComp refA' == consRefA') then do
367367 validateAllRef
368368 Right False
369- else if altA' == consAltA' || (strandCheck && revComp altA' == consAltA') then do
369+ else if altA' == consAltA' || (strandCheck && revComp altA' == consAltA') then do
370370 validateAllAlt
371371 Right False
372372 else if refA' == consAltA' || (strandCheck && revComp refA' == consAltA') then do
@@ -398,8 +398,8 @@ getConsensusAlleles snpEntries =
398398 let allAlleles = concat $ [[r, a] | EigenstratSnpEntry _ _ _ _ r a <- snpEntries]
399399 uniqueAlleles = nub . filter (not . isMissing) $ allAlleles
400400 in case uniqueAlleles of
401- [] -> (' N' , ' N' )
402- [r] -> (' N' , r)
401+ [] -> (' N' , ' N' )
402+ [r] -> (' N' , r)
403403 (ref : alt : _) -> (ref, alt) -- at this point we don't care if we have more than two alleles,
404404 -- this will be checked in the Recoding step.
405405
0 commit comments