Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 705,917 G→A 19.3% intergenic (+27/‑176) 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,9170GA19.3% 67.7 / 12.1 31intergenic (+27/‑176)nagE/glnSN‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (12/13);  new base A (3/3);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

CGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTTA  >  NC_000913/705770‑706063
                                                                                                                                                   |                                                                                                                                                  
cGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGa                                                                                                                                                   <  4:120656/149‑1 (MQ=255)
  cccGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGc                                                                                                                                                 <  1:129856/149‑1 (MQ=255)
     gCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                              <  1:552108/149‑1 (MQ=255)
                                     atCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                              >  1:658202/1‑117 (MQ=255)
                                     atCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                              <  2:658202/117‑1 (MQ=255)
                                      tCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGc                                                                                                                                                 >  2:980083/1‑113 (MQ=255)
                                      tCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGc                                                                                                                                                 <  1:980083/113‑1 (MQ=255)
                                                     gCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGc                                                                                                                                                 <  1:458307/98‑1 (MQ=255)
                                                     gCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGc                                                                                                                                                 >  2:458307/1‑98 (MQ=255)
                                                              ttAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                              >  1:1341789/1‑92 (MQ=255)
                                                              ttAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                              <  2:1341789/92‑1 (MQ=255)
                                                                aaaGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                              <  4:1349301/90‑1 (MQ=255)
                                                                aaaGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                              >  3:1349301/1‑90 (MQ=255)
                                                                aaaGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTCAAGCGTc                                                                                                                                              >  1:850178/1‑90 (MQ=255)
                                                                aaaGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTCAAGCGTc                                                                                                                                              <  2:850178/90‑1 (MQ=255)
                                                                         ggCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTtctc                                                                                                            >  4:305984/1‑115 (MQ=255)
                                                                         ggCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTtctc                                                                                                            <  3:305984/115‑1 (MQ=255)
                                                                             atatTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGAtc                                                                                                                 >  2:12155/1‑106 (MQ=255)
                                                                             atatTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGAtc                                                                                                                 <  1:12155/106‑1 (MQ=255)
                                                                                          tCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                              >  4:572319/1‑64 (MQ=255)
                                                                                          tCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                              <  3:572319/64‑1 (MQ=255)
                                                                                                              tCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCtt                                                                      <  2:819111/116‑1 (MQ=255)
                                                                                                              tCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCtt                                                                      >  1:819111/1‑116 (MQ=255)
                                                                                                                aaaaaGTAATCTGCTTTATGCCTGATGCGACGCTCAAGCGTCGCATc                                                                                                                                         >  2:208713/1‑47 (MQ=21)
                                                                                                                aaaaaGTAATCTGCTTTATGCCTGATGCGACGCTCAAGCGTCGCATc                                                                                                                                         <  1:208713/47‑1 (MQ=25)
                                                                                                                    aGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGcc                               >  3:800288/1‑149 (MQ=255)
                                                                                                                                     cTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGAt              <  4:800288/149‑1 (MQ=255)
                                                                                                                                            gACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaaaaa                                                >  3:670751/1‑108 (MQ=255)
                                                                                                                                            gACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaaaaa                                                <  4:670751/108‑1 (MQ=255)
                                                                                                                                                 ttGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTta  >  4:499514/1‑149 (MQ=255)
                                                                                                                                                 ttGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTta  >  4:761421/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTTA  >  NC_000913/705770‑706063

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