Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,497,705 T→G 76.4% R20R (CGT→CGG gldA → glycerol dehydrogenase, NAD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,497,7050TG76.4% 98.2 / 12.7 38R20R (CGT→CGGgldAglycerol dehydrogenase, NAD
Reads supporting (aligned to +/- strand):  ref base T (3/6);  new base G (20/9);  total (23/15)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GAGCAATTATGGACCGCATTATTCAATCACCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGTCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAGTGGTGGGTGACAAATTTGTTTTAGG  >  W3110S.gb/3497638‑3497770
                                                                   |                                                                 
gAGCAATTATGGACCGCATTATTCAATCACCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGGCTg                                                                >  1:2253193/1‑71 (MQ=255)
 aGCAATTATGGACCGCATTATTCAATCACCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGGCTgg                                                               >  1:2682442/1‑71 (MQ=255)
 aGCAATTATGGACCGCATTATTCAATCACCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGGCTgg                                                               >  1:1501097/1‑71 (MQ=255)
 aGCAATTATGGACCGCATTATTCAATCACCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGGCTg                                                                >  1:1179938/1‑70 (MQ=255)
 aGCAATTATGGACCGCATTATTCAATCACCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGGCTg                                                                <  1:1658108/70‑1 (MQ=255)
    aaTTATGGACCGCATTATTCAATCACCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGTCTGGGCg                                                            <  1:445732/71‑1 (MQ=255)
    aaTTATGGACCGCATTATTCAATCACCGGGTAAATACATACAGGGCGCTGATGTGATTAATCGTCTGGGCg                                                            <  1:1697967/71‑1 (MQ=255)
        aTGGACCGCATTATTCAATCACCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGGCTGGGCGAATa                                                        >  1:1381138/1‑71 (MQ=255)
         tGGACCGCATTATTCAATCACCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGTCTGGGCGAAt                                                         <  1:2352044/69‑1 (MQ=255)
            aCCGCATTATTCAATCACCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGGCTGGGCGAATACCTg                                                    <  1:981087/71‑1 (MQ=255)
            aCCGCATTATTCAATCACCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGGCTGGGCGAATACCTg                                                    <  1:1887743/71‑1 (MQ=255)
             ccGCATTATTCAATCACCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGTCTGGGCg                                                            <  1:2126163/62‑1 (MQ=255)
                    attCAATCACCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGGCTGGGCGAATACCTGAAGCCGCt                                            >  1:307477/1‑71 (MQ=255)
                     ttCAATCACCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGGCTGGGCGAATACCTGAAGCCGc                                             >  1:1875081/1‑69 (MQ=255)
                            aCCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGGCTGGGCGaa                                                          >  1:2174965/1‑49 (MQ=255)
                                  tAAATACATCCAGGGCGCTGATGTGATTAATCGTCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGt                              <  1:711325/71‑1 (MQ=255)
                                  tAAATACATCCAGGGCGCTGATGTGATTAATCGGCTGGGCTAAt                                                         <  1:1980318/44‑1 (MQ=38)
                                  tAAATACATCCAGGGCGCTGATGTGATTAATCGGCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGt                              <  1:798812/71‑1 (MQ=255)
                                    aaTACATCCAGGGCGCTGATGTGATTAATCGGCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGtt                             <  1:1013399/70‑1 (MQ=255)
                                    aaTACATCCAGGGCGCTGATGTGATTAATCGGCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGt                              >  1:2297025/1‑69 (MQ=255)
                                     aTACATCCAGGGCGCTGATGTGATTAATCGGCTGGGCg                                                            <  1:1494248/38‑1 (MQ=255)
                                      tACATCCAGGGCGCTGATGTGATTAATCGGCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAgt                          >  1:398531/1‑71 (MQ=255)
                                      tACATCCAGGGCGCTGATGTGATTAATCGGCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAgt                          <  1:1757910/71‑1 (MQ=255)
                                       aCATCCAGGGCGCTGATGTGATTAATCGGCTGGGCGaa                                                          >  1:1457434/1‑38 (MQ=255)
                                       aCATCCAGGGCGCTGATGTGATTAATCGGCTGGGCGa                                                           >  1:1724982/1‑37 (MQ=255)
                                            catggCGCTGATGTGATTAATCGGCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAGTggtggg                    >  1:1699536/4‑71 (MQ=255)
                                             aGGGCGCTGATGTGATTAATCGTCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAGTggtgggt                   <  1:391217/71‑1 (MQ=255)
                                              gggCGCTGATGTGATTAATCGTCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAGTggtgggt                   >  1:831186/1‑70 (MQ=255)
                                              gggCGCTGATGTGATTAATCGTCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAGTggtggg                    >  1:626971/1‑69 (MQ=255)
                                                  gcTGATGTGATTAATCGGCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAGTGGTGGGTGACa               >  1:1084986/1‑70 (MQ=255)
                                                        gtgATTAATCGTCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAGTGGTGGGTGACAAATTTGt        >  1:1403352/1‑71 (MQ=255)
                                                        gtgATTAATCGGCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAgtgg                        >  1:450509/1‑55 (MQ=255)
                                                         tgATTAATCGGCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAGTGGTGGGTGACAAATTTGtt       >  1:947120/1‑71 (MQ=255)
                                                           aTTAATCGGCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAgtg                         >  1:2532139/1‑51 (MQ=255)
                                                           aTTAATCGGCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAgtg                         >  1:2528812/1‑51 (MQ=255)
                                                              aaTCGGCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAGTGGTGGGTGACAAATTTGTTTTa    >  1:1535111/1‑69 (MQ=255)
                                                              aaTCGGCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAGTGGTGGGTGACAAATTTGTTTTAgg  >  1:2121848/1‑71 (MQ=255)
                                                                tCGGCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGtt                             <  1:1503438/42‑1 (MQ=255)
                                                                   |                                                                 
GAGCAATTATGGACCGCATTATTCAATCACCGGGTAAATACATCCAGGGCGCTGATGTGATTAATCGTCTGGGCGAATACCTGAAGCCGCTGGCAGAACGCTGGTTAGTGGTGGGTGACAAATTTGTTTTAGG  >  W3110S.gb/3497638‑3497770

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: