Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,419,077 G→A 33.7% L296L (CTA→TTA)  ccmH ← heme lyase, CcmH subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,419,0770GA33.7% 53.0 / 34.9 68L296L (CTA→TTA) ccmHheme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base G (19/26);  new base A (11/12);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTG  >  minE/1419010‑1419141
                                                                   |                                                                
atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                                  >  1:1076911/1‑68 (MQ=255)
atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                                  >  1:1110800/1‑68 (MQ=255)
atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                                  >  1:1007424/1‑68 (MQ=255)
atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                                  >  1:2691490/1‑68 (MQ=255)
atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                                  >  1:894459/1‑68 (MQ=255)
atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                                  >  1:1627222/1‑68 (MQ=255)
atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                                  >  1:763046/1‑68 (MQ=255)
atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                                  >  1:1581006/1‑68 (MQ=255)
atcTCCAACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                                  >  1:437335/1‑68 (MQ=255)
     ccACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACg                                                            >  1:89383/1‑69 (MQ=255)
     acACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACg                                                            >  1:1861820/2‑69 (MQ=255)
       aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCTCACGGa                                                          <  1:2825431/69‑1 (MQ=255)
       aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGa                                                          <  1:829395/69‑1 (MQ=255)
       aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGa                                                          <  1:2925769/69‑1 (MQ=255)
       aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGa                                                          <  1:3079513/69‑1 (MQ=255)
       aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGa                                                          <  1:2408556/69‑1 (MQ=255)
       aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGa                                                          <  1:293472/69‑1 (MQ=255)
       aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGa                                                          <  1:2796787/69‑1 (MQ=255)
       aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACAAGGa                                                          <  1:3421798/69‑1 (MQ=255)
          cgGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGatat                                                       <  1:2292541/69‑1 (MQ=255)
          cgGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGatat                                                       <  1:2128819/69‑1 (MQ=255)
          cgGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGatat                                                       <  1:1480012/69‑1 (MQ=255)
                cGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTa                                                  <  1:2789326/68‑1 (MQ=255)
                        ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                                         >  1:2247670/1‑69 (MQ=255)
                        ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                                         >  1:2192357/1‑69 (MQ=255)
                        ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                                         >  1:2483798/1‑69 (MQ=255)
                        ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                                         >  1:214881/1‑69 (MQ=255)
                        ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                                         >  1:2906538/1‑69 (MQ=255)
                        ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                                         >  1:2992264/1‑69 (MQ=255)
                        ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                                         >  1:3058070/1‑69 (MQ=255)
                        ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                                         >  1:3082894/1‑69 (MQ=255)
                        ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                                         >  1:1314563/1‑69 (MQ=255)
                        ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                                         >  1:1181174/1‑69 (MQ=255)
                        acAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGt                                          >  1:3058074/2‑68 (MQ=255)
                                     gctCAAAGGCATTAAACGCATACATGCTTAGCACACGGa                                                          <  1:810974/39‑1 (MQ=255)
                                           aGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGc                       <  1:2565719/68‑1 (MQ=255)
                                           aGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGc                       <  1:2225789/68‑1 (MQ=255)
                                               aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAGCTGGCGCAGCAGTTc                  <  1:1596765/69‑1 (MQ=255)
                                               aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc                  <  1:979732/69‑1 (MQ=255)
                                               aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc                  <  1:174010/69‑1 (MQ=255)
                                               aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc                  <  1:1809885/69‑1 (MQ=255)
                                               aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc                  <  1:3035326/69‑1 (MQ=255)
                                               aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc                  <  1:1663608/69‑1 (MQ=255)
                                               aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc                  <  1:1992957/69‑1 (MQ=255)
                                               aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc                  <  1:1649500/69‑1 (MQ=255)
                                               aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc                  <  1:512668/69‑1 (MQ=255)
                                               aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc                  <  1:2186001/69‑1 (MQ=255)
                                               aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc                  <  1:523255/69‑1 (MQ=255)
                                               aTTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTc                  <  1:2647015/69‑1 (MQ=255)
                                               aTTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTc                  <  1:1748453/69‑1 (MQ=255)
                                               aTTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTc                  <  1:3386092/69‑1 (MQ=255)
                                                ttAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc                  <  1:1947206/68‑1 (MQ=255)
                                                           cATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGc                       <  1:434157/52‑1 (MQ=255)
                                                           cATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGc                       <  1:1555349/52‑1 (MQ=255)
                                                           cATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGc                       <  1:1447052/52‑1 (MQ=255)
                                                           cATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGc                       <  1:314654/52‑1 (MQ=255)
                                                           cATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGc                       <  1:1669495/52‑1 (MQ=255)
                                                           cATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGc                       <  1:1848111/52‑1 (MQ=255)
                                                           cATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGc                       <  1:2473216/52‑1 (MQ=255)
                                                           cATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGc                       <  1:2524498/52‑1 (MQ=255)
                                                           cATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCg      >  1:316850/1‑69 (MQ=255)
                                                           cATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCg      >  1:1411843/1‑69 (MQ=255)
                                                               cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTg  >  1:499543/1‑69 (MQ=255)
                                                               cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTg  >  1:2548540/1‑69 (MQ=255)
                                                               cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTg  >  1:94097/1‑69 (MQ=255)
                                                               cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTg  >  1:141326/1‑69 (MQ=255)
                                                               cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTg  >  1:867780/1‑69 (MQ=255)
                                                               cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTg  >  1:2413882/1‑69 (MQ=255)
                                                               cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTg  >  1:2355667/1‑69 (MQ=255)
                                                               cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTg  >  1:2894764/1‑69 (MQ=255)
                                                               cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTg  >  1:2018350/1‑69 (MQ=255)
                                                               cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTg  >  1:297154/1‑69 (MQ=255)
                                                               cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTg  >  1:407612/1‑69 (MQ=255)
                                                               cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTg  >  1:396513/1‑69 (MQ=255)
                                                               cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTg  >  1:3385750/1‑69 (MQ=255)
                                                               cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTg  >  1:3149463/1‑69 (MQ=255)
                                                               cTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGAAGTTCACCGCCGAGGCGGTTg  >  1:2578731/1‑69 (MQ=255)
                                                                   |                                                                
ATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTG  >  minE/1419010‑1419141

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: ‑