Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,606,282 T→G 57.6% intergenic (+3/‑1) amiA → / → hemF N‑acetylmuramoyl‑l‑alanine amidase I/coproporphyrinogen III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,606,2820TG57.6% 35.7 / 57.8 66intergenic (+3/‑1)amiA/hemFN‑acetylmuramoyl‑l‑alanine amidase I/coproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base T (2/26);  new base G (22/16);  total (24/42)
Fisher's exact test for biased strand distribution p-value = 2.01e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

GGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCT  >  minE/1606220‑1606350
                                                              |                                                                    
ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                                <  1:1881221/69‑1 (MQ=255)
ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                                <  1:2051070/69‑1 (MQ=255)
ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                                <  1:2134287/69‑1 (MQ=255)
ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                                <  1:1748171/69‑1 (MQ=255)
ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                                <  1:2135826/69‑1 (MQ=255)
ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                                <  1:1533701/69‑1 (MQ=255)
ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                                <  1:2358139/69‑1 (MQ=255)
  cGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa                                                                <  1:1301254/67‑1 (MQ=255)
   gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                             <  1:1104405/69‑1 (MQ=255)
   gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                             <  1:2373057/69‑1 (MQ=255)
   gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                             <  1:2807089/69‑1 (MQ=255)
   gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                             <  1:1072763/69‑1 (MQ=255)
    tGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                             <  1:2874272/68‑1 (MQ=255)
         aGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGca                                                       <  1:1524868/69‑1 (MQ=255)
         aGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGca                                                       <  1:2061367/69‑1 (MQ=255)
           ttATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGca                                                       <  1:494087/67‑1 (MQ=255)
            tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACg                                                         <  1:673727/64‑1 (MQ=255)
            tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACg                                                         <  1:2049228/64‑1 (MQ=255)
            tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACg                                                         <  1:1976084/64‑1 (MQ=255)
            tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCacac                                                    <  1:2186453/69‑1 (MQ=255)
            tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCacac                                                    <  1:7288/69‑1 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGtgatga                                                                  >  1:965430/1‑46 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAGACCCGACGCACACCAGGTTaa                                            >  1:1631841/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGcaca                                                     >  1:1276619/1‑59 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:249986/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:3328248/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:3350073/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:3365461/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:3428131/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:1989972/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:283026/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:2815214/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:2734313/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:2703380/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:2401693/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:666148/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:1529633/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:1175700/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:1765700/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:2272388/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:1611271/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:2203143/1‑68 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTa                                             >  1:2378217/1‑67 (MQ=255)
                     tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGCTAAGTGATGAAACCCGACGCACACCAGGTTaa                                            >  1:878926/1‑68 (MQ=255)
                      ggTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAAc                                          <  1:2904586/69‑1 (MQ=255)
                          cGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCg                                                            >  1:2788395/1‑47 (MQ=255)
                             cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:364641/68‑1 (MQ=255)
                             cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:124013/68‑1 (MQ=255)
                             cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:1178808/68‑1 (MQ=255)
                             cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:3433238/68‑1 (MQ=255)
                             cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:784282/68‑1 (MQ=255)
                             cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:612957/68‑1 (MQ=255)
                             cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:932091/68‑1 (MQ=255)
                             cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:1315863/68‑1 (MQ=255)
                             cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:1951301/68‑1 (MQ=255)
                             cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:3000247/68‑1 (MQ=255)
                             cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:135215/68‑1 (MQ=255)
                             cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:1844939/68‑1 (MQ=255)
                             cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:2521342/68‑1 (MQ=255)
                             cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:1764234/68‑1 (MQ=255)
                             cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:1413607/68‑1 (MQ=255)
                              aaCCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                    <  1:2040196/67‑1 (MQ=255)
                                                  aaaGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAt              >  1:2735172/1‑69 (MQ=255)
                                                           aGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagca     <  1:633066/69‑1 (MQ=255)
                                                           aGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagca     <  1:901962/69‑1 (MQ=255)
                                                           aGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagca     <  1:885422/69‑1 (MQ=255)
                                                           aGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagca     <  1:3422448/69‑1 (MQ=255)
                                                              tATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt  <  1:2264820/69‑1 (MQ=255)
                                                              tATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt  <  1:2167927/69‑1 (MQ=255)
                                                              tATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt  <  1:3012266/69‑1 (MQ=255)
                                                              tATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt  <  1:535753/69‑1 (MQ=255)
                                                              tATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt  <  1:420518/69‑1 (MQ=255)
                                                              tATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt  <  1:2831933/69‑1 (MQ=255)
                                                              tATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt  <  1:3001603/69‑1 (MQ=255)
                                                              |                                                                    
GGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCT  >  minE/1606220‑1606350

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: