Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	224,625	T→A	noncoding (855/1542 nt)	rrsH →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	224,625	0	T	A	92.6%	62.2 / ‑6.5	27	noncoding (855/1542 nt)	rrsH	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (16/9);  minor base G (0/2);  total (16/11)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCA  >  NC_000913/224535‑224871                                                                                           |                                                                                                                                                                                                                                                       gAAAGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGAAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAAGGCGTTAAGTCGACCAACTGGGGAGTACGGCCGCAAGGGTAAAACTCAAATGAATGGACGGGGGACCGCACAAACGGTGGAGCAAGTGGGTTACTTCGATGCAACGCGAAGAACATTACCTGGACTTGAAATCAACAGAACTTTccca                                                                                          >  1:352849/1‑247 (MQ=11)   aaGCGTGGGGAGCAAACAGGCTTCGATACCCTGGTGGTCCACGCCGTAAACGATGTCGCCTTGGCGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCACCGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGCGAt                                                                                      >  2:215399/1‑251 (MQ=11)   aaGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGAt                                                                                      <  1:128850/251‑1 (MQ=21)      cgtagggAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTCCGGCCGCAAGGTTAAAACTCACATGAATTGACGGGGGCCCGCACCAGCGGTGGAGCCTGTGGTTTAATTCGATGCAACGCGAAGACCCTTACATGGTCTTGACACCCACAGAACTTTCCAgagaaggt                                                                                   >  2:448819/5‑247 (MQ=11)          gggAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGCCTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCCAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTgg                                                                               >  1:56005/1‑251 (MQ=255)                aacaaggatataaTACCCTGGTATTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCTTGGCTGCCGGATCTAACGCGTTAAGTCGACCGCCTGGTGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACTGTGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGATCCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCtt                                                                         <  2:587094/239‑1 (MQ=14)                aaaCAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCGGAGCTGGATTGGTGCCtt                                                                         >  2:405393/1‑251 (MQ=255)                  acagtataaGATACCCTGGTCTCCCACGCCGTAAACGATGTCGACTTGGAGTTTGTGCCCTTGAGGCGTGGCTACCGGATCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCGTCg                                                                       <  2:26238/243‑1 (MQ=14)                   cAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTATGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCgg                                                                      <  1:555795/251‑1 (MQ=255)                         tAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGACCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACGTTCCAGAGATGGATTGGTGCCTTCGGGAACtg                                                                >  2:439896/1‑251 (MQ=255)                          aGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACtgt                                                               >  1:457257/1‑251 (MQ=255)                               cccTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACATGGTCTTGACATCCACAGAACTTTCCACAGATGGATTGGAGCCTTCGGGAACtgggcgac                                                          >  2:9719/1‑245 (MQ=255)                                ccTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTCCGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGAc                                                          >  1:123455/1‑250 (MQ=255)                                  tGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCaga                                                                                         <  1:13322/217‑1 (MQ=25)                                  tGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAAACTTACCTGGTCTTGACATCCACAGAACTTTCCaga                                                                                         >  2:13322/1‑217 (MQ=25)                                      aGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAACTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGGCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGGGAGAAAGGtgct                                                   >  1:71888/1‑251 (MQ=255)                                         ccACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGtgctgc                                                 >  1:489010/1‑250 (MQ=255)                                             gCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGATATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGc                                            <  1:66285/251‑1 (MQ=255)                                                    aCGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTgtcgtc                                     <  2:119808/251‑1 (MQ=255)                                                        tGTCGACTTGGCGGTTGTGCCCTTGAGGCGTGGCTAACGGAGCTCACGCGTTAAGTCGACCGCCTGGGGAGTCCGGCCGACAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCAGCATGGCTGTCGTCacct                                 >  2:28701/1‑248 (MQ=14)                                                              ctttgAGGTTGTGCCCTTGAGGCGTGGCCACCGCAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGCCGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGCTGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGCCAGGTGCTGCATGGCTGTCGTCAGCTCGtgtt                           >  2:97028/5‑251 (MQ=14)                                                              cTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGTGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGGAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGGGATGGATTGGTGCCTGCTGGAACTGTGAGACAGGTGCGGACTGGCTGTCGTCAGCTCGtgtt                           >  2:98414/1‑251 (MQ=14)                                                                  gagtgtgtgCCCGTGAGGCGTGGCTGCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACGCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAGCCTTACCAGGTCTTGACATCCACAGAACTGTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGGCGTCAGCTCGTGTtgtgg                       <  2:389993/246‑2 (MQ=14)                                                                          gCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTg                <  2:109604/250‑1 (MQ=255)                                                                            ccTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGt             >  1:457587/1‑251 (MQ=255)                                                                               ggAGGCGTGGCTACCGGAGCGAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGGGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGtt            <  2:316334/248‑1 (MQ=255)                                                                                       ggCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCa  <  1:254518/251‑1 (MQ=255)                                                                                           |                                                                                                                                                                                                                                                       GAAAGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCA  >  NC_000913/224535‑224871

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