Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,461,157 G→A 62.5% pseudogene (798/4605 nt)
pseudogene (798/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,1570GA62.5% 22.6 / 28.6 48pseudogene (798/4605 nt)
pseudogene (798/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 G (7/11);  new base A (21/9);  total (28/20)
Fisher's exact test for biased strand distribution p-value = 6.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCA  >  minE/1461126‑1461223
                               |                                                                  
ttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                >  1:1103645/1‑68 (MQ=255)
ttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                >  1:913356/1‑68 (MQ=255)
ttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                >  1:2941412/1‑68 (MQ=255)
          ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                     <  1:746087/69‑1 (MQ=255)
          ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                     <  1:449701/69‑1 (MQ=255)
          ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                     <  1:3168823/69‑1 (MQ=255)
          ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                     <  1:3039060/69‑1 (MQ=255)
          ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                     <  1:151370/69‑1 (MQ=255)
          ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                     <  1:1540428/69‑1 (MQ=255)
          ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                     <  1:1682500/69‑1 (MQ=255)
           gTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                     <  1:1374790/68‑1 (MQ=255)
            tCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                     <  1:2887718/67‑1 (MQ=255)
              cacaCCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCgtcgtc                 >  1:1137598/1‑69 (MQ=255)
              cacaCCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCgtcgtc                 >  1:3074596/1‑69 (MQ=255)
              cacaCCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCgtcgtc                 >  1:2941958/1‑69 (MQ=255)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:985412/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:1213874/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:714551/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:669198/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:479860/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:1291267/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:3367665/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:1301618/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:3166815/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:3158919/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:1480407/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:1567560/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:1782258/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:1841359/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:1884499/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:210695/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaa                                   >  1:3045739/1‑41 (MQ=37)
                        gCTCGTTACCAGACACTTTGCTAAGCGCCACTGTATCGGaaa                                  >  1:486377/1‑42 (MQ=25)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                <  1:1636331/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                <  1:931713/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                <  1:1716169/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                <  1:480789/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                <  1:2089951/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                <  1:403017/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                <  1:21149/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                <  1:2555913/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                <  1:2733755/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                <  1:3077862/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                <  1:2994308/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACccacca  >  1:519329/1‑68 (MQ=255)
                               |                                                                  
TTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCA  >  minE/1461126‑1461223

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: