Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 423,196 A→G 100% N144N (AAT→AAC)  nagD ← UMP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE423,1960AG77.3% 44.2 / 7.0 35N144N (AAT→AAC) nagDUMP phosphatase
Reads supporting (aligned to +/- strand):  ref base A (7/1);  new base G (9/18);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 1.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGT  >  minE/423151‑423254
                                             |                                                          
caACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCTGATTGCTGGCGATAAAACGTg                                         <  1:901136/65‑1 (MQ=255)
            aGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGcgacga                          >  1:1362888/1‑68 (MQ=255)
            aGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGcgacga                          >  1:2372855/1‑68 (MQ=255)
            aGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGcgacga                          >  1:2151048/1‑68 (MQ=255)
            aGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGagacga                          >  1:1392934/1‑68 (MQ=255)
                                    ggTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGt  >  1:1405713/1‑68 (MQ=255)
                                    ggTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGt  >  1:815268/1‑68 (MQ=255)
                                    ggTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGt  >  1:2179084/1‑68 (MQ=255)
                                     gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat       >  1:3129417/4‑62 (MQ=255)
                                     gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat       >  1:3115491/4‑62 (MQ=255)
                                     gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat       >  1:2644389/4‑62 (MQ=255)
                                     gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat       >  1:2437252/4‑62 (MQ=255)
                                     gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat       >  1:864941/4‑62 (MQ=255)
                                     gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat       >  1:2050943/4‑62 (MQ=255)
                                     gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat       >  1:1615966/4‑62 (MQ=255)
                                     gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat       >  1:916785/4‑62 (MQ=255)
                                     gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat       >  1:1313863/4‑62 (MQ=255)
                                          cGGGTTGGTTGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:3157537/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:926658/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:699879/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:568845/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:3210528/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:1084625/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:998982/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:28295/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:2791785/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:2651558/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:2287086/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:2069954/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:1689854/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:1592616/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:1544053/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:109947/60‑1 (MQ=255)
                                          cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCCTTATGcatcat    <  1:3096913/60‑1 (MQ=255)
                                          cGGGTTGGT‑GCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat    <  1:3297173/59‑1 (MQ=38)
                                             |                                                          
CAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGT  >  minE/423151‑423254

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