Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,914,848 T→C 28.0% 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,8480TC28.0% 32.9 / 15.9 25E208E (GAA→GAGptsPfused PEP‑protein phosphotransferase (enzyme I) of PTS system
Reads supporting (aligned to +/- strand):  ref base T (8/10);  new base C (2/5);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

AAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCA  >  minE/1914782‑1914916
                                                                  |                                                                    
aaaGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCa                                                                  <  1:1694603/71‑1 (MQ=255)
   gCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCagag                                                               >  1:1529660/1‑71 (MQ=255)
         gCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGa                                                          <  1:271773/70‑1 (MQ=255)
            gTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTCCAGAGCCgg                                                           <  1:1643726/66‑1 (MQ=255)
             tAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCa                                                      <  1:25986/70‑1 (MQ=255)
             tAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAg                                                     >  1:824516/1‑71 (MQ=255)
              aGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTAAGTCGTTCGCGTTCCAGAGCCGGATCCAGc                                                    >  1:134384/1‑71 (MQ=255)
                 ggcgAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGtt                                                 <  1:1946880/71‑1 (MQ=255)
                 ggcgAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTCCAGAGCCGGATCCAGCGtt                                                 <  1:1279057/71‑1 (MQ=255)
                 ggcgAAACTCGTTTGCCGATTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTCCAGAGCCGGATCCAGCGtt                                                 <  1:465902/71‑1 (MQ=255)
                  gctaaaCTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGtt                                                 <  1:1851587/67‑1 (MQ=255)
                  gcgAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTcc                                                                   <  1:215583/52‑1 (MQ=255)
                  gcgAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGtt                                                 <  1:2163648/70‑1 (MQ=255)
                  gcgAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTCCAGAGCCGGATCCAGCGtt                                                 >  1:1359854/1‑70 (MQ=255)
                        cTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGtt                                                 <  1:310882/64‑1 (MQ=255)
                        cTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGCTCCAGAGCCGGATCCAGCGTTGATGCCt                                          <  1:846160/71‑1 (MQ=255)
                                        tcCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCAt                          >  1:211775/1‑71 (MQ=255)
                                                  cGGTCAGTCGTTCGCGCTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAAc                >  1:228413/1‑71 (MQ=255)
                                                         tCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCAt         >  1:2193118/1‑71 (MQ=255)
                                                         tCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCAt         >  1:1530771/1‑71 (MQ=255)
                                                           gTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGa                                        >  1:1225489/1‑38 (MQ=255)
                                                           gTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCa          <  1:2047309/68‑1 (MQ=255)
                                                            ttCGCGCTCCAGAGCCGGATCCAGCGTTGATGCCt                                          <  1:61790/35‑1 (MQ=255)
                                                              cgcgTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTg             >  1:1686784/1‑62 (MQ=255)
                                                                cgTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCa  <  1:1873263/71‑1 (MQ=255)
                                                                  |                                                                    
AAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCA  >  minE/1914782‑1914916

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: