Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,891,664 Δ1 bp 100% intergenic (‑82/+409) gudP ← / ← ygdL predicted D‑glucarate transporter/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,891,6640T.79.2% 55.9 / 16.8 24intergenic (‑82/+409)gudP/ygdLpredicted D‑glucarate transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base . (19/0);  total (21/3)
Fisher's exact test for biased strand distribution p-value = 4.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CTGGACTGCGACATCGCCAGGAAATCAGAGGTGACGTAGGGTGTTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTATATCTGGCAATAGCAGTATA  >  minE/1891621‑1891730
                                           |                                                                  
ctGGACTGCGACATCGCCAGGAAATCAGAGGTGACGTAGGGTGTTTTTTGCCGTTTTTATAGGTCGTTc                                           >  1:3083601/1‑69 (MQ=255)
       gCGACATCGCCAGGAAATCAGAGGTGACGTAGGGTGTTTTTTGCCGTTTTTATAGGTCGTTCGCCGAAt                                    >  1:2057321/1‑69 (MQ=255)
                   ggAAATCAGAGGTGACGTAGGGTGTTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGtt                        <  1:1693847/69‑1 (MQ=255)
                   ggAAATCAGAGGTGACGTAGGGTGTTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGtt                        <  1:2220755/69‑1 (MQ=255)
                      aaTCAGAGGTGACGTAGGGTGTTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTtat                     <  1:450582/69‑1 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCgt                           >  1:1664836/1‑61 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTtata                    >  1:1633722/1‑68 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:2767454/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:976882/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:771952/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:735436/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:700531/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:525547/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:3120676/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:3084036/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:3012862/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:1165277/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:2679102/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:2393358/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:2098551/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:1664933/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:1583148/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTatat                   >  1:1185871/1‑69 (MQ=255)
                       aTCAGAGGTGACGTAGGGTGTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTCTTTatat                   >  1:763219/1‑69 (MQ=255)
                                         tgtTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTATATCTGGCAATAGCAGtat   >  1:2866848/1‑68 (MQ=255)
                                          gtTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTATATCTGGCAATAGCAGtata  >  1:1457259/1‑68 (MQ=255)
                                           |                                                                  
CTGGACTGCGACATCGCCAGGAAATCAGAGGTGACGTAGGGTGTTTTTTGCCGTTTTTATAGGTCGTTCGCCGAATACGGCGCGTGTTTATATCTGGCAATAGCAGTATA  >  minE/1891621‑1891730

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: