Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,695,660 C→A 20.1% R117L (CGC→CTC)  yfhQ ← predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,695,6600CA20.1% 198.0 / 19.6 64R117L (CGC→CTC) yfhQpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (36/15);  new base A (9/4);  total (45/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCTTTTCAGGCCGCATTCGC  >  minE/1695590‑1695729
                                                                      |                                                                     
gTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGAcgc                                                                       >  1:373955/1‑71 (MQ=255)
   ttCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCt                                                                    >  1:212631/1‑71 (MQ=255)
   ttCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCt                                                                    >  1:593691/1‑71 (MQ=255)
   ttCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCt                                                                    >  1:2717506/1‑71 (MQ=255)
        ggTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGAGCTCGCGa                                                               >  1:1214368/1‑71 (MQ=255)
            aGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCa                                                            <  1:482930/70‑1 (MQ=255)
            aGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCa                                                            <  1:449104/70‑1 (MQ=255)
            aGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCa                                                            >  1:1655328/1‑70 (MQ=255)
              ccGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAac                                                         <  1:3126556/71‑1 (MQ=255)
              ccGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAac                                                         <  1:2410907/71‑1 (MQ=255)
                  aaTCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGAGCTCGCGACCAAacac                                                       >  1:1760758/1‑69 (MQ=255)
                  aaTAGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAac                                                         >  1:1656881/1‑67 (MQ=255)
                   aTCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAg                                                    >  1:417782/1‑71 (MQ=255)
                   aTCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAg                                                    <  1:1735357/71‑1 (MQ=255)
                      gcgACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGa                                                               >  1:744710/1‑57 (MQ=255)
                      gcgACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAaca                                                        >  1:169989/1‑64 (MQ=255)
                      gcgACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAaca                                                        >  1:1952189/1‑64 (MQ=255)
                      gcgACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGAGCTCGCGACCaaa                                                          >  1:2050216/1‑62 (MQ=255)
                           aTGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGAGCTcgc                                                                 <  1:1400477/50‑1 (MQ=255)
                           aTGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGAGCTcgc                                                                 <  1:2903716/50‑1 (MQ=255)
                           aTGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGAGCTcgc                                                                 <  1:2276534/50‑1 (MQ=255)
                           aTGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGAGCTcgc                                                                 <  1:2276275/50‑1 (MQ=255)
                            tGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCgg                                           >  1:2781341/1‑71 (MQ=255)
                            tGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCgg                                           >  1:535829/1‑71 (MQ=255)
                            agATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCgg                                           >  1:1567463/2‑71 (MQ=255)
                             gATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCggg                                          >  1:976124/1‑71 (MQ=255)
                             gATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCggg                                          >  1:1039168/1‑71 (MQ=255)
                             gATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGAGCTCGCGa                                                               >  1:610229/1‑50 (MQ=255)
                               tAATGGAATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCgggg                                         >  1:4483/1‑70 (MQ=255)
                                 aTGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCaaa                                                          >  1:727016/1‑51 (MQ=255)
                                 aTGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCaaa                                                          >  1:120518/1‑51 (MQ=255)
                                   ggCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCaaa                                                          >  1:210856/1‑49 (MQ=255)
                                   ggCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCaaa                                                          >  1:969232/1‑49 (MQ=255)
                                   ggCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCaaa                                                          >  1:1308527/1‑49 (MQ=255)
                                   ggCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTAttt                                    >  1:1714670/1‑71 (MQ=255)
                                      aTTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGcc                                 >  1:2761636/1‑71 (MQ=255)
                                       tttCTGCAACTCTTCATTGGTCAAGCCGACGCGCTcgc                                                                 <  1:1231023/38‑1 (MQ=255)
                                       tttCTGCAACTCTTCATTGGTCAAGCCGACGCGCTcgc                                                                 <  1:487869/38‑1 (MQ=255)
                                         tctcCAACTCTTCATTGGTCAAGCCGACGAGCTCGCGACCAAACACCAg                                                    >  1:735747/1‑49 (MQ=38)
                                         tCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGa                                                               >  1:2170875/1‑38 (MQ=255)
                                         tCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGAc                                                              >  1:2087690/1‑39 (MQ=255)
                                         tCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGAc                                                              >  1:1088859/1‑39 (MQ=255)
                                         tCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCt                              >  1:782215/1‑71 (MQ=255)
                                         tCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCt                              >  1:3198238/1‑71 (MQ=255)
                                           tGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCa                                                            <  1:2618067/39‑1 (MQ=255)
                                           tGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCa                                                            <  1:2614086/39‑1 (MQ=255)
                                           tGCAACTCTTCATTGGTCAAGCCGACGAGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCtt                             >  1:815283/1‑70 (MQ=255)
                                               aCTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGc                         <  1:10592/70‑1 (MQ=255)
                                                ctctTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGa                       >  1:713897/1‑71 (MQ=255)
                                                    tCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCt                   >  1:64808/1‑71 (MQ=255)
                                                     cATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGc                    <  1:1578125/69‑1 (MQ=255)
                                                      aTTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCtt                  <  1:2135139/70‑1 (MQ=255)
                                                       ttGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCtttt                <  1:1202306/71‑1 (MQ=255)
                                                       ttGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCtttt                <  1:889647/71‑1 (MQ=255)
                                                       ttGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCtttt                <  1:1076948/71‑1 (MQ=255)
                                                         ggTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCtttt                >  1:2341105/1‑69 (MQ=255)
                                                           tCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCTTTTCAgg            >  1:211700/1‑71 (MQ=255)
                                                           tCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCTTTTCAgg            >  1:2822608/1‑71 (MQ=255)
                                                            cAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCTTTTCAGGc           >  1:3160091/1‑71 (MQ=255)
                                                             aaGCCGACGAGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCa                           >  1:811498/1‑54 (MQ=255)
                                                             aaGCCGACGAGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCa                           >  1:1001368/1‑54 (MQ=255)
                                                             aaGCCGACGAGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCa                           >  1:2169880/1‑54 (MQ=255)
                                                              aGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCTTTTCAGGCCg         >  1:300603/1‑71 (MQ=255)
                                                               gCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCTTTTCAGGCCg         >  1:1566040/1‑70 (MQ=255)
                                                                ccGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCTTTTCAGGCCGCa       >  1:525311/1‑71 (MQ=255)
                                                                     gcgcTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCTTTTCAGGCCGCATTcgc  >  1:1986052/1‑71 (MQ=255)
                                                                     gcgcTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCTTTTCAGGCCGCATTcgc  >  1:434225/1‑71 (MQ=255)
                                                                     gcgcTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCTTTTCAGGCCGCATTcgc  >  1:3199907/1‑71 (MQ=255)
                                                                     gcgcTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCTTTTCAGGCCGCATTcgc  >  1:1334698/1‑71 (MQ=255)
                                                                      |                                                                     
GTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACGCTTTTCAGGCCGCATTCGC  >  minE/1695590‑1695729

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: