Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 699,398 A→G 25.0% intergenic (+50/‑36) PP_16SD → / → PP_t05 16S ribosomal RNA/tRNA‑Ile

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947699,3980AG25.0% 23.3 / 4.9 16intergenic (+50/‑36)PP_16SD/PP_t0516S ribosomal RNA/tRNA‑Ile
Reads supporting (aligned to +/- strand):  ref base A (5/7);  new base G (2/2);  total (7/9)
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

GGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATCTGCCCAGACCTACCATTACATGGTTCAGCCGTAGAATACGG  >  NC_002947/699293‑699536
                                                                                                         |                                                                                                                                          
gggTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGTCTGTa                                                                                                >  2:316254/1‑150 (MQ=32)
         cGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGTCTGTAGCTCAGTTg                                                                                       <  1:123687/150‑1 (MQ=32)
           tAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGTCTGTAGCTCAGTTGGt                                                                                     <  1:316254/150‑1 (MQ=32)
               aaGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTACCGCTACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                                                                 >  1:54963/1‑150 (MQ=17)
                       cGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCAccc                                                                         >  1:317635/1‑150 (MQ=32)
                               aaCCTGCGGCTGGATCACCTCCTTAATCGGCAACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAg                                                                 <  2:266440/150‑1 (MQ=17)
                                                 ctcctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCa                                               >  2:576731/1‑150 (MQ=32)
                                                   cctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCaaa                                             <  1:148721/150‑1 (MQ=32)
                                                   cctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCATGCTAAGGGTGAGGTCGGCAGTTCaaa                                             <  2:317635/150‑1 (MQ=17)
                                                            acgacATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATCTGCCCAg                                    <  2:430563/150‑1 (MQ=32)
                                                                      ccTGCTGATGAGCTCCCACACGAATTGCTTGATTCGTTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGTCTGTAGCACAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATCTGCCCAGACCTACCAtt                          >  1:38878/1‑150 (MQ=17)
                                                                      ccTGCTGATGAGCTCCCACACGAATTGCTTGATTCGTTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGACTGTTGAACTGTTGGTTAGAGCGCACCCCTGATACGGGTGAGGTCGGCAGTTCCAATCTGACCAGAGCTACCAtt                          >  1:38872/1‑150 (MQ=11)
                                                                           tgatgaGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAg                                                           <  1:538254/112‑1 (MQ=32)
                                                                           tgatgaGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAg                                                           >  2:538254/1‑112 (MQ=33)
                                                                                              ttGCTTGATTCGTTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATCTGCCCAGACCTACCATTACATGGTTCAGCCGTAGAATACgg  <  2:38878/150‑1 (MQ=17)
                                                                                              tgtcTTGAGTCGTTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGATCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGTTCGGCAGTTCAAATCTGACCAGACCTACCCTTAAATGGTTCAGCCGTAGAATACgg  <  2:38872/147‑1 (MQ=11)
                                                                                                         |                                                                                                                                          
GGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCTACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATCTGCCCAGACCTACCATTACATGGTTCAGCCGTAGAATACGG  >  NC_002947/699293‑699536

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