Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,914,848 T→C 56.8% E208E (GAA→GAG ptsP ← fused PEP‑protein phosphotransferase (enzyme I) of PTS system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,914,8480TC56.8% 16.5 / 43.6 37E208E (GAA→GAGptsPfused PEP‑protein phosphotransferase (enzyme I) of PTS system
Reads supporting (aligned to +/- strand):  ref base T (7/9);  new base C (15/6);  total (22/15)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGC  >  minE/1914780‑1914914
                                                                    |                                                                  
gcAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTc                                                                  <  1:2218858/71‑1 (MQ=255)
                aGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAg                                                   <  1:270757/70‑1 (MQ=255)
                aGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTCCAGAGc                                                            >  1:469718/1‑61 (MQ=255)
                  cggcgAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGt                                                <  1:2930083/71‑1 (MQ=255)
                   ggcgAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGtt                                               <  1:157663/71‑1 (MQ=255)
                   ggcgAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGtt                                               <  1:2601684/71‑1 (MQ=255)
                   ggcgAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGtt                                               <  1:2357075/71‑1 (MQ=255)
                   ggcgAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTCCAGAGCCGGATCCAGCGtt                                               <  1:1371325/71‑1 (MQ=255)
                   ggcgAAACTCGTTTGCAGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTCCAGAGCCGGATCCAGCGtt                                               <  1:3176555/71‑1 (MQ=255)
                     cgAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTCCAGAGCCGGATCCAGCGTTGa                                             >  1:2208824/1‑71 (MQ=255)
                      gAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTCCAGAGCCGGATCCAGCGTTGAt                                            >  1:3546388/1‑71 (MQ=255)
                       aaaCTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGtt                                               <  1:1393958/67‑1 (MQ=255)
                       aaaCTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGtt                                               <  1:4233492/67‑1 (MQ=255)
                         aCTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGa                                             >  1:807336/1‑67 (MQ=255)
                         aCTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGa                                             >  1:4158726/1‑67 (MQ=255)
                         aCTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTCCAGAGCCGGATCCAGCGtt                                               <  1:2673209/65‑1 (MQ=255)
                             gTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGtt                                               >  1:578220/1‑61 (MQ=255)
                             gTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGtt                                               >  1:3909709/1‑61 (MQ=255)
                             gTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTcc                                                                 >  1:2351428/1‑43 (MQ=255)
                             gTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTcc                                                                 >  1:2536815/1‑43 (MQ=255)
                             gTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTcc                                                                 >  1:945773/1‑43 (MQ=255)
                             gTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTcc                                                                 >  1:66743/1‑43 (MQ=255)
                             gTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTcc                                                                 >  1:565610/1‑43 (MQ=255)
                                    gCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTg                              >  1:903593/1‑71 (MQ=255)
                                      ttcttcCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGtt                                               <  1:2518895/52‑1 (MQ=255)
                                        cttcCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTc                           >  1:101830/1‑70 (MQ=255)
                                        cttcCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTc                           >  1:935547/1‑70 (MQ=255)
                                                    cGGTCAGTCGTTCGCGCTCCAGAGCCGGATCCAGCGTTGa                                             >  1:47729/1‑40 (MQ=255)
                                                    cGGTCAGTCGTTCGCGCTCCAGAGCCGGATCCAGCGTTGa                                             >  1:763849/1‑40 (MQ=255)
                                                    cGGTCAGTCGTTCGCGCTCCAGAGCCGGATCCAGCGTTGATGCCTGa                                      >  1:3521917/1‑47 (MQ=255)
                                                    cGGTCAGTCGTTCGCGCTCCAGAGCCGGATCCAGCGTTGATGCCTGa                                      >  1:1969137/1‑47 (MQ=255)
                                                        cAGTCGTTCGCGCTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTgg          >  1:2549104/1‑71 (MQ=255)
                                                          gTCGTTCGCGCTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCa        <  1:3234189/71‑1 (MQ=255)
                                                           tCGTTCGCGCTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTa                      <  1:2876402/56‑1 (MQ=255)
                                                           tCGTTCGCGCTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTa                      <  1:949908/56‑1 (MQ=255)
                                                                cgcgCTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTg   >  1:386556/1‑70 (MQ=255)
                                                                cgcgCTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGc  >  1:2095325/1‑71 (MQ=255)
                                                                    |                                                                  
GCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGC  >  minE/1914780‑1914914

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: