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	96.8%	80.8 / NA	31	noncoding (855/1542 nt)	rrsH	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (17/13);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCAGGTGCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGC  >  NC_000913/224525‑224870                                                                                                     |                                                                                                                                                                                                                                                       gCTCAGGTGCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGa                                                                                                  >  1:348866/1‑250 (MQ=14)           gAAAGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGATTAATTCGATGCAA‑CGCGAGGAACCTTACCCGGGCTTGACATCCACAAAACTTTcccagg                                                                                       >  2:157233/1‑247 (MQ=11)            aaaGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAgaga                                                                                      <  1:78173/251‑1 (MQ=21)             aagggTGGGGAGCAAACAGGATTAGATACCCGGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGTGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGAt                                                                                     <  1:519217/247‑1 (MQ=255)              aggtgggggAGAAAACAGGATTAGATACCCTGGTAGTCTGCGACGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTCTCCAGAGATg                                                                                    <  2:350972/246‑1 (MQ=11)                  tGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCCCGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGAtt                                                                                >  1:215952/1‑251 (MQ=37)                  tGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACATGGAGGTTGTACCCTTGCGGCGTGGCTACCGGAGCTAAAGCGTTAAGTCGACCGCCTGGGGAGGACGGCCGCCAGGTTAAAACTCAAATGAATTGACGGAGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑AGGGCAGAACCTTCCCTGGTCTTGAGATCCACAGACATTTCCAGAGATGGAtt                                                                                >  2:395845/1‑251 (MQ=11)                       aGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGc                                                                           >  1:10395/1‑251 (MQ=39)                          aaaCAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCtt                                                                        >  1:137904/1‑251 (MQ=255)                          aaaCAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCtt                                                                        <  1:152767/251‑1 (MQ=255)                           aaCAGGATTAGATACCCGGGTAGGCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTGGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCTCAGAACTTTCCAGAGATGGATTGGTGCCTTc                                                                       <  1:42731/251‑1 (MQ=255)                                gATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGa                                                                   >  1:152371/1‑250 (MQ=255)                                      aTACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCCAATGAATTGACGGGGGCCCGCACAAGCGGGGGGGCCTGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCCAAGAACTTTCAAGAGGTGGGTTGGTGCCTTAGGGCACtgggg                                                            >  2:495554/1‑248 (MQ=11)                                       tACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACtgtg                                                             >  1:29824/1‑249 (MQ=255)                                        aCCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGCGCTAACGCGTTAAGGCGACCGCCTGGGGAGTCCGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGGGGTTTAATTCGTTGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGACCTTTCCAGAGATGGACTGGTGCCTTCGGGAACTGTga                                                            >  2:141448/1‑249 (MQ=14)                                             ggTAGGCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAgg                                                      <  1:532425/250‑1 (MQ=255)                                                 gTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGtgctg                                                 >  1:481612/1‑251 (MQ=255)                                                    cACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCAt                                              >  2:560396/1‑251 (MQ=255)                                                    cACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCGTGAGGCGGGGCTACCGGATCTAACGCGTTAAGTCAACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCAt                                              <  2:29824/251‑1 (MQ=255)                                                      cGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATgg                                            >  1:71087/1‑251 (MQ=255)                                                      cGCCGTAAACGATGTCGACTTGGAGGTGTTGCCCTTGAGGCGTGGTTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGATCTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATgg                                            <  2:519808/251‑1 (MQ=255)                                                            aaaCTATGTCGACTTGGAGTTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTTATACAGGTTCTGCAtgtctgtc                                       <  2:137904/250‑1 (MQ=14)                                                             aaCGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCgg                                                                     <  1:200616/219‑1 (MQ=255)                                                             aaCGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGAGGAGCATGTGGGTTAATTCGATGCCA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATAGGTGCCTTCgg                                                                     >  2:200616/1‑219 (MQ=255)                                                                gATGTCGACGTGGAGGTTGTGCCCGTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGGCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGGGAGACAGGTTCTGCATGGCTGTCGTCAg                                  <  2:394209/251‑1 (MQ=255)                                                                        cgtggtggTTGGGCCCTGGAGGGGGGTCTTCCGGAGCTAACGCGTTAAGGCGACCGCCTTGGTATTACGGCCGCAAGGTTAAAACTCAAATGAAGTGTCGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCTCAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGtgtt                          <  2:486874/245‑1 (MQ=14)                                                                             aGGTTGTGCCCTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGAGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGa                      >  1:251709/1‑250 (MQ=255)                                                                                 ggtgCCCTGGGGGGGTGGCTACGGGAGCGAACGCGTGAAGTCGACCGGCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATg                  <  2:170793/249‑1 (MQ=14)                                                                                      ccTTGAGGCGTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATTCATTGACGGGGGCCCGCACAAGCGGCGGCGCATGTGGTTTAATTCGA‑GCAAACGCGAAGAACCTTACCTGGTCTTGACCTCCACAGAACTTTCCTGAGATGGTTTGGAGCCTTCGGGGACTGTGAAACGGGTGCTGCATGGCTGTAGTCAGCTCGTGTTGTGAGATGTTGGGt            >  2:400204/1‑251 (MQ=14)                                                                                               gTGGCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTCATTCGCTGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCg   >  1:334823/1‑251 (MQ=255)                                                                                                 ggCTACCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGc  <  1:261961/250‑1 (MQ=255)                                                                                                     |                                                                                                                                                                                                                                                       GCTCAGGTGCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAA‑CGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGC  >  NC_000913/224525‑224870

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