Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 705,916 T→C 20.7% 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,9160TC20.7% 63.2 / 11.8 29intergenic (+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 (10/13);  new base C (3/3);  total (13/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 = 4.14e-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                                                                                                                                              <  2:658202/117‑1 (MQ=255)
                                     atCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                              >  1:658202/1‑117 (MQ=255)
                                      tCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGc                                                                                                                                                 <  1:980083/113‑1 (MQ=255)
                                      tCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGc                                                                                                                                                 >  2:980083/1‑113 (MQ=255)
                                        gacgaTTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCtt                                                                                                                                                     <  1:719359/107‑1 (MQ=255)
                                        gacgaTTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCtt                                                                                                                                                     >  2:719359/1‑107 (MQ=255)
                                                     gCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGc                                                                                                                                                 <  1:458307/98‑1 (MQ=255)
                                                     gCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGc                                                                                                                                                 >  2:458307/1‑98 (MQ=255)
                                                       ttGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCtt                                                                                                                                                     <  1:418345/92‑1 (MQ=255)
                                                       ttGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCtt                                                                                                                                                     >  2:418345/1‑92 (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                                                                                                                                              <  2:850178/90‑1 (MQ=255)
                                                                aaaGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTCAAGCGTc                                                                                                                                              >  1:850178/1‑90 (MQ=255)
                                                                         ggCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTtctc                                                                                                            <  3:305984/115‑1 (MQ=255)
                                                                         ggCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTtctc                                                                                                            >  4:305984/1‑115 (MQ=255)
                                                                             atatTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGAtc                                                                                                                 >  2:12155/1‑106 (MQ=255)
                                                                             atatTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGAtc                                                                                                                 <  1:12155/106‑1 (MQ=255)
                                                                                          tCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                              <  3:572319/64‑1 (MQ=255)
                                                                                          tCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                              >  4:572319/1‑64 (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 < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: