Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 705,916 T→C 11.8% intergenic (+26/‑177) nagE → / → glnS N‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components/glutamyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913705,9160TC11.8% 121.8 / 5.2 42intergenic (+26/‑177)nagE/glnSN‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (17/20);  new base C (2/3);  total (19/23)
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

CGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTT  >  NC_000913/705770‑706062
                                                                                                                                                  |                                                                                                                                                  
cGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGa                                                                                                                                                  <  1:603702/149‑1 (MQ=255)
  cccGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGc                                                                                                                                                <  1:1371159/149‑1 (MQ=255)
  cccGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGc                                                                                                                                                <  1:235533/149‑1 (MQ=255)
     gCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                             <  4:583149/149‑1 (MQ=255)
         gatgatTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTTATGCGACGCTTGAGCGTCGCAt                                                                                                                                         >  2:1186963/1‑149 (MQ=255)
                                gcaATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGt                                                                                                                                              >  3:1751060/1‑121 (MQ=255)
                                gcaATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGt                                                                                                                                              <  4:1751060/121‑1 (MQ=255)
                                     atCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGAtcttc                                                                                                             <  1:1186963/149‑1 (MQ=255)
                                             tttCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCGGATGCGACGCTCAAGCGTCGc                                                                                                                                           >  4:1308715/1‑111 (MQ=255)
                                             tttCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCGGATGCGACGCTCAAGCGTCGc                                                                                                                                           <  3:1308715/111‑1 (MQ=255)
                                                         gATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCtt                                                                                                   <  3:1373874/139‑1 (MQ=255)
                                                         gATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTtcat                                                                                         >  2:868852/1‑149 (MQ=255)
                                                         gATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTACTCTGCCGCtt                                                                                                   >  4:1373874/1‑139 (MQ=255)
                                                              ttAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGc                                                                                                                                                <  2:19405/89‑1 (MQ=255)
                                                              ttAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGc                                                                                                                                                >  1:19405/1‑89 (MQ=255)
                                                               tAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCatcaatca                                                                                   <  1:868852/149‑1 (MQ=255)
                                                                         ggCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGGTGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATTAATCATCCCCATAAt                                                                         >  3:1191257/1‑149 (MQ=255)
                                                                             atatTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                             <  3:706417/77‑1 (MQ=255)
                                                                             atatTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCtt                                                                     >  4:521760/1‑149 (MQ=255)
                                                                             atatTGTGGCGGGTCAAACACCGCTGTATGAAAACAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                             >  4:706417/1‑77 (MQ=255)
                                                                                    ggCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGAt                                                              <  3:521760/149‑1 (MQ=255)
                                                                                             aaCACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCtt                                                                     >  4:1447028/1‑133 (MQ=255)
                                                                                             aaCACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCtt                                                                     <  3:1447028/133‑1 (MQ=255)
                                                                                                acCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                             >  1:404091/1‑58 (MQ=255)
                                                                                                acCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                             <  2:404091/58‑1 (MQ=255)
                                                                                                acCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                             >  4:1674081/1‑58 (MQ=255)
                                                                                                acCGCTGTATGAAATCAAAAAGTAATCTGCCTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                             <  3:1674081/58‑1 (MQ=255)
                                                                                                   gCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTa                                                                 <  4:889139/131‑1 (MQ=255)
                                                                                                   gCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTa                                                                 >  3:889139/1‑131 (MQ=255)
                                                                                                           aaaTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGCAGAGAAGATCTCCACCGCttg                                                                                                  <  4:20505/90‑2 (MQ=37)
                                                                                                                aaaaaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTaa                                          <  2:528411/141‑1 (MQ=255)
                                                                                                                aaaaaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTaa                                          >  1:528411/1‑141 (MQ=255)
                                                                                                                         tCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTctct                                                                                                          <  2:990147/68‑1 (MQ=255)
                                                                                                                         tCTGCTTTATGACTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTctct                                                                                                          >  1:990147/1‑68 (MQ=255)
                                                                                                                                    tcTGATTCGACGCTCAAGCGTCGCATCAAACAATGACAAGCGGTGGAGATCTTCCCTGCCGTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGa              <  1:651406/148‑1 (MQ=255)
                                                                                                                                    ccTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTCCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTCAGATTATCAATTTTAAAAAACTAACAGTTGTCCGCCTGTCCCGCTTATAAGa              <  4:1191257/149‑1 (MQ=255)
                                                                                                                                            gACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGt                                                                   >  4:308891/1‑88 (MQ=255)
                                                                                                                                            gACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGt                                                                   <  3:308891/88‑1 (MQ=255)
                                                                                                                                            gACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGc      >  4:1518235/1‑149 (MQ=255)
                                                                                                                                            gACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAAc                                         >  1:199331/1‑114 (MQ=255)
                                                                                                                                            gACGCTCAAGCGTCGAATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAAc                                         <  2:199331/114‑1 (MQ=255)
                                                                                                                                                cTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGtt  >  3:875093/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTT  >  NC_000913/705770‑706062

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: