Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,914,596 T→G 11.7% K292N (AAA→AAC 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,5960TG11.7% 130.0 / 14.8 60K292N (AAA→AACptsPfused PEP‑protein phosphotransferase (enzyme I) of PTS system
Reads supporting (aligned to +/- strand):  ref base T (27/26);  new base G (5/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

CCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCC  >  minE/1914527‑1914662
                                                                     |                                                                  
cccTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCtt                                                                   <  1:2738593/71‑1 (MQ=255)
 ccTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCgtt                                                                  <  1:645265/71‑1 (MQ=255)
 ccTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCgtt                                                                  >  1:1949734/1‑71 (MQ=255)
  cTTGATTAGCGTCATCAAGATGAAACAGCAATCTCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCttt                                                                  >  1:3011878/1‑70 (MQ=255)
  cTTGATGACGGTAATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCttt                                                                  >  1:10264/1‑70 (MQ=255)
    tGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCttt                                                                  >  1:1036195/1‑68 (MQ=255)
     gATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCttt                                                                  >  1:126892/1‑67 (MQ=255)
      aTTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCGTTgaga                                                              <  1:610399/70‑1 (MQ=255)
         aGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGt                                                          >  1:3316572/1‑71 (MQ=255)
         aGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGt                                                          >  1:2987585/1‑71 (MQ=255)
          gCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGt                                                          >  1:519013/1‑70 (MQ=255)
           cGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTa                                                        >  1:2670947/1‑71 (MQ=255)
              catcaAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCa                                                     >  1:1399041/1‑71 (MQ=255)
               atcaAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCa                                                     <  1:2023173/70‑1 (MQ=255)
                tcaAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTga                                                                <  1:1342895/58‑1 (MQ=255)
                  aaGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTAt                                                       >  1:2948455/1‑65 (MQ=255)
                  aaGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACtt                                                  <  1:710778/70‑1 (MQ=255)
                  aaGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACt                                                   >  1:3510027/1‑69 (MQ=255)
                   aGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAg                                                >  1:3441845/1‑71 (MQ=255)
                      tGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAgc                                               >  1:2692536/1‑69 (MQ=255)
                         aaCAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCgg                                           >  1:919427/1‑70 (MQ=255)
                         aaCAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCgg                                           >  1:1951227/1‑70 (MQ=255)
                          aCAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGc                                          >  1:1666671/1‑70 (MQ=255)
                          aCAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGc                                          >  1:2627759/1‑70 (MQ=255)
                          aCAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCa                                         <  1:2194515/71‑1 (MQ=255)
                          aCAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCa                                         <  1:3032739/71‑1 (MQ=255)
                              cAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGAt                                                             >  1:3233661/1‑47 (MQ=255)
                              cAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTAt                                                       >  1:822162/1‑53 (MQ=255)
                              cAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACt                                     >  1:1225366/1‑71 (MQ=255)
                               aaTCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTg                                    <  1:3291589/71‑1 (MQ=255)
                               aaTCGCTGACCCAGCGCACGTAAATCGCCAGCCAGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTg                                    >  1:342685/1‑71 (MQ=255)
                                aTCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTgaga                                                              <  1:455698/44‑1 (MQ=255)
                                aTCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATa                                                            <  1:136011/46‑1 (MQ=255)
                                 tCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGtt                                  <  1:665325/71‑1 (MQ=255)
                                 tCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGtt                                  <  1:1336536/71‑1 (MQ=255)
                                 tCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGtt                                  <  1:336762/71‑1 (MQ=255)
                                 tCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCGTTGAGATAGTTATCACTTAGCGCGGCAAACTGtt                                  <  1:2672754/71‑1 (MQ=255)
                                  cGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTc                                 >  1:3059825/1‑71 (MQ=255)
                                  cGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCGTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTc                                 >  1:3867702/1‑71 (MQ=255)
                                  cGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCGTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTc                                 >  1:328773/1‑71 (MQ=255)
                                   gCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCg                                <  1:2473509/71‑1 (MQ=255)
                                   gCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCGTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCg                                >  1:2289466/1‑71 (MQ=255)
                                    cTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCAc                                                    >  1:625755/1‑50 (MQ=255)
                                     tGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCTCGGCAAACTGTTCgg                               >  1:1670026/1‑70 (MQ=255)
                                     tGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCgga                              >  1:1130362/1‑70 (MQ=255)
                                      gACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGc                              >  1:1620572/1‑70 (MQ=255)
                                       aCCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCa                             <  1:243539/70‑1 (MQ=255)
                                           aGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCa                                                     <  1:1549206/42‑1 (MQ=255)
                                             cgcACGTAAATCGCCAGCCCGTTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAAttttt                      <  1:597973/71‑1 (MQ=255)
                                             cgcACGTAAATCGCCAGCCCGCTCTTTGAGATATTTAt                                                       <  1:1910067/38‑1 (MQ=255)
                                                aCGTAAATCGCCAGCCCGCTCGTTGAGATAGTTATCACTTAGCGCGGCAAAc                                      >  1:1795266/1‑52 (MQ=255)
                                                  gTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACtt                                                  >  1:3381337/1‑38 (MQ=255)
                                                   tAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAgcg                                              >  1:1696927/1‑41 (MQ=255)
                                                   tAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCa                             <  1:3839441/58‑1 (MQ=255)
                                                   tAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAAttttt                      <  1:394930/65‑1 (MQ=255)
                                                       tCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGt            <  1:1566837/71‑1 (MQ=255)
                                                         gCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCa                             <  1:3405037/52‑1 (MQ=255)
                                                         gCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCa                             <  1:1001236/52‑1 (MQ=255)
                                                         gCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGttt          >  1:3200049/1‑71 (MQ=255)
                                                            aGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTa       <  1:970479/71‑1 (MQ=255)
                                                             gCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAAtttt                       <  1:3188807/54‑1 (MQ=255)
                                                             gCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAAtttt                       <  1:3192888/54‑1 (MQ=255)
                                                             gCCCGCTCGTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAAtttt                       >  1:2879989/1‑54 (MQ=255)
                                                               ccGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATg                 <  1:3679620/58‑1 (MQ=255)
                                                                cGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCatat                          <  1:358518/48‑1 (MQ=255)
                                                                 gCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGcc  >  1:1456552/1‑71 (MQ=255)
                                                                 gCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGcc  >  1:2117551/1‑71 (MQ=255)
                                                                   tctTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAAttttt                      <  1:294600/49‑1 (MQ=255)
                                                                     |                                                                  
CCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCC  >  minE/1914527‑1914662

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: