Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,167,122 A→G N377S (AAC→AGC)  gtsD → mannose/glucose ABC transporter ATP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,167,1220AG93.8% 48.2 / ‑3.8 16N377S (AAC→AGC) gtsDmannose/glucose ABC transporter ATP binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (8/7);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GGCGCCAGATGTGGCGTGCCGGGTGGGTGACACCCTCAACCTGCAATTCGACCCGGCCCGGGTGCTGCTGTTCGACGCCGCCAATGGCGAGCGCCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAACGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCATCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGGCAGC  >  NC_002947/1166991‑1167250
                                                                                                                                   |                                                                                                                                
ggCGCCAGATGTGGCGTGCCGGGTGGGTGACACCCTCAACCTGCAATTCGACCCGGCCCGGGTGCTGCTGTTCGACGCCGCCAATGGCGAGCGCCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAGCGTg                                                                                                                              <  2:167703/136‑1 (MQ=255)
 gcgcCAGATGTGGCGTGCCGGGTGGGTGACACCCTCAACCTGCAATTCCACCCGGCCCGGGTGCTGCTGTTCGACGCCGCCAATGGCGAGCGCCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAGCGTgg                                                                                                                             <  2:264346/136‑1 (MQ=255)
                tGCCGGGTGGGTGACACCCTCAACCTGCAATTCGACCCGGCCCGGGTGCTGCTGTTCGACGCCGCCAATGGCGAGCGCCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAGCGTGGCTCACTTCAAAGGcc                                                                                                              >  2:257670/1‑136 (MQ=255)
                  ccGGGTGGGTGACACCCTCAACCTGCAATTCGACCCGGCCCGGGTGCTGCTGTTCGACGCCGCCAATGGCGAGCGCCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAGCGTGGCTCACTTCAAAGGCCGt                                                                                                            <  2:228728/136‑1 (MQ=255)
                                    cAACCTGCAATTCGACCCGGCCCGGGTGCTGCTGTTCGACGCCGCCAATGGCGCGCGCCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAGCGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATc                                                                                          >  1:100810/1‑136 (MQ=255)
                                                      ggcccggGTGCTGCTGTTCGACGCCGCCAATGGCGAGCGCCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAGCGTGGCTCACTTCAAAGGCCGt                                                                                                            <  1:47997/100‑1 (MQ=255)
                                                      ggcccggGTGCTGCTGTTCGACGCCGCCAATGGCGAGCGCCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAGCGTGGCTCACTTCAAAGGCCGt                                                                                                            >  2:47997/1‑100 (MQ=255)
                                                      ggcccggGTGCTGCTGTTCGACGCCGCCAATGGCGAGCGCCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAGCGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACg                                                                        <  1:264562/136‑1 (MQ=255)
                                                           gggTGCTGCTGTTCGACGCCGCCAATGGCGAGCGCCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAGCGTGGCGCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGCCGCAAAg                                                                   >  1:21836/1‑136 (MQ=255)
                                                                    tgTTCGACGCCGCCAATGGCGAGCGCCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAGCGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAAc                                                          >  2:277327/1‑136 (MQ=255)
                                                                                                           gTACTGGCGCAACCGCAAAGGACAGCGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCATCAAGACACTGGGATCGTTGGCcttgc                   <  2:100810/136‑1 (MQ=255)
                                                                                                               tGGCGCAACCGCAAAGGACAGCGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCATCAAGACACTGGGATCGTTGGCCttgcttgc               >  2:154864/1‑136 (MQ=255)
                                                                                                                    cAACCGCAAAGGACAGCGTGGCTCACTTCAAAGGCCGTTGACTCGTCt                                                                                                  >  1:186782/1‑48 (MQ=255)
                                                                                                                    cAACCGCAAAGGACAGCGTGGCTCACTTCAAAGGCCGTTGACTCGTCt                                                                                                  <  2:186782/48‑1 (MQ=255)
                                                                                                                         gCAAAGGACAGCGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCATCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGgc     >  2:131316/1‑136 (MQ=255)
                                                                                                                            aaGGACAACGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCATCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGgcagc  <  1:192341/136‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                
GGCGCCAGATGTGGCGTGCCGGGTGGGTGACACCCTCAACCTGCAATTCGACCCGGCCCGGGTGCTGCTGTTCGACGCCGCCAATGGCGAGCGCCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAACGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCATCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGGCAGC  >  NC_002947/1166991‑1167250

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: