Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 756,089 A→G 57.7% P72P (CCA→CCG lolD → outer membrane‑specific lipoprotein transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE756,0890AG57.7% 17.6 / 42.6 59P72P (CCA→CCGlolDouter membrane‑specific lipoprotein transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (11/14);  new base G (19/15);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TGCACCTGCTGGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCAAACTGTCGTCGGCAGCGAAAGCTGAACTGCGCAACCAGAAGCTGGGCTTTATTTATC  >  minE/756028‑756153
                                                             |                                                                
tgcACCTGTTAGGCTGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:2553855/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:1642026/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:2502614/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:2480684/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:2403138/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:3322109/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:46349/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:2733827/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:2128401/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:329439/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:2021065/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:1367228/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:3289808/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:1157410/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:3021656/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:3009144/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:1559039/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATCTTTAACGGTCAACCGATGAGCa                                                           >  1:1479791/1‑69 (MQ=255)
tgcACCTGTTAGGCGGGATGGATACGCCAACCTCCGGCGATGTGATCTTTAAAGGTCAACCGATGAGCa                                                           >  1:1810906/1‑69 (MQ=255)
 gcACCTGCTGGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCaa                                                          <  1:1526423/69‑1 (MQ=255)
 gcACCTGCTGGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCaa                                                          <  1:700827/69‑1 (MQ=255)
 gcACCTGCTGGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCaa                                                          <  1:1771860/69‑1 (MQ=255)
 gcACCTGCTGGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCaa                                                          <  1:924552/69‑1 (MQ=255)
 gcACCTGCTGGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCaa                                                          <  1:2672396/69‑1 (MQ=255)
 gcACCTGCTGGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCaa                                                          <  1:1614892/69‑1 (MQ=255)
 gcACCTGCTGGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCaa                                                          <  1:1591694/69‑1 (MQ=255)
 gcACCTGCTGGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCaa                                                          <  1:1089512/69‑1 (MQ=255)
    ccTGCTGGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCAAACt                                                       >  1:382598/1‑69 (MQ=255)
    ccTGCTGGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCAAACt                                                       >  1:639051/1‑69 (MQ=255)
       gctgGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCAAACTgtc                                                    >  1:1074869/1‑69 (MQ=255)
       gctgGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCAAACTgtc                                                    >  1:2604047/1‑69 (MQ=255)
       gctgGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCAAACTgtc                                                    >  1:2270415/1‑69 (MQ=255)
       gctgGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCAAACTgtc                                                    >  1:2244984/1‑69 (MQ=255)
       gctgGGCGCGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCAAACTgtc                                                    >  1:95558/1‑69 (MQ=255)
                       aCGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCAAACTGTCGTCGGCAGCGAAAGc                                     >  1:342403/1‑68 (MQ=255)
                        cGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCaa                                                          >  1:2377204/1‑46 (MQ=255)
                        cGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCaa                                                          >  1:472316/1‑46 (MQ=255)
                        cGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCaa                                                          >  1:857116/1‑46 (MQ=255)
                          ccAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCAAACTGTCGTCGGCAGCGAAAGCTGa                                  <  1:2549299/68‑1 (MQ=255)
                                                    cggTCAACCGATGAGCAAACTGTCGTCGGCGGCCAAAGCGGAACTGCGCAACCAGAAGCTGGGAttta        <  1:308108/67‑1 (MQ=255)
                                                    cggTCAACCGATGAGCAAACTGTCGTCGGCGGCCAAAGCGGAACTGCGCAACCAGAAGCTGGGAttta        <  1:997947/67‑1 (MQ=255)
                                                    cggTCAACCGATGAGCAAACTGTCGTCGGCGGCCAAAGCGGAACTGCGCAACCAGAAGCTGGGAttta        <  1:1303492/67‑1 (MQ=255)
                                                    cggTCAACCGATGAGCAAACTGTCGTCGGCGGCCAAAGCGGAACTGCGCAACCAGAAGCTGGGAttta        <  1:894667/67‑1 (MQ=255)
                                                    cggTCAACCGATGAGCAAACTGTCGTCGGCGGCCAAAGCGGAACTGCGCAACCAGAAGCTGGGAttta        <  1:1359508/67‑1 (MQ=255)
                                                    cggTCAACCGATGAGCAAACTGTCGTCGGCGGCCAAAGCGGAACTGCGCAACCAGAAGCTGGGAttta        <  1:803276/67‑1 (MQ=255)
                                                    cggTCAACCGATGAGCAAACTGTCGTCGGCGGCCAAAGCGGAACTGCGCAACCAGAAGCTGGGAttta        <  1:2065434/67‑1 (MQ=255)
                                                    cggTCAACCGATGAGCAAACTGTCGTCGGCGGCCAAAGCGGAACTGCGCAACCAGAAGCTGGGAttta        <  1:2262893/67‑1 (MQ=255)
                                                    cggTCAACCGATGAGCAAACTGTCGTCGGCGGCCAAAGCGGAACTGCGCAACCAGAAGCTGGGAttta        <  1:2727273/67‑1 (MQ=255)
                                                    cggTCAACCGATGAGCAAACTGTCGTCGGCGGCCAAAGCGGAACTGCGCAACCAGAAGCTGGGAttta        <  1:3298396/67‑1 (MQ=255)
                                                    cggTCAACCGATGAGCAAACTGTCGTCGGCGGCCAAAGCGGAACTGCGCAACCAGAAGCTGGGAttta        <  1:2743736/67‑1 (MQ=255)
                                                    cggTCAACCGATGAGCAAACTGTCGTCGGCGGCCAAAGCGGAACTGCGCAACCAGAAGCTGGGAttta        <  1:3136977/67‑1 (MQ=255)
                                                     ggTCAACCGATGAGCAAACTGTCGTCGGCGGCCAAAGCGGAACTGCGCAACCAGAAGCTGGGAttta        <  1:2263352/67‑1 (MQ=255)
                                                     ggTCAACCGATGAGCAAACTGTCGTCGGCGGCCAAAGCGGAACTGCGCAACCAGAAGCTGGGAttta        <  1:3308368/67‑1 (MQ=255)
                                                       tCAACCGATGAGCAAACTGTCGTCGGCGGCCAAAGCGGAACTGCGCAACCAGAAGCTGGGAttta        <  1:118272/65‑1 (MQ=255)
                                                         aGCCAATGAGCAAACTGTCGTCGGCAGCGAAAGCTGAACTGCGCAACCAGAAGCTGGGCTTTATTTATc  <  1:1026450/69‑1 (MQ=255)
                                                         aGCCAATGAGCAAACTGTCGTCGGCAGCGAAAGCTGAACTGCGCAACCAGAAGCTGGGCTTTATTTATc  <  1:49581/69‑1 (MQ=255)
                                                         aGCCAATGAGCAAACTGTCGTCGGCAGCGAAAGCTGAACTGCGCAACCAGAAGCTGGGCTTTATTTATc  <  1:1587271/69‑1 (MQ=255)
                                                         aGCCAATGAGCAAACTGTCGTCGGCAGCGAAAGCTGAACTGCGCAACCAGAAGCTGGGCTTTATTTATc  <  1:1524356/69‑1 (MQ=255)
                                                         aGCCAATGAGCAAACTGTCGTCGGCAGCGAAAGCTGAACTGCGCAACCAGAAGCTGGGCTTTATTTATc  <  1:929011/69‑1 (MQ=255)
                                                             |                                                                
TGCACCTGCTGGGCGGGCTGGATACGCCAACCTCCGGCGATGTGATTTTTAATGGTCAGCCAATGAGCAAACTGTCGTCGGCAGCGAAAGCTGAACTGCGCAACCAGAAGCTGGGCTTTATTTATC  >  minE/756028‑756153

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: