Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,016,291 T→C 50.0% 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,2910TC50.0% ‑5.7 / 16.9 14V52A (GTA→GCA) rmfribosome modulation factor
Reads supporting (aligned to +/- strand):  ref base T (5/2);  new base C (2/5);  total (7/7)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

TATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTT  >  W3110S.gb/1016227‑1016353
                                                                |                                                              
tATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGgtagt                                                             <  1:1098309/68‑1 (MQ=255)
tATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGgtagt                                                             <  1:1276143/68‑1 (MQ=255)
tATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGgtagt                                                             <  1:2539796/68‑1 (MQ=255)
tATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGgcagta                                                            <  1:105964/69‑1 (MQ=255)
tATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGgcagta                                                            <  1:2472717/69‑1 (MQ=255)
    aGACTCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGgcagta                                                            <  1:1154915/65‑1 (MQ=255)
           gAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGAtt                                                 >  1:2454269/1‑69 (MQ=255)
           gAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGAt                                                  >  1:1295281/1‑68 (MQ=255)
           gAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGAt                                                  >  1:516370/1‑68 (MQ=255)
           gAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGAt                                                  >  1:1718197/1‑68 (MQ=255)
               cAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGt                                             <  1:13568/69‑1 (MQ=255)
                  aGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGtctc                                          >  1:3037202/1‑69 (MQ=255)
                   ggTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGgtagta                                                            >  1:3384069/1‑50 (MQ=255)
                   ggTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGgtagt                                                             >  1:1959733/1‑49 (MQ=255)
                       acaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACATGgcagt                                                             <  1:1799421/45‑1 (MQ=255)
                       acaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGgcagt                                                             <  1:3664548/45‑1 (MQ=255)
                             gCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGaa                               <  1:1584493/69‑1 (MQ=255)
                                gggAGGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAACGGCCTGATTCTGTCTCTTTAAAAAGAAAcc                            <  1:781769/69‑1 (MQ=255)
                                   aGGCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTc                          >  1:623419/1‑68 (MQ=255)
                                    ggCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCg                        <  1:1047858/69‑1 (MQ=255)
                                    ggCTGGCGAGAAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCg                        <  1:1474910/69‑1 (MQ=255)
                                            agaAGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCgg                >  1:1943402/1‑69 (MQ=255)
                                               aGCCATGGCGGACAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAgg             <  1:2163356/69‑1 (MQ=255)
                                                          aCAGGGCAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCtttt  <  1:3053549/69‑1 (MQ=255)
                                                                |                                                              
TATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTT  >  W3110S.gb/1016227‑1016353

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: