Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,031,888 A→T 19.4% K229M (AAG→ATG) ‡ dppA‑IV → dipeptide transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,031,8880AT19.4% 66.1 / 12.6 31K229M (AAG→ATG) ‡dppA‑IVdipeptide transport protein
Reads supporting (aligned to +/- strand):  ref base A (13/12);  new base T (3/3);  total (16/15)
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

CGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGTGCCAGGTTGCCCTGTCGCCCAAGCCGTTGGACATTGCCGAAGCGGGCCAGGACGGTAACCTCAAGG  >  NC_002947/6031742‑6032031
                                                                                                                                                  |                                                                                                                                               
cGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGcc                                                                                                                                              >  1:359820/1‑150 (MQ=255)
       tACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCTTGCcggcgaag                                                                                                                                       <  2:179357/150‑4 (MQ=255)
          gCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGAc                                                                                                                                    <  1:79309/150‑1 (MQ=255)
                    tgctCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCt                                                                                                                          >  2:65001/1‑150 (MQ=255)
                      ctCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTc                                                                                                                        <  2:392407/150‑1 (MQ=255)
                       tCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCAAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGAGCCGCTGATCTTCg                                                                                                                       >  2:484257/1‑150 (MQ=255)
                                  acacCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCTTGCcggcgaa                                                                                                                                        >  2:368358/1‑120 (MQ=255)
                                  acacCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCTTGCcggcgaa                                                                                                                                        <  1:368358/122‑3 (MQ=255)
                                        gAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTAACGCGCCAACCCGGACTACATCGCCGGCAAGCCGGCGGTCGAc                                                                                                                                    <  2:293675/120‑1 (MQ=255)
                                        gAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTAACGCGCCAACCCGGACTACATCGCCGGCAAGCCGGCGGTCGAc                                                                                                                                    >  1:293675/1‑120 (MQ=255)
                                                            gATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCTTGCcggcgaag                                                                                                                                       >  1:453648/1‑94 (MQ=39)
                                                            gATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCTTGCcggcgaag                                                                                                                                       <  2:453648/97‑4 (MQ=39)
                                                             aTCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTc                                                                                 >  1:188363/1‑150 (MQ=255)
                                                                               gtgtTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCTAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGAGCCAg                                                               >  1:71177/1‑150 (MQ=255)
                                                                               gtgtTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGTGCCAg                                                               <  1:65001/150‑1 (MQ=255)
                                                                                     cAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGTGCCAGGTTGcc                                                         <  2:295987/150‑1 (MQ=255)
                                                                                       gCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGTGCCAGGTTGCCCt                                                       <  2:359820/150‑1 (MQ=255)
                                                                                               aGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGc                                                                                                       >  2:154240/1‑94 (MQ=255)
                                                                                               aGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGc                                                                                                       <  1:154240/94‑1 (MQ=255)
                                                                                                       gCCGTGGTGCGCTACTACGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGTGCCAGGTTGCCCTGTCGCCCAAGCCGTTg                                       >  2:45309/1‑150 (MQ=255)
                                                                                                       gCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTc                                                                                                                        <  1:83364/69‑1 (MQ=255)
                                                                                                       gCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTc                                                                                                                        >  2:83364/1‑69 (MQ=255)
                                                                                                               gcgcTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGTGCCAGGTTGCCCTGTCGCCCAAGCCGTTGGACATTGc                               <  2:66144/150‑1 (MQ=255)
                                                                                                                 gcTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGTGCCAGGTTGCCCTGTCGCCCAAGCCGTTg                                       >  2:289558/1‑140 (MQ=255)
                                                                                                                 gcTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGTGCCAGGTTGCCCTGTCGCCCAAGCCGTTg                                       <  1:289558/140‑1 (MQ=255)
                                                                                                                   tACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGTGCCAGGTTGCCCTGTCGCCCAAGCCGTTGGACATTGCCGaa                           >  2:236099/1‑150 (MQ=255)
                                                                                                                               ccGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGTGCCAGGTTGCCCTGTCGCCCAAGCCGTTGGACATTGCCGAAGCGGGCCAGGAc               >  1:249034/1‑150 (MQ=255)
                                                                                                                                cGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGTGCCAGGTTGCCCTGTCGCCCAAGCCGTTGGACATTGCCGAAGCGGGCCAGGACg              <  1:236099/150‑1 (MQ=255)
                                                                                                                                       ctTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGTGCCAGGTTGCCCTGTCGCCCAAGCCGTTGGACATTGCCGAAGCGGGCCAGGACGGTAACCt       >  1:54165/1‑150 (MQ=255)
                                                                                                                                          cGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGTGCCAGGTTGCCCTGTCGCCCAAGCCGTTGGACATTGCCGAAGCGGGCCAGGACGGTAACCTCaa    <  2:54165/150‑1 (MQ=255)
                                                                                                                                            ccGGCTTGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGTGCCAGGTTGCCCTGTCGCCCAAGCCGTTGGACATTGCCGAAGCGGGCCAGGACGGTAACCTCAAgg  >  2:108477/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CGCAGAATACGCTGACAAGCTGCTCAAGGCCGGCACACCGGAAAAGCTCAACAACCAACCGATCGGCACCGGGCCGTTCGTGTTCCAGCGTTTCCAGAAAGACGCCGTGGTGCGCTACCGCGCCAACCCGGACTACTTCGCCGGCAAGCCGGCGGTCGACCCGCTGATCTTCGCCATCACCACCGATGCCAACGTGCGCCTGCAGAAGCTCAAGCGTGGCGAGTGCCAGGTTGCCCTGTCGCCCAAGCCGTTGGACATTGCCGAAGCGGGCCAGGACGGTAACCTCAAGG  >  NC_002947/6031742‑6032031

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