Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,016,291 T→C 100% V52A (GTA→GCA)  rmf → ribosome modulation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,016,2910TC97.0% 106.4 / ‑3.9 33V52A (GTA→GCA) rmfribosome modulation factor
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (8/24);  total (8/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATAC  >  W3110S.gb/1016227‑1016358
                                                                |                                                                   
tATCAGACGCTGAATCAAAGGTCACAAT‑‑‑‑‑‑‑‑GGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTg                                                         <  1:3445416/69‑1 (MQ=255)
tATCAGACGCTGAATCAAAGGTCACAAT‑‑‑‑‑‑‑‑GGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTg                                                         <  1:2489916/69‑1 (MQ=255)
tATCAGACGCTGAATCAAAGGTCACAAT‑‑‑‑‑‑‑‑GGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTg                                                         <  1:2003973/69‑1 (MQ=255)
tATCAGACGCTGAATCAAAGGTCACAAT‑‑‑‑‑‑‑‑GGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTg                                                         <  1:1960987/69‑1 (MQ=255)
 aTCAGACGCTGAATCAAAGGTCACAAT‑‑‑‑‑‑‑‑GGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTg                                                         <  1:769437/68‑1 (MQ=255)
  tCAGACGCTGAATCAAAGGTCACAAT‑‑‑‑‑‑‑‑GGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTg                                                         <  1:2332383/67‑1 (MQ=255)
    aGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGcc                                                           >  1:2465533/1‑71 (MQ=255)
             aTCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACATGGCAGTAATGGCCTGATTCTGt                                                  <  1:3716273/71‑1 (MQ=255)
             aTCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGt                                                  <  1:1082603/71‑1 (MQ=255)
             aTCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGt                                                  <  1:898163/71‑1 (MQ=255)
             aTCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGt                                                  <  1:1915268/71‑1 (MQ=255)
             aTCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGt                                                  <  1:1318895/71‑1 (MQ=255)
                      cacaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGgcag                                                                   >  1:809793/1‑45 (MQ=255)
                      cacaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTaa                                          >  1:1911775/1‑70 (MQ=255)
                      cacaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTaa                                          >  1:1866678/1‑70 (MQ=255)
                          aTGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTTTCTCTTTAAAAAGa                                     <  1:2109506/71‑1 (MQ=255)
                           tGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGaa                                    <  1:2177511/71‑1 (MQ=255)
                           tGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGaa                                    <  1:1408985/71‑1 (MQ=255)
                           tGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGaa                                    <  1:2749810/71‑1 (MQ=255)
                              cTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGgcagt                                                                  <  1:2654643/38‑1 (MQ=255)
                                  gAGGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCg                             <  1:300169/71‑1 (MQ=255)
                                     gCTGGCGAGAAGGCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGa                                     <  1:2319999/60‑1 (MQ=255)
                                     gCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGa                                     <  1:2015452/60‑1 (MQ=255)
                                     gCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGa                                     <  1:3367439/60‑1 (MQ=255)
                                     gCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGa                                     <  1:421585/60‑1 (MQ=255)
                                     gCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGa                                     <  1:1818717/60‑1 (MQ=255)
                                         gCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTg                        <  1:908745/69‑1 (MQ=255)
                                                   aTGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGaa                                    <  1:1141936/47‑1 (MQ=255)
                                                      gCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGa                                     >  1:3305689/1‑43 (MQ=255)
                                                      gCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGa                                     >  1:3087077/1‑43 (MQ=255)
                                                      gCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGa                                     >  1:3079557/1‑43 (MQ=255)
                                                      gCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGa                                     >  1:21432/1‑43 (MQ=255)
                                                      gCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGa                                     >  1:1330820/1‑43 (MQ=255)
                                                          aCAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGa     <  1:972096/71‑1 (MQ=255)
                                                             gggCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATAc  <  1:2754299/71‑1 (MQ=255)
                                                                |                                                                   
TATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATAC  >  W3110S.gb/1016227‑1016358

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: