Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,012,376 G→A 56.5% A500T (GCC→ACC)  rbbA → fused ribosome‑associated ATPases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,012,3760GA56.5% 2.9 / 26.9 23A500T (GCC→ACC) rbbAfused ribosome‑associated ATPases
Reads supporting (aligned to +/- strand):  ref base G (8/2);  new base A (9/4);  total (17/6)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

TCACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAA  >  W3110S.gb/4012310‑4012434
                                                                  |                                                          
tcACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAaccgc                                                        >  1:890898/1‑71 (MQ=255)
        gCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTg                                                 >  1:894102/1‑70 (MQ=255)
         cACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCt                                                  >  1:340224/1‑68 (MQ=255)
         cACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTgg                                                >  1:84071/1‑70 (MQ=255)
            gCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGaa                                              >  1:1277372/1‑69 (MQ=255)
              cGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAgg                                          <  1:471703/71‑1 (MQ=255)
                 aaaaGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCa                                        >  1:620850/1‑70 (MQ=255)
                 aaaaGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCAt                                       >  1:169886/1‑71 (MQ=255)
                  aaaGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCAtt                                      >  1:203961/1‑71 (MQ=255)
                   aaGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCAttt                                     <  1:467614/71‑1 (MQ=255)
                       gCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCg                                 >  1:223178/1‑71 (MQ=255)
                       gCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCg                                 >  1:682673/1‑71 (MQ=255)
                           gCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTa                             >  1:1304627/1‑71 (MQ=255)
                               gCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTg                         >  1:1555503/1‑71 (MQ=255)
                                         gCAGGAACTGGTTGAGAAACGCGTAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCgg               >  1:1410387/1‑71 (MQ=255)
                                         gCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCAttt                                     <  1:1572104/49‑1 (MQ=255)
                                          cAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCAttt                                     <  1:829111/48‑1 (MQ=255)
                                                cTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcagggca        <  1:654684/71‑1 (MQ=255)
                                                  ggatGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcagggca        >  1:1109834/4‑69 (MQ=255)
                                                  ggTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcagggca        >  1:1218690/1‑69 (MQ=255)
                                                  ggTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcagggca        >  1:1211537/1‑69 (MQ=255)
                                                   gTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCagag     >  1:1656875/1‑71 (MQ=255)
                                                     tGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCgg               <  1:751822/59‑1 (MQ=255)
                                                      gagaAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCaa  >  1:1692417/1‑71 (MQ=255)
                                                                  |                                                          
TCACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAA  >  W3110S.gb/4012310‑4012434

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: