Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,461,181 G→T 49.2% pseudogene (774/4605 nt)
pseudogene (774/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,1810GT49.2% 8.0 / 67.8 63pseudogene (774/4605 nt)
pseudogene (774/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 (8/24);  new base T (21/10);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 9.67e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

TTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAG  >  minE/1461126‑1461246
                                                       |                                                                 
ttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                                       >  1:1103645/1‑68 (MQ=255)
ttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                                       >  1:913356/1‑68 (MQ=255)
ttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                                       >  1:2941412/1‑68 (MQ=255)
          ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                                            <  1:3039060/69‑1 (MQ=255)
          ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                                            <  1:3168823/69‑1 (MQ=255)
          ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                                            <  1:746087/69‑1 (MQ=255)
          ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                                            <  1:1682500/69‑1 (MQ=255)
          ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                                            <  1:151370/69‑1 (MQ=255)
          ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                                            <  1:1540428/69‑1 (MQ=255)
          ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                                            <  1:449701/69‑1 (MQ=255)
           gTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                                            <  1:1374790/68‑1 (MQ=255)
            tCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                                            <  1:2887718/67‑1 (MQ=255)
              cacaCCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCgtcgtc                                        >  1:3074596/1‑69 (MQ=255)
              cacaCCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCgtcgtc                                        >  1:2941958/1‑69 (MQ=255)
              cacaCCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCgtcgtc                                        >  1:1137598/1‑69 (MQ=255)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                                         >  1:2831582/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                                         >  1:1098081/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                                         >  1:2710983/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                                         >  1:3158919/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                                         >  1:3166815/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                                         >  1:3367665/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                                         >  1:479860/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                                         >  1:669198/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                                         >  1:714551/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                                         >  1:985412/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                                         >  1:1567560/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:1480407/1‑42 (MQ=37)
                        gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa                                                         >  1:210695/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)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                                       <  1:480789/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                                       <  1:2555913/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                                       <  1:403017/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                                       <  1:1636331/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                                       <  1:1716169/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                                       <  1:3077862/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                                       <  1:931713/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                                       <  1:2994308/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                                       <  1:2089951/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                                       <  1:2733755/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca                                                       <  1:21149/38‑1 (MQ=255)
                              tGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACccacca                         >  1:519329/1‑68 (MQ=255)
                               accTGTCACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt                                            <  1:2854323/47‑1 (MQ=38)
                                   gACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACcatccatc                   <  1:2819340/69‑1 (MQ=255)
                                   gACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATc                   <  1:927964/69‑1 (MQ=255)
                                   gACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATc                   <  1:1381820/69‑1 (MQ=255)
                                        tttGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTc              <  1:1003611/69‑1 (MQ=255)
                                        tttGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTc              <  1:2053717/69‑1 (MQ=255)
                                        tttGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTc              <  1:2064975/69‑1 (MQ=255)
                                        tttGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTc              <  1:2512996/69‑1 (MQ=255)
                                        tttGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTc              <  1:2668998/69‑1 (MQ=255)
                                        tttGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTc              <  1:593831/69‑1 (MQ=255)
                                         ttGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTc              <  1:2795753/68‑1 (MQ=255)
                                         ttGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTc              <  1:3421980/68‑1 (MQ=255)
                                               aGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAgg       >  1:2891816/1‑69 (MQ=255)
                                                    cACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAg  <  1:1938640/69‑1 (MQ=255)
                                                    cACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAg  <  1:3035392/69‑1 (MQ=255)
                                                       |                                                                 
TTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAG  >  minE/1461126‑1461246

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: