Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 381,851 A→C 57.8% G88G (GGT→GGG yaiF ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb381,8510AC57.8% 52.3 / 40.7 45G88G (GGT→GGGyaiFhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (7/12);  new base C (11/15);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGT  >  W3110S.gb/381786‑381913
                                                                 |                                                              
aGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGccc                                                              <  1:1664958/68‑1 (MQ=255)
     tCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAAt                                                        <  1:1114483/69‑1 (MQ=255)
           aaTATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACgg                                                   <  1:2546276/68‑1 (MQ=255)
                  gggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACgg                                                   <  1:2528409/61‑1 (MQ=255)
                  gggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACgg                                                   <  1:719755/61‑1 (MQ=255)
                      aaaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACaaa                                       <  1:1188890/69‑1 (MQ=255)
                      aaaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACaaa                                       <  1:3464454/69‑1 (MQ=255)
                       aaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCAt                                             <  1:1297183/62‑1 (MQ=255)
                       aaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCAt                                             <  1:1623520/62‑1 (MQ=255)
                        aTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAc                                                <  1:3152681/58‑1 (MQ=255)
                        aTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAc                                                <  1:2144648/58‑1 (MQ=255)
                        aTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAc                                                <  1:2727807/58‑1 (MQ=255)
                         tCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACa                                         >  1:1921326/1‑64 (MQ=255)
                         tCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATAc                                    >  1:3280129/1‑69 (MQ=255)
                          ccAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATaa                                                      <  1:3133366/50‑1 (MQ=255)
                          ccAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATaa                                                      <  1:28877/50‑1 (MQ=255)
                               tACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAAt                                                        <  1:2906754/43‑1 (MQ=255)
                                         gCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGAt                    <  1:2466641/69‑1 (MQ=255)
                                              ttGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaacaa               >  1:726669/1‑69 (MQ=255)
                                              ttGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaacaa               >  1:2262473/1‑69 (MQ=255)
                                              ttGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaaca                >  1:1164230/1‑68 (MQ=255)
                                              ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaaa                                       >  1:2150325/1‑45 (MQ=255)
                                              ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaaa                                       <  1:3651870/45‑1 (MQ=255)
                                              ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaaa                                       >  1:2873664/1‑45 (MQ=255)
                                              ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaaa                                       <  1:2092138/45‑1 (MQ=255)
                                              ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaaa                                       <  1:1365182/45‑1 (MQ=255)
                                              ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaaa                                       <  1:1382826/45‑1 (MQ=255)
                                              ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaaa                                       >  1:3035310/1‑45 (MQ=255)
                                              ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaaa                                       >  1:3188342/1‑45 (MQ=255)
                                              ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaaa                                       <  1:2269334/45‑1 (MQ=255)
                                               tGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACa                                         >  1:1300003/1‑42 (MQ=255)
                                               tGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACa                                         <  1:3556008/42‑1 (MQ=255)
                                               tGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACggg                                 <  1:3137654/50‑1 (MQ=255)
                                                   ccAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCa          <  1:3089969/69‑1 (MQ=255)
                                                   ccAATCACTACCGGACCCTGAATAACGTCACCATGACAAATACGTGTATTCGCGCCGATAACAACTTc           <  1:2025822/68‑1 (MQ=255)
                                                   ccAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTc           <  1:3028777/68‑1 (MQ=255)
                                                         aCTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCa    >  1:3147983/1‑69 (MQ=255)
                                                         aCTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCa    >  1:2082669/1‑69 (MQ=255)
                                                         aCTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCa    >  1:1352096/1‑69 (MQ=255)
                                                          cTACCGGCCCCTGAATAACGGCACCATGACAAATACggg                                 >  1:2097167/1‑39 (MQ=255)
                                                          cTACCGGCCCCTGAATAACGGCACCATGACAAATACggg                                 >  1:2126460/1‑39 (MQ=255)
                                                          cTACCGGCCCCTGAATAACGGCACCATGACAAATACggg                                 >  1:803094/1‑39 (MQ=255)
                                                          cTACCGGCCCCTGAATAACGGCACCATGACAAATACgg                                  >  1:3652081/1‑38 (MQ=255)
                                                            aCCGGACACTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGt  >  1:1716342/1‑68 (MQ=255)
                                                             ccGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaacaa               <  1:2598610/54‑1 (MQ=255)
                                                             ccGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaacaa               <  1:48517/54‑1 (MQ=255)
                                                                 |                                                              
AGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGT  >  W3110S.gb/381786‑381913

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: