Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,914,611 A→G 34.4% S287S (AGT→AGC 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,6110AG34.4% 39.0 / 30.7 58S287S (AGT→AGCptsPfused PEP‑protein phosphotransferase (enzyme I) of PTS system
Reads supporting (aligned to +/- strand):  ref base A (11/27);  new base G (8/12);  total (19/39)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCC  >  minE/1914545‑1914679
                                                                  |                                                                    
aaGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCa                                                                      >  1:2275962/1‑67 (MQ=255)
aaGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCa                                                                      >  1:2072805/1‑67 (MQ=255)
aaGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCa                                                                      >  1:1971602/1‑67 (MQ=255)
    tGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAgc                                                                >  1:1731419/1‑69 (MQ=255)
                 gCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCGCTTAgc                                                                >  1:3275657/1‑56 (MQ=255)
                 gCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGt                                                    >  1:572573/1‑68 (MQ=255)
                 gCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGt                                                    >  1:1987450/1‑68 (MQ=255)
                 gCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGt                                                    >  1:2692936/1‑68 (MQ=255)
                 gCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGt                                                    >  1:3036258/1‑68 (MQ=255)
                 gCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGt                                                    >  1:1856020/1‑68 (MQ=255)
                 gCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCGCTTAGCGAGGCAAACTGt                                                    >  1:1999085/1‑68 (MQ=255)
                 gCTGACCCAGCGCACGTAAATCGACAGCCCGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGt                                                    >  1:2284664/1‑68 (MQ=255)
                              aCGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTc                                      <  1:1555980/69‑1 (MQ=255)
                               cGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTc                                      <  1:2061514/68‑1 (MQ=255)
                                          aGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGtttt                          <  1:423068/69‑1 (MQ=255)
                                          aGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGtttt                          <  1:842448/69‑1 (MQ=255)
                                           gCCCGCTCTTTGAGATAGTTATCACTTAGCGCGTCAAACTGTTCGTCAAATTTTTCAATGACCGttttt                         <  1:118358/69‑1 (MQ=255)
                                           gCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGttttt                         <  1:2767215/69‑1 (MQ=255)
                                           gCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGttttt                         <  1:3130079/69‑1 (MQ=255)
                                           gCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGttttt                         <  1:3270777/69‑1 (MQ=255)
                                           gCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGttttt                         <  1:74822/69‑1 (MQ=255)
                                           gCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGttttt                         <  1:1121727/69‑1 (MQ=255)
                                           gCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGttttt                         <  1:713966/69‑1 (MQ=255)
                                           gCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGttttt                         <  1:1215559/69‑1 (MQ=255)
                                           gCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGttttt                         <  1:34602/69‑1 (MQ=255)
                                           gCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGttttt                         <  1:1251809/69‑1 (MQ=255)
                                           gCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGttttt                         <  1:3369390/69‑1 (MQ=255)
                                           gCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGtttt                          <  1:256696/68‑1 (MQ=255)
                                           gCCCGCTCTTTGAGATAGTTATCACTTAGCGAGGCAAACTGTTCGGCAAATTTTTCAATGACCGttttt                         <  1:1443663/69‑1 (MQ=255)
                                              cGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGtt                                                   <  1:2931727/40‑1 (MQ=255)
                                              cGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGtt                                                   <  1:2956791/40‑1 (MQ=255)
                                              cGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGtt                                                   <  1:3355200/40‑1 (MQ=255)
                                              cGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGtt                                                   <  1:688501/40‑1 (MQ=255)
                                              cGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGtt                                                   <  1:2679122/40‑1 (MQ=255)
                                              cGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGtt                                                   <  1:3289292/40‑1 (MQ=255)
                                              cGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGtt                                                   <  1:3054509/40‑1 (MQ=255)
                                              cGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGtt                                                   <  1:103226/40‑1 (MQ=255)
                                              cGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGtt                                                   <  1:253221/40‑1 (MQ=255)
                                              cGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGtt                                                   <  1:22585/40‑1 (MQ=255)
                                              cGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGtt                                                   <  1:1197252/40‑1 (MQ=255)
                                              cGCTCTTTGAGATAGTTATCGCTTAGCGCGGCAAACTGtt                                                   <  1:1056302/40‑1 (MQ=255)
                                                 tctTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTc                                      >  1:906230/1‑50 (MQ=255)
                                                 tctTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTc                                      >  1:902270/1‑50 (MQ=255)
                                                 tctTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTc                                      >  1:856675/1‑50 (MQ=255)
                                                 tctTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTc                                      >  1:852689/1‑50 (MQ=255)
                                                 tctTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTc                                      >  1:3119082/1‑50 (MQ=255)
                                                 tctTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGcc                   <  1:1264213/69‑1 (MQ=255)
                                                  ctTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGt                             <  1:2493802/58‑1 (MQ=255)
                                                     tGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCa                 <  1:3043700/67‑1 (MQ=255)
                                                      gagaTAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGtt                            <  1:187218/55‑1 (MQ=255)
                                                      gagaTAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGtt                            <  1:670914/55‑1 (MQ=255)
                                                      gagaTAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGtt                            <  1:2905091/55‑1 (MQ=255)
                                                      gagaTAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGtt                            <  1:1495847/55‑1 (MQ=255)
                                                      gagaTAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGtt                            <  1:2891084/55‑1 (MQ=255)
                                                      gagaTAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGtt                            <  1:2062433/55‑1 (MQ=255)
                                                      gagaTAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGtt                            <  1:486387/55‑1 (MQ=255)
                                                              tATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCg      >  1:1013370/1‑69 (MQ=255)
                                                              tATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCg      >  1:3063237/1‑69 (MQ=255)
                                                              tATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCg      >  1:520255/1‑69 (MQ=255)
                                                               aTCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGa     >  1:125066/1‑69 (MQ=255)
                                                                 cACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGc   >  1:1753285/1‑69 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAAt                                   >  1:74136/1‑36 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAAt                                   >  1:1487229/1‑36 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAAt                                   >  1:1177239/1‑36 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAAt                                   >  1:1968235/1‑36 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAAt                                   >  1:1901160/1‑36 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAAt                                   >  1:1555253/1‑36 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAAt                                   >  1:125011/1‑36 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAAt                                   >  1:3008963/1‑36 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGATcc  <  1:1167117/69‑1 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGcc  <  1:1825223/69‑1 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGcc  <  1:926456/69‑1 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGcc  <  1:1929874/69‑1 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGcc  <  1:2949119/69‑1 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGcc  <  1:843379/69‑1 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGcc  <  1:2974762/69‑1 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGcc  <  1:761428/69‑1 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGcc  <  1:1691436/69‑1 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGcc  <  1:1517066/69‑1 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGcc  <  1:539763/69‑1 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGcc  <  1:3315362/69‑1 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGcc  <  1:503617/69‑1 (MQ=255)
                                                                  aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGcc  <  1:3409831/69‑1 (MQ=255)
                                                                  |                                                                    
AAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCC  >  minE/1914545‑1914679

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: