Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,898,828 (T)5→4 coding (130/1878 nt) bglF ← PTS beta‑glucoside transporter subunit IIABC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,898,8240T.100.0% 85.0 / NA 21coding (134/1878 nt)bglFPTS beta‑glucoside transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (12/9);  total (12/9)

CCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTT  >  NZ_CP009273/3898741‑3898909
                                                                                   |                                                                                     
ccAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttcag                                                                                 >  1:59298/1‑87 (MQ=255)
  aGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTa                                                                              <  2:319499/90‑1 (MQ=255)
                 acaTGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCgctt                                                               <  2:59298/90‑1 (MQ=255)
                          ccTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCAt                                                      <  1:243816/90‑1 (MQ=255)
                            tATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATcc                                                    >  2:50970/1‑90 (MQ=255)
                            tATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATcc                                                    >  2:103368/1‑90 (MQ=255)
                               gACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCttt                                                 >  1:94539/1‑90 (MQ=255)
                                             ggCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTaa                                   >  1:324101/1‑90 (MQ=255)
                                                      tttCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgc                          <  1:312958/90‑1 (MQ=255)
                                                       ttCCACCACCATAATAATACCGGGGGTCTTTTCAGTACct                                                                           >  1:117932/1‑40 (MQ=39)
                                                       ttCCACCACCATAATAATACCGGGGGTCTTTTCAGTACct                                                                           <  2:117932/40‑1 (MQ=255)
                                                       ttCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcg                         >  1:194348/1‑90 (MQ=255)
                                                          caccaccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa                      <  2:241284/90‑1 (MQ=255)
                                                               ccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCAt                 >  2:207917/1‑90 (MQ=255)
                                                               ccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCAt                 >  1:160516/1‑90 (MQ=255)
                                                                caTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATc                <  2:194348/90‑1 (MQ=255)
                                                                   aataataCCGGGGGTCTTTTCAGTACCTCTGCTTGCgctt                                                               <  1:286237/40‑1 (MQ=255)
                                                                   aataataCCGGGGGTCTTTTCAGTACCTCTGCTTGCgctt                                                               >  2:286237/1‑40 (MQ=255)
                                                                      aataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa                      <  2:207747/78‑1 (MQ=255)
                                                                      aataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa                      >  1:207747/1‑78 (MQ=255)
                                                                            gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg    >  1:299832/1‑90 (MQ=255)
                                                                              gggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGtt  >  1:76850/1‑90 (MQ=255)
                                                                                   |                                                                                     
CCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTT  >  NZ_CP009273/3898741‑3898909

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: