Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	223,855	T→A	30.0%	noncoding (85/1542 nt)	rrsH →	16S ribosomal RNA of rrnH operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	223,855	0	T	A	30.0%	17.5 / 15.7	20	noncoding (85/1542 nt)	rrsH	16S ribosomal RNA of rrnH operon
Reads supporting (aligned to +/- strand):  ref base T (8/6);  new base A (2/4);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-02

CAAGTCTCAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGT  >  NC_000913/223709‑223888                                                                                                                                                   |                                  cAAGTCTCAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAgcttgc                                 >  1:908503/1‑149 (MQ=32)       tcAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCAAGcttcttc                           <  3:1183685/149‑2 (MQ=12)         aaGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGc                         >  2:1270174/1‑149 (MQ=32)         aaGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCAAGcttcttcct                         >  4:1127163/1‑146 (MQ=255)         aaGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCAAGcttcttcct                         <  3:1401339/149‑4 (MQ=255)         aaGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCAAGcttcttcct                         <  3:1127163/149‑4 (MQ=255)                aaCACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAg                  >  2:1189893/1‑149 (MQ=32)                    cGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGc              >  2:722254/1‑149 (MQ=32)                     gTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGc              <  1:1294054/148‑1 (MQ=32)                     gTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGc              >  2:1294054/1‑148 (MQ=32)                     gTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCg             <  1:1270174/149‑1 (MQ=32)                       aaTTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGa           <  1:547429/149‑1 (MQ=32)                        aTTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGAc          >  2:96047/1‑149 (MQ=32)                            aTTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGt      <  1:722254/149‑1 (MQ=32)                               aCAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCAAGcttctt                            <  3:1133602/124‑1 (MQ=12)                               aCAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCAAGcttctt                            >  4:1133602/1‑124 (MQ=12)                                cAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTg                <  3:769341/135‑1 (MQ=31)                                cAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTg                >  4:769341/1‑135 (MQ=31)                                cAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAACCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGt  >  3:287913/1‑149 (MQ=14)                                 aaaGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAATACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGc              <  2:1462219/136‑1 (MQ=11)                                 aaaGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAATACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGc              >  1:1462219/1‑136 (MQ=11)                                                                                                                                                   |                                  CAAGTCTCAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGT  >  NC_000913/223709‑223888

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