Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,701,075 T→C 30.4% L107L (TTA→TTG rbsC ← D‑ribose transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,701,0750TC30.4% 58.9 / 35.9 46L107L (TTA→TTGrbsCD‑ribose transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (10/22);  new base C (10/4);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 2.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

GCGATAAACGCCTGGACGCGACCTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATGCCGACGATAGATGCAGCAACTGCGC  >  W3110S.gb/3701010‑3701145
                                                                 |                                                                      
gCGATAAACGCCTGGACGCGACCTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCCAAcgc                                                                   >  1:2248864/1‑71 (MQ=255)
gCGATAAACGCCTGGACGCGACCTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCCAAcgc                                                                   >  1:1945899/1‑71 (MQ=255)
 cGATAAACGCCTGGGCGCGACCTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCTAAcgc                                                                   >  1:4327877/1‑70 (MQ=255)
     aaaCGCCTGGACGCGACCTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGagag                                                              >  1:1886819/1‑71 (MQ=255)
     aaaCGCCTGGACGCGACCTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCCAACGCga                                                                 >  1:895918/1‑68 (MQ=255)
     aaaCGCCTGGACGCGACCTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCCAACGCGagag                                                              >  1:2853903/1‑71 (MQ=255)
        cGCCTGGACGCGACCTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGagcag                                                           <  1:307502/71‑1 (MQ=255)
        cGCCTGGACGCGACCTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGagcag                                                           <  1:1761896/71‑1 (MQ=255)
              gACGCGACCTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGagcagcagc                                                       <  1:2438725/69‑1 (MQ=255)
              gACGCGACCTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCg                                                      >  1:3259539/1‑70 (MQ=255)
              gACGCGACCTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGa                                                     >  1:2626271/1‑71 (MQ=255)
              gACGCGACCTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGa                                                     <  1:1355139/71‑1 (MQ=255)
              gACGCGACCTTTCGCTACAATCACCCCGGTAACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGa                                                     >  1:1845479/1‑71 (MQ=255)
                   gACCTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGc                                                 >  1:4500258/1‑70 (MQ=255)
                      cTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCTCCCAACGCGAGAGCAGCAGCGACAGccacc                                             <  1:336625/71‑1 (MQ=255)
                      cTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCTa                                                                       <  1:2763347/45‑1 (MQ=255)
                       tttCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCCAACGCGAGAGCAGCAGCGa                                                     >  1:3051819/1‑62 (MQ=255)
                        ttCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGagcagcag                                                        >  1:3383618/1‑58 (MQ=255)
                             tACAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCga                                                                 >  1:3593241/1‑44 (MQ=255)
                               cAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCAtt                                     <  1:2854155/70‑1 (MQ=255)
                               cAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTg                                    <  1:835413/71‑1 (MQ=255)
                               cAATCACCCCGGTTACCGCACCAATTGCGGCACCCAACGCGAGAGCAGCAGCGACAGCCACCAGCGCAtt                                     <  1:4419102/70‑1 (MQ=255)
                                 aTCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGAc                                  <  1:4446267/71‑1 (MQ=255)
                                  tCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACt                                 <  1:2948581/71‑1 (MQ=255)
                                  tCACCCCGGTTACCGCACCAATTGCGGCACCCAACGCGAGAGCAGCAGCGCCAGCCACCAGCGCATTGa                                   >  1:1295865/1‑69 (MQ=255)
                                    aCCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTc                               <  1:2170428/71‑1 (MQ=255)
                                    aCCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTc                               <  1:3720405/71‑1 (MQ=255)
                                    aCCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTc                               <  1:3221744/71‑1 (MQ=255)
                                    aCCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTc                               <  1:1606238/71‑1 (MQ=255)
                                        cGGTTACCGCACCAATTGCGGCACCCAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTc                               <  1:883910/67‑1 (MQ=255)
                                           ttACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTc                               <  1:758628/64‑1 (MQ=255)
                                             aCCGCACCAATTGCGGCACCCAACGCGAGagcagca                                                         >  1:2588609/1‑36 (MQ=255)
                                             aCCGCACCAATTGCGGCACCCAACGCGAGAGCAGCAGCGACAGc                                                 <  1:2749176/44‑1 (MQ=255)
                                                     aaTTGCGGCACCTAACGCGAGAGCAGCAGCGACAGccacc                                             >  1:4269707/1‑40 (MQ=255)
                                                     aaTTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACt                                 <  1:580920/52‑1 (MQ=255)
                                                       ttGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACt                                 <  1:201101/50‑1 (MQ=255)
                                                       ttGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATGCCGACGATAGATgc             <  1:3018083/70‑1 (MQ=255)
                                                        tGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATGCCGACGATAGATgcag           >  1:2580193/1‑71 (MQ=255)
                                                         gCGGCACCCAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATGCCGACGATAGATgcag           >  1:3859240/1‑70 (MQ=255)
                                                         gCGGCACCCAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATGCCGACGATAGATgcag           >  1:2582747/1‑70 (MQ=255)
                                                             cACCCAACGCGAGAGCAGAAGCGACAGCCACCAGCGCATTGACTTCaa                             >  1:887317/1‑48 (MQ=38)
                                                               ccTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATGCCGACGATAGATgcagca         <  1:109284/66‑1 (MQ=255)
                                                               ccTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATGCCGACGATAGATGCAGCAACTgc    <  1:1893667/71‑1 (MQ=255)
                                                                cTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTc                               <  1:4375767/43‑1 (MQ=255)
                                                                cTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTc                               <  1:3600293/43‑1 (MQ=255)
                                                                cTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTc                               <  1:1906354/43‑1 (MQ=255)
                                                                 tAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCaa                             <  1:3981423/44‑1 (MQ=255)
                                                                 tAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATGCCGACGATAGATGCAGCAACTgcgc  >  1:1034713/1‑71 (MQ=255)
                                                                 tAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATGACGACGATAGa                >  1:3631604/1‑57 (MQ=255)
                                                                 |                                                                      
GCGATAAACGCCTGGACGCGACCTTTCGCTACAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATGCCGACGATAGATGCAGCAACTGCGC  >  W3110S.gb/3701010‑3701145

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