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.98.2% 224.5 / ‑2.7 54coding (134/1878 nt)bglFPTS beta‑glucoside transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (35/18);  minor base C (0/1);  total (35/19)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG  >  NZ_CP009273/3898735‑3898907
                                                                                         |                                                                                   
tAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCt                                                                                     >  1:523073/1‑90 (MQ=255)
tAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCt                                                                                     <  2:697427/90‑1 (MQ=255)
tAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCt                                                                                     <  2:602433/90‑1 (MQ=255)
 aaCCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCtt                                                                                    >  2:528750/1‑90 (MQ=255)
    cGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttc                                                                                 <  2:216208/90‑2 (MQ=255)
                aTCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCt                                                                    >  2:353688/1‑90 (MQ=255)
                 tCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCtt                                                                   >  1:11512/1‑90 (MQ=255)
                  cGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTg                                                                  <  1:506890/90‑1 (MQ=255)
                      cacaTGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCgct                                                              >  1:441162/1‑90 (MQ=255)
                              tACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTc                                                      >  1:367721/1‑90 (MQ=255)
                               aCCTATGACCACCTGAATCTGGCCACCGCTTGCCCCCCCCATAACAAAACCGGGGGCCTTTTCAGTACCTCTTCCTGCGCTTTGCTTTCa                                                     <  1:20612/90‑1 (MQ=255)
                                ccTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCAt                                                    <  1:480866/90‑1 (MQ=255)
                                  tATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATcc                                                  >  2:429500/1‑90 (MQ=255)
                                  tATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATcc                                                  >  1:349103/1‑90 (MQ=255)
                                  tATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATcc                                                  >  2:208391/1‑90 (MQ=255)
                                    tGACCACCTTAAACTGGCCACCGCTTTCCCCCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtt                                                <  1:212078/90‑1 (MQ=255)
                                    tGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttcag                                                                               <  2:552769/60‑4 (MQ=255)
                                    tGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttcag                                                                               >  1:552769/1‑57 (MQ=255)
                                          ccTGAAACTGGCCACCGCTTTCCCCCCCCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaattt                                          <  1:560205/90‑1 (MQ=255)
                                          ccTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaattt                                          >  2:557186/1‑90 (MQ=255)
                                            tGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaatttaa                                        <  1:755173/90‑1 (MQ=255)
                                                   ggCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTaa                                 >  2:644449/1‑90 (MQ=255)
                                                    gCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAATACGTAAt                                >  2:389861/1‑90 (MQ=255)
                                                    gCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                <  2:214108/90‑1 (MQ=255)
                                                       aCCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCCTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATcgc                             >  1:371908/1‑90 (MQ=255)
                                                             ttCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcg                       >  2:501000/1‑90 (MQ=255)
                                                             ttCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcg                       >  2:758048/1‑90 (MQ=255)
                                                              tCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcgc                      <  2:335678/90‑1 (MQ=255)
                                                                caccaccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa                    <  2:198672/90‑1 (MQ=255)
                                                                caccaccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa                    <  2:687617/90‑1 (MQ=255)
                                                                  ccaccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATg                  >  1:642571/1‑90 (MQ=255)
                                                                  ccaccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATg                  >  1:445148/1‑90 (MQ=255)
                                                                     ccaTAATAATACTGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCAt               >  1:38717/1‑90 (MQ=255)
                                                                     ccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCAt               >  1:282836/1‑90 (MQ=255)
                                                                       ataataataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                <  1:478179/71‑1 (MQ=255)
                                                                       ataataataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                >  2:478179/1‑71 (MQ=255)
                                                                       ataataataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCa             >  2:126309/1‑90 (MQ=255)
                                                                             ataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca       <  2:99797/90‑1 (MQ=255)
                                                                             ataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca       <  2:65102/90‑1 (MQ=255)
                                                                             ataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca       <  1:398296/90‑1 (MQ=255)
                                                                             ataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca       >  1:436501/1‑90 (MQ=255)
                                                                              taCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcac      <  2:164225/90‑1 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca     >  1:106900/1‑90 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca     >  1:461601/1‑90 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca     >  2:481398/1‑90 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca     >  1:648199/1‑90 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcac      <  2:721822/89‑1 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcac      >  1:721822/1‑89 (MQ=255)
                                                                                 cGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAAt   >  1:570104/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAAc                                     <  2:147713/55‑1 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAAc                                     >  1:147713/1‑55 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcgc                      >  1:309936/1‑70 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcgc                      <  2:309936/70‑1 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  2:213286/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  <  2:461601/90‑1 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  1:85262/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  2:33521/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  1:27466/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  1:225466/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  2:711233/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  2:302990/1‑90 (MQ=255)
                                                                                         |                                                                                   
TAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG  >  NZ_CP009273/3898735‑3898907

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