Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,177,633 A→G 42.4% D162D (GAT→GAC arcB ← hybrid sensory histidine kinase in two‑component regulatory system with ArcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,177,6330AG42.4% 16.6 / 49.2 66D162D (GAT→GACarcBhybrid sensory histidine kinase in two‑component regulatory system with ArcA
Reads supporting (aligned to +/- strand):  ref base A (18/20);  new base G (11/17);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTACAGCCGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTTGTGTCTCT  >  minE/2177575‑2177687
                                                          |                                                      
ggTTACAGCCGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCaaggaagg                                              <  1:730423/69‑1 (MQ=255)
ggTTACAGCCGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCaaggaagg                                              <  1:65360/69‑1 (MQ=255)
ggTTACAGCCGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCaaggaagg                                              <  1:827035/69‑1 (MQ=255)
ggTTACAGCCGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCaaggaagg                                              <  1:2978874/69‑1 (MQ=255)
ggTTACAGCCGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCaaggaagg                                              <  1:2832909/69‑1 (MQ=255)
ggTTACAGCCGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCaaggaagg                                              <  1:2258311/69‑1 (MQ=255)
ggTTACAGCCGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCaaggaagg                                              <  1:1823646/69‑1 (MQ=255)
ggTTACAGCCGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCaaggaagg                                              <  1:837850/69‑1 (MQ=255)
 gTTACAGCCGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCaaggaagg                                              <  1:2489530/68‑1 (MQ=255)
    aCAGCCGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAggaaggaa                                            >  1:619428/1‑67 (MQ=255)
    aCAGCCGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAAc                                           >  1:196527/1‑68 (MQ=255)
    aCAGCCGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTAGGGTGAAGCATCAAGGAAGGAAc                                           >  1:2517176/1‑68 (MQ=255)
    aCAACCGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAggaaggaa                                            >  1:3263247/1‑67 (MQ=255)
         cGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGa                                     <  1:861806/69‑1 (MQ=255)
         cGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGa                                     <  1:2562171/69‑1 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCa                                                        >  1:1926377/1‑43 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGaa                                    >  1:3079641/1‑63 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                <  1:2455472/67‑1 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                >  1:2421247/1‑67 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                <  1:744917/67‑1 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                >  1:364704/1‑67 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                >  1:1101742/1‑67 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                >  1:2779768/1‑67 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                <  1:2802299/67‑1 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                >  1:2804375/1‑67 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                >  1:3038456/1‑67 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                >  1:3208582/1‑67 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                <  1:241127/67‑1 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                <  1:1787703/67‑1 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                >  1:1692263/1‑67 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                >  1:1924469/1‑67 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                >  1:1971489/1‑67 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                >  1:2064573/1‑67 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                >  1:1604980/1‑67 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                >  1:564110/1‑67 (MQ=255)
                ctctTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAg                                <  1:2371504/67‑1 (MQ=255)
                               aCGATAAAAAACCAGGTCGGGTGAAGCATCaaggaagg                                              <  1:2134847/38‑1 (MQ=255)
                               aCGATAAAAAACCAGGTCGGGTGAAGCATCaaggaagg                                              <  1:613991/38‑1 (MQ=255)
                               aCGATAAAAAACCAGGTCGGGTGAAGCATCaaggaagg                                              <  1:1531206/38‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:3307725/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:652663/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:3257149/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:1248099/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:738704/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:957731/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:3175429/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:828018/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:2804072/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:2734492/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:2559586/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:2338574/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:2263033/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:221281/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:2099648/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:1971056/68‑1 (MQ=255)
                                       aaaCCAGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTtgt        <  1:1657570/68‑1 (MQ=255)
                                            aGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTg               >  1:847118/1‑56 (MQ=255)
                                            aGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAAt            >  1:2267433/1‑59 (MQ=255)
                                            aGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTTGTGTctct  >  1:2687690/1‑69 (MQ=255)
                                            aGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTTGTGTctct  >  1:53458/1‑69 (MQ=255)
                                            aGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTTGTGTctct  >  1:385249/1‑69 (MQ=255)
                                            aGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTTGTGTctct  >  1:3049562/1‑69 (MQ=255)
                                            aGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTTGTGTctct  >  1:304684/1‑69 (MQ=255)
                                            aGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTTGTGTctct  >  1:1672298/1‑69 (MQ=255)
                                            aGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTTGTGTctct  >  1:1583861/1‑69 (MQ=255)
                                            aGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTTGTGTctct  >  1:1581642/1‑69 (MQ=255)
                                            aGGTCGGGTGAAGCGTCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTTGTGTctct  >  1:1215270/1‑69 (MQ=255)
                                                          |                                                      
GGTTACAGCCGGAAAACTCTTTATCTTCGTTACGATAAAAAACCAGGTCGGGTGAAGCATCAAGGAAGGAACGTAAGAATGAGGATTGCTGCTCGAGCTGAATTTGTGTCTCT  >  minE/2177575‑2177687

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