Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,461,154 T→G 100% pseudogene (801/4605 nt)
pseudogene (801/4104 nt)
yfaS ←
yfaS ←
ECK2220:JW2222+JW2221:b4500; hypothetical protein
ECK2220:JW2222:b2228; hypothetical protein, N‑ter fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,461,1540TG83.4% 73.0 / 6.0 36pseudogene (801/4605 nt)
pseudogene (801/4104 nt)
yfaS
yfaS
ECK2220:JW2222+JW2221:b4500; hypothetical protein
ECK2220:JW2222:b2228; hypothetical protein, N‑ter fragment
Reads supporting (aligned to +/- strand):  ref base T (6/0);  new base G (21/9);  total (27/9)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCACAAGATCTCTGAGCCGGGCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTC  >  minE/1461088‑1461208
                                                                  |                                                      
ccACAAGATCTCTGAGCCGGGCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCttt                                                      >  1:1433959/1‑69 (MQ=255)
                                      ttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                 >  1:1103645/1‑68 (MQ=255)
                                      ttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                 >  1:2941412/1‑68 (MQ=255)
                                      ttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                 >  1:913356/1‑68 (MQ=255)
                                                ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt      <  1:1682500/69‑1 (MQ=255)
                                                ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt      <  1:3039060/69‑1 (MQ=255)
                                                ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt      <  1:3168823/69‑1 (MQ=255)
                                                ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt      <  1:449701/69‑1 (MQ=255)
                                                ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt      <  1:1540428/69‑1 (MQ=255)
                                                ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt      <  1:151370/69‑1 (MQ=255)
                                                ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt      <  1:746087/69‑1 (MQ=255)
                                                 gTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt      <  1:1374790/68‑1 (MQ=255)
                                                  tCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt      <  1:2887718/67‑1 (MQ=255)
                                                    cacaCCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCgtcgtc  >  1:2941958/1‑69 (MQ=255)
                                                    cacaCCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCgtcgtc  >  1:1137598/1‑69 (MQ=255)
                                                    cacaCCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCgtcgtc  >  1:3074596/1‑69 (MQ=255)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:714551/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:669198/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:985412/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:479860/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:3367665/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:3166815/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:3158919/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:1098081/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:2831582/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:2710983/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:210695/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:1884499/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:1841359/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:1782258/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:1567560/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:1480407/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:1301618/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:1291267/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:1213874/1‑42 (MQ=37)
                                                              gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaa                    >  1:3045739/1‑41 (MQ=37)
                                                              gCTCGTTACCAGACACTTTGCTAAGCGCCACTGTATCGGaaa                   >  1:486377/1‑42 (MQ=25)
                                                                  |                                                      
CCACAAGATCTCTGAGCCGGGCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTC  >  minE/1461088‑1461208

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: