Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 815,460 G→C 68.8% L85L (CTC→CTG lolB ← chaperone for lipoproteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE815,4600GC68.8% 52.1 / 30.8 48L85L (CTC→CTGlolBchaperone for lipoproteins
Reads supporting (aligned to +/- strand):  ref base G (13/2);  new base C (13/20);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 4.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

GCACGTTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTT  >  minE/815405‑815524
                                                       |                                                                
gCACGTTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGt                                                     >  1:1265410/1‑69 (MQ=255)
gCACGTTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGt                                                     >  1:1016237/1‑69 (MQ=255)
gCACGTTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGt                                                     >  1:1803770/1‑69 (MQ=255)
     ttACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGa                                                <  1:2289519/69‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCGGCAGACGGt                                                     <  1:1715640/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACgtt                                                     <  1:1416633/59‑3 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:3390457/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:338517/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:3318612/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:889591/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:319090/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:3019269/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:3007592/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:2979524/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:2831615/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:2689561/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:2245080/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:2176474/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:1700571/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:1739841/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:2236271/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:2109841/59‑1 (MQ=255)
          cGGTTGAGCATTCAGCTCCAGTCCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:1674755/59‑1 (MQ=255)
          cGGTTGAGCATTAAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTCAGCAGCAGACGGt                                                     <  1:3270408/59‑1 (MQ=255)
               gAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGAt                                               >  1:1812197/1‑60 (MQ=255)
               gAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGAt                                               >  1:2941799/1‑60 (MQ=255)
               gAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGATAGCGAt                                               >  1:825712/1‑60 (MQ=255)
                      cAGCTCAAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTtgct                                >  1:3041605/1‑68 (MQ=255)
                           ccAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAg                           >  1:818682/1‑68 (MQ=255)
                           ccAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAg                           >  1:2773564/1‑68 (MQ=255)
                           ccAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAg                           >  1:692747/1‑68 (MQ=255)
                           ccAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAg                           >  1:3420893/1‑68 (MQ=255)
                              gTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGaaaa                       <  1:2688819/69‑1 (MQ=255)
                                      ctgcCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCg               >  1:1315182/1‑69 (MQ=255)
                                               gggTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACtttt      >  1:2112430/1‑69 (MQ=255)
                                                   tAGTCAGCAGCAGACGGTAGCGATCCTGGCCAGTTTGCTGCCAGAAGAAACGGGCGTACACTTTttgtt  >  1:2347956/1‑69 (MQ=255)
                                                   tAGTCAGCAGCAGACGGTAGCGATCCTGGCCAGTTTGCTGCCAGAAGAAACGGGCGTACACTTTttgtt  >  1:954651/1‑69 (MQ=255)
                                                   tAGTCAGCAGCAGACGGTAGCGATCCTGGCCAGTTTGCTGCCAGAAGAAACGGGCGTACACTTTttgtt  >  1:919106/1‑69 (MQ=255)
                                                   tAGTCAGCAGCAGACGGTAGCGATCCTGGCCAGTTTGCTGCCAGAAGAAACGGGCGTACACTTTttgtt  >  1:2310204/1‑69 (MQ=255)
                                                   tAGTCAGCAGCAGACGGTAGCGATCCTGGCCAGTTTGCTGCCAGAAGAAACGGGCGTACACTTTttgtt  >  1:456991/1‑69 (MQ=255)
                                                   tAGTCAGCAGCAGACGGTAGCGATCCTGGCCAGTTTGCTGCCAGAAGAAACGGGCGTACACTTTttgtt  >  1:3391411/1‑69 (MQ=255)
                                                   tAGTCAGCAGCAGACGGTAGCGATCCTGGCCAGTTTGCTGCCAGAAGAAACGGGCGTACACTTTttgtt  >  1:292410/1‑69 (MQ=255)
                                                   tAGTCAGCAGCAGACGGTAGCGATCCTGGCCAGTTTGCTGCCAGAAGAAACGGGCGTACACTTTttgtt  >  1:2134872/1‑69 (MQ=255)
                                                   tAGTCAGCAGCAGACGGTAGCGATCCTGGCCAGTTTGCTGCCAGAAGAAACGGGCGTACACTTTttgtt  >  1:3120113/1‑69 (MQ=255)
                                                   tAGTCAGCAGCAGACGGTAGCGATCCTGGCCAGTTTGCTGCCAGAAGAAACGGGCGTACACTTTttgtt  >  1:2859625/1‑69 (MQ=255)
                                                   tAGTCAGCAGCAGACGGTAGCGATCCTGGCCAGTTTGCTGCCAGAAGAAACGGGCGTACACTTTttgtt  >  1:1925989/1‑69 (MQ=255)
                                                   tAGTCAGCAGCAGACGGTAGCGATCCTGGCCAGTTTGCTGCCAGAAGAAACGGGCGTACACTTTttgtt  >  1:2910288/1‑69 (MQ=255)
                                                   tAGTCAGCAGCAGACGGTAGCGATCCTGGCCAGTTTGCTGCCAGAAGAAACGGGCGTACACTTTttgt   >  1:1926437/1‑68 (MQ=255)
                                                       |                                                                
GCACGTTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTT  >  minE/815405‑815524

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