Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,966,215 T→G 56.6% 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
*W3110S.gb2,966,2150TG56.6% 23.8 / 18.9 23K292N (AAA→AACptsPfused PEP‑protein phosphotransferase (enzyme I) of PTS system
Reads supporting (aligned to +/- strand):  ref base T (3/7);  new base G (2/11);  total (5/18)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

CTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTT  >  W3110S.gb/2966148‑2966272
                                                                   |                                                         
cTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCttt                                                         >  1:1125176/1‑70 (MQ=255)
   gATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCttt                                                         >  1:978558/1‑67 (MQ=255)
    aTTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTg                                                        <  1:1310400/67‑1 (MQ=255)
    aTTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCGTTGAGAt                                                    <  1:24938/71‑1 (MQ=255)
              tcaAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCGTTGAGATAGTTATCACt                                          <  1:1146289/71‑1 (MQ=255)
               caAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCGTTGAGATAGTTATCACt                                          >  1:1531407/1‑70 (MQ=255)
                  gATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCGTTGAGAt                                                    <  1:993779/57‑1 (MQ=255)
                  gATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCGTTGAGATAGTTATCACTTAgc                                      <  1:766522/71‑1 (MQ=255)
                   aTGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCGTTGAGATAGTTATCACTTAgcg                                     <  1:1524232/71‑1 (MQ=255)
                              aTCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGt                          <  1:665122/71‑1 (MQ=255)
                               tCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTa                                               <  1:31182/49‑1 (MQ=255)
                               tCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCa                                <  1:236233/64‑1 (MQ=255)
                                       ccAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCaaa                  <  1:1606418/70‑1 (MQ=255)
                                       ccAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCaaa                  <  1:172820/70‑1 (MQ=255)
                                            gcACGTAAATCGCCAGCCCGCTCTTTGAGATATTTATCACTTAGCGCGGCAAACTGTTCGGCAAAttttt             >  1:1067226/1‑70 (MQ=255)
                                            gcACGTAAATCGCCAGCCCGCTCGTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTc            <  1:443924/71‑1 (MQ=255)
                                             cACGTAAATCGCCAGCCCGCTCGTTGAGATAGTTATCACTTAGCGCGGCAAACTGtt                         <  1:1600622/57‑1 (MQ=255)
                                             cACGTAAATCGCCAGCCCGCTCGTTGAGATAGTTATCACTTAGCGCGGCAAACTGtt                         <  1:1699849/57‑1 (MQ=255)
                                             cACGTAAATCGCCAGCCCGCTCGTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCa           <  1:1494214/71‑1 (MQ=255)
                                                 tAAATCGCCAGCCCGCTCGTTGAGATAGTTATCACtt                                         <  1:1310817/37‑1 (MQ=255)
                                                  aaaTCGCCAGCCCGCTCGTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGAc      >  1:1245273/1‑71 (MQ=255)
                                                     tCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAgcg                                     <  1:530400/37‑1 (MQ=255)
                                                     tCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACt                            >  1:507237/1‑46 (MQ=255)
                                                     tCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACt                            >  1:877564/1‑46 (MQ=255)
                                                       gCCAGCCCGCTCGTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGtt  <  1:355168/70‑1 (MQ=255)
                                                                   |                                                         
CTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTT  >  W3110S.gb/2966148‑2966272

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: