Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 804,403 T→G 19.7% C238W (TGT→TGG)  U712_03910 → PTS system N‑acetylglucosamine‑specific EIICB component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881804,4030TG19.7% 40.8 / 3.1 39C238W (TGT→TGG) U712_03910PTS system N‑acetylglucosamine‑specific EIICB component
Reads supporting (aligned to +/- strand):  ref base T (14/17);  new base G (4/4);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATATTTTCTGGTTCCAATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGTTTGCATTTATGTTTTTATCTCCG  >  CP006881/804265‑804535
                                                                                                                                          |                                                                                                                                    
cAATATTTTCTGGTTCCTATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCgtgt                                                                                                                                      <  3:7607/139‑1 (MQ=255)
cAATATTATCTGGATCCAATTCGGAGAATATAACCGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCgtgt                                                                                                                                      <  4:375450/139‑1 (MQ=255)
        tCTGGTTCCAATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTAt                                                                                                                              >  7:18869/1‑139 (MQ=255)
        tCTGGTTCCAATTCGGAGAATATAACGGAGTGACAGCGGATTTGGCGCGCTTTTTCGCGAAAGATCCGAATGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGAGTCAggggat                                                                                                                              >  1:128752/1‑135 (MQ=255)
          tGGTTCCAATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATgg                                                                                                                            >  7:129926/1‑139 (MQ=255)
           ggTTCCAATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGGTATTGt                                                                                                                           >  6:315703/1‑139 (MQ=255)
             ttCCAATTCGGAGAATATAACGGAGTGTCAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGGCTGGCTATGGTtg                                                                                                                         >  1:76447/1‑139 (MQ=255)
             ttCCAATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTtg                                                                                                                         >  3:293222/1‑139 (MQ=255)
                         gaATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGc                                                                                                             <  8:283238/139‑1 (MQ=255)
                          aaTATACCGGAGTGACCGTGGATTTGGCGAGCTGTTTCGCGAAAGATCCGACTGCGGGCACATAACTGACCGGATGTTTCCCAATCATGCTGATCGGTCTTCCGGCGGCGTGTCTGGCTATGGTGGTGACTGCGAAGc                                                                                                             <  2:277806/138‑1 (MQ=255)
                             ataACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTGTTGCTATGGTTGTGACTGCGAAGCCGTc                                                                                                         >  4:374420/1‑139 (MQ=255)
                              tatcGGATTGACTTGTGATGTGTCGAGCTTTTTCGCGTAAGATCCGACTGCGGGCACATATATGACAGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTTCGAAGCCGTc                                                                                                         <  2:128752/135‑1 (MQ=255)
                                    aGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTGTGGCTATGGTTGTGACTGCGAAGCCGTCTAAACgg                                                                                                  >  2:117832/1‑139 (MQ=255)
                                    aGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGGCTGGCTATGGTTGTGACTGCGAAGCCGTCGAAACgg                                                                                                  >  7:364631/1‑139 (MQ=255)
                                      tGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTTTGGTTGTGACTGCGAAGCCGTCTAAACGGa                                                                                                 >  4:346019/1‑138 (MQ=255)
                                      tGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTTACTGCGAAGCCGGCAAAACGGa                                                                                                 >  6:2039/1‑138 (MQ=255)
                                                tttGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCGATGGTTGTGCCTGCGAAGCCGTCTAAACGGAAAGCGACGGCt                                                                                      >  4:194729/1‑139 (MQ=255)
                                                     cgcgcTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAat                                                                                 <  2:76447/139‑1 (MQ=255)
                                                                 cgAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGAtt                                                                     <  6:185815/139‑1 (MQ=255)
                                                                 cgAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGGGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGAtt                                                                     <  8:18869/139‑1 (MQ=255)
                                                                 cgAAAGATCCGACTGCGGGCACATAAATGCCCGGATTTTTCCAAAGCATGATGTTCGGGCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGAtt                                                                     <  4:373883/139‑1 (MQ=255)
                                                                           gACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGATTTGCTTTAACg                                                           <  8:355466/139‑1 (MQ=255)
                                                                           gACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGGCTTCCGGCGGCGGGGCTGGCTATGGTGGTGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGATTTGCTTTAACg                                                           <  1:114066/139‑1 (MQ=255)
                                                                                cGGGCACATATATGCCCGGATTTTTCCCAATCATGAGGTGCGGTCTGCCGGCGGGGTGTCTGGCTATGGTTGTGACGGCGAAGCCGGCTAAACGGAAAGCGACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGtt                                                      <  2:96743/139‑1 (MQ=255)
                                                                                   gCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTAt                                                   <  8:364631/139‑1 (MQ=255)
                                                                                      ctaatatGTCGGGAGTTTTCCCATTCATGATGTTCGGTCTTCCGGGTGCGTGTCTGGCTATGGTTGTGCCTGCGTAGCCGGCTAAACGTAAAGCGACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATcac                                                <  5:33764/136‑1 (MQ=255)
                                                                                      caTATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATcac                                                <  4:293222/139‑1 (MQ=255)
                                                                                      caTATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGGCTGGCTATGGTTGTGACTGCGAAGCCGGCTAAACGGAACGCGACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATcac                                                <  4:66248/139‑1 (MQ=255)
                                                                                          tatGACCGGATTTTTCCCAATCATGATGTTCGTGCTTCCGGCGGCGGGTCTGGCTATGGTTGTGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAgg                                             <  5:315703/138‑1 (MQ=255)
                                                                                          tatGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAgg                                             <  8:113863/138‑1 (MQ=255)
                                                                                          taagACCTGATTTTTCCCAATCATGAGGTTCTGTCTTCCGGCGGCGTGTCTGTCTATGGTTGGGCCTGCTAAGCCGTAAAAACGGAAAGCGACGGCTGTAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAgg                                             <  3:244211/135‑1 (MQ=255)
                                                                                               ccGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTGTGGCTATGGTTGTGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTAc                                       >  8:316649/1‑139 (MQ=255)
                                                                                                 ggATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGGGGTGGGATGGTTGTGACTGCGAAGCCGGCGAAACGGAAAGCGACGGCGGGAATGATGTTTGGATTTGGGTTAACCGGGTTGATGACAGGGATTTCgg                                     >  8:89762/1‑139 (MQ=255)
                                                                                                  gATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACg                                      >  8:325360/1‑137 (MQ=255)
                                                                                                  gATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGAATGCGAAGCCGTCAAAACGGAAATCGACCGCGGGAATGAAGATTGGATTTGCTTTAACGGCGTTTTTAAAAGGGAATaac                                      >  6:318290/1‑135 (MQ=255)
                                                                                                   aTTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGGGTGGCGATGGTTGTGACAGGGAACCCGTGTAAACCGAAAGCGGCGGCCGGTATGATGATTGGATTTGGGTTAACGGGGGTTTTAGCCGGGATTTCggcg                                   >  6:363153/1‑137 (MQ=255)
                                                                                                                      ttCGGTCTTCCGGCGGCGTGTCTGGCTGTGGTTGTGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGTTTGCATTTa                >  3:90304/1‑139 (MQ=255)
                                                                                                                      tccGGTCTTCCGGCGGCGGGGCTGGCTTTGGTTTGTCCTGCGAAGCGGTCTAAACGGAACTCGACGGCTGGAATGATGATTGGATTTGCTTTAACGGCTTTTATCCCAGGGATTACGGAGACGATTGAGTTTGCAttt                 <  6:133139/136‑1 (MQ=255)
                                                                                                                                  gcggcgTGTCTGGCTTTGGTGTTGTCTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGTTTGCATTTATGTTTTTAtctc    <  3:374420/139‑1 (MQ=255)
                                                                                                                                  gcggcgTGTCTGGCTATGGTTGTGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGTTTGCATTTATGTTTTTAtctc    <  3:346019/139‑1 (MQ=255)
                                                                                                                                    gtcgTGGCTGGCTAGGGTTGTGACTTCGATGCCGTCAAAACGGAAAGCGACGGCTGGAATGAAGATTGGAATTGCTTTAAAGGCGTTTATCACAGTGATTACGGAGCCGATTGAGTTTGCATTTATGTTTTTATCTCCg  <  4:83837/137‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                    
CAATATTTTCTGGTTCCAATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGTCTAAACGGAAAGCGACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGTTTGCATTTATGTTTTTATCTCCG  >  CP006881/804265‑804535

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