Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	698,371	C→T	20.0%	noncoding (548/1522 nt)	PP_16SD →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	698,371	0	C	T	20.0%	36.7 / 3.8	20	noncoding (548/1522 nt)	PP_16SD	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (7/9);  new base T (2/2);  total (9/11)
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

TGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGGATGTGAAAGCCCCGG  >  NC_002947/698229‑698418                                                                                                                                               |                                                tGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGGGCAAGCGTTAATCGGAATTACTGGGCGTaacg                                           >  2:272905/1‑147 (MQ=17)  gTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGTGTAAAgc                                          >  6:301570/1‑149 (MQ=17)  gTAAAGCACTTTAAGTTGGGAGCAAGCGCAGTAAGTTAATACCTTGCAGTCTTGAAGTAACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGAGCAAGCGTTAATCGGAATTACTGGGTGTaaacc                                          >  6:296729/1‑147 (MQ=11)   tAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAgcg                                         >  8:222904/1‑149 (MQ=32)   tAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAgcg                                         >  3:148529/1‑149 (MQ=32)     aaGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCTGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAgcgc                                        <  4:131163/148‑1 (MQ=17)           tttAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAgg                                 >  5:295958/1‑149 (MQ=32)           tttAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAgg                                 >  1:312949/1‑149 (MQ=32)             tAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAgg                                 <  6:295958/147‑1 (MQ=32)               aGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAggtgg                              <  3:314711/148‑1 (MQ=32)                gTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGtt                            >  6:199873/1‑149 (MQ=32)                      aggaAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTaa                      <  5:16223/149‑1 (MQ=32)                      aggaAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTaa                      <  3:291325/149‑1 (MQ=32)                         aaGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGtt                   <  2:172881/149‑1 (MQ=32)                         aaGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGtt                   <  7:222904/149‑1 (MQ=32)                           gggCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTgg                 <  1:66017/149‑1 (MQ=32)                            ggCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGGa                <  7:106632/149‑1 (MQ=32)                             gCAGTAAGTTAATACGTTGCTGTCTTGACGTTAGAGAGAGAATAAGCACCGGCTAACTCTGTGCCAGGAGCCGCGGTTATTCAGAGGGTGCAAGCGTTAATCGGTATTACTGGGTGTAAAGCTCGCGTAGGTGGTTTGTTAAGTTGGaa               <  5:296729/149‑2 (MQ=11)                             gCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGTGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGGaa               <  5:301570/149‑2 (MQ=17)                                          aCCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGGATGTGAAAGCCCCgg  >  4:161774/1‑149 (MQ=32)                                                                                                                                               |                                                TGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGGATGTGAAAGCCCCGG  >  NC_002947/698229‑698418

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