Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,353,829 G→A 12.9% T26M (ACG→ATG)  yhcC ← predicted Fe‑S oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,353,8290GA12.9% 84.2 / 4.5 31T26M (ACG→ATG) yhcCpredicted Fe‑S oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (9/18);  new base A (2/2);  total (11/20)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

ACAGAATGTGCAGCCACCACGCCCGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCTTATGCACCTTTTGCCCATAACGACGGGTGAGATCACCACCAAACATATTGACTAATTTCTGT  >  W3110S.gb/3353762‑3353897
                                                                   |                                                                    
aCAGAATGTGCAGCCACCACGCCCGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGt                                                                     >  1:109944/1‑69 (MQ=255)
 cAGAATGTGCAGCCACCACGCCCGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCa                                                                   >  1:552781/1‑70 (MQ=255)
 cAGAATGTGCAGCCACCACGCCCGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCa                                                                   <  1:150753/70‑1 (MQ=255)
   gAATGTGCAGCCACCACGCCCGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCt                                                                <  1:839039/71‑1 (MQ=255)
     aTGTGCAGCCACCACGCCCGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCtt                                                               <  1:256900/70‑1 (MQ=255)
      tgtgCAGCCACCACGCCCGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCTTAt                                                             >  1:477527/1‑71 (MQ=255)
        tgCAGCCACCACGCCCGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCt                                                                <  1:1637191/66‑1 (MQ=255)
               accacGCCCGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCTTATGCACCtttt                                                    <  1:1184993/71‑1 (MQ=255)
               accacGCCCGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCATCAGCTTATGCACCtttt                                                    <  1:1743527/71‑1 (MQ=255)
                  acGCCCGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCTTATGCACCTTTTGcc                                                 >  1:766987/1‑71 (MQ=255)
                  acGCCCGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCTTATGCACCTTTTGcc                                                 <  1:1654685/71‑1 (MQ=255)
                    gCCCGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGGCAGCTTATGCACCTTTTGCCCa                                               <  1:1765199/71‑1 (MQ=255)
                     cccGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCtt                                                               >  1:526862/1‑54 (MQ=255)
                     cccGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCATCAGCTTATGCACCTTTTGCCCAt                                              <  1:1207812/71‑1 (MQ=255)
                            gTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCTTATGCACCt                                                       <  1:349580/55‑1 (MQ=255)
                            gTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCTTATGCACCTTTTGCCCATAacgacg                                       <  1:1382686/71‑1 (MQ=255)
                             tACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCTTATGCACCTTTTGCCCATAACGACgg                                      <  1:1160808/71‑1 (MQ=255)
                                 gTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCTTATGCACCTTTTGCCCATAACGACGGGt                                    >  1:433508/1‑69 (MQ=255)
                                  tCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCTTATGCACCTTTTGCCCATAACGACGGGTgag                                 <  1:531214/71‑1 (MQ=255)
                                       gTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCTTATGCACCTTTTGCCCATAACGACGGGTGAGATcac                            <  1:1662450/71‑1 (MQ=255)
                                                gCTAAAACCGCCGTGCAGCGTCAGCTTATGCACCTTTTGCCCAt                                              <  1:785931/44‑1 (MQ=255)
                                                 cTAAAACCGCCGTGCAGCATCAGCTTATGCACCTTTTGcc                                                 >  1:1282173/1‑40 (MQ=255)
                                                  tAAAACCGCCGTGCAGCGTCAGCTTATGCACCTTTTGCCCATaa                                            >  1:632441/1‑44 (MQ=255)
                                                  tAAAACCGCCGTGCAGCATCAGCTTATGCACCTTTTGCCCATAACGACGGGTg                                   >  1:1281357/1‑53 (MQ=255)
                                                    aaaCCGCCGTGCAGCGTCAGCTTATGCACCTTTTGCCCATAACGACGGGTGAGATCACCACCAAACata                 <  1:473683/69‑1 (MQ=255)
                                                    aaaCCGCCGTGCAGCGTCAGCTTATGCACCTTTTGCCCATAACGACGGGTGAGATCACCACCAAACata                 <  1:66565/69‑1 (MQ=255)
                                                      aCCGCCGTGCAGCGTCAGCTTATGCACCTTTTGCCCATAACGACGGGTGAGATCACCACCAAACATATTGa             >  1:836225/1‑71 (MQ=255)
                                                              gCAGCGTCAGCTTATGCACCTTTTGCCCATAACGACGGGTGAGATCACCACCAAACATATTGACTAAt        >  1:584737/1‑68 (MQ=255)
                                                                 gCGTCAGCTTATGCACCTTTTGCCCATAACGACGGGTGAGATCACCACCAAACATAtt               <  1:1104849/58‑1 (MQ=255)
                                                                 gCGTCAGCTTATGCACCTTTTGCCCATAACGACGGGTGAGATCACCACCAAACATATTGACTAATTTCTGt  <  1:1685304/71‑1 (MQ=255)
                                                                 gCGTCAGCTTATGCACCTTTTGCCCATAACGACGGGTGAGATCACCACCAAACATATTGACTAATTTCTGt  <  1:1030375/71‑1 (MQ=255)
                                                                   |                                                                    
ACAGAATGTGCAGCCACCACGCCCGATGGTACCGTCACGGTTAGGGCAGCTAAAACCGCCGTGCAGCGTCAGCTTATGCACCTTTTGCCCATAACGACGGGTGAGATCACCACCAAACATATTGACTAATTTCTGT  >  W3110S.gb/3353762‑3353897

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: