Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	226,403	0	C	T	70.2%	21.5 / 40.4	47	noncoding (645/2904 nt)	rrlH	23S ribosomal RNA of rrnH operon
Reads supporting (aligned to +/- strand):  ref base C (6/8);  new base T (17/16);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01

GAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTGATCTAGC  >  NC_000913/226258‑226449                                                                                                                                                  |                                               gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                           <  2:83049/150‑1 (MQ=2) gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                           <  1:498876/150‑1 (MQ=2) gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                           <  1:337310/150‑1 (MQ=2) gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                           <  2:300204/150‑1 (MQ=2) gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                           <  1:299368/150‑1 (MQ=2) gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCg                                          <  2:217682/151‑1 (MQ=2) gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCg                                          <  1:159190/151‑1 (MQ=2) gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTACATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                           <  1:252778/150‑1 (MQ=2)  aaataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCga                                         <  2:290311/151‑2 (MQ=2)  aaataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAATGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCg                                          <  1:189694/150‑1 (MQ=2)   aataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                           <  1:113336/145‑1 (MQ=2)   aataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                           >  2:113320/4‑148 (MQ=2)   aataGAACCTGAAACCGTATGCGTACAAGCAGTAGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                           <  1:203396/145‑1 (MQ=2)   aataGAACCTGAAACCGTATGCGTACAAGCAGTAGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                           >  2:203362/4‑148 (MQ=2)    ataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCgagt                                       <  2:36496/149‑4 (MQ=2)    ataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGttt                                       >  1:88249/3‑150 (MQ=2)     taGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCTAAAGCGAGTCTTAATAGGGCGttt                                       >  1:213761/2‑149 (MQ=2)      aGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGttt                                       >  1:717664/1‑148 (MQ=2)       gAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAg                                     <  1:441554/151‑1 (MQ=31)          ccTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTg                                  >  1:126576/1‑151 (MQ=31)           cTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGc                                 >  2:710348/1‑151 (MQ=31)           cTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCg                                          >  2:78873/1‑141 (MQ=2)           cTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCg                                          <  1:78887/141‑1 (MQ=2)               aaCCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAggg                             <  2:594864/151‑1 (MQ=31)                  cgtatgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGtat                          >  2:374067/7‑150 (MQ=2)                   gtgtACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGtat                          >  2:698099/1‑150 (MQ=31)                     gtacgtacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGtat                          >  1:593911/1‑148 (MQ=31)                     atgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGa                       <  2:717562/147‑1 (MQ=2)                      tgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGAcc                     >  2:17004/3‑151 (MQ=2)                      tgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCAAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGAcc                     >  1:162371/3‑151 (MQ=2)                      tgcgtacAAGCAGTGGGAGACTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGAcc                     >  2:282960/3‑151 (MQ=2)                       gcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTATCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGAccc                    >  1:277986/2‑151 (MQ=2)                        cgtacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGAccc                    <  2:641733/151‑1 (MQ=31)                         gtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCg                   >  2:452353/1‑150 (MQ=2)                         gtacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCg                   >  1:746970/1‑151 (MQ=31)                           acAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGGGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGct                                >  1:116501/1‑134 (MQ=2)                           acAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGGGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGct                                <  2:116483/135‑2 (MQ=2)                            cAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCTGAAAc               >  1:13286/1‑151 (MQ=2)                            cAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCGAAAc               >  2:225066/1‑151 (MQ=2)                            cAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCGAAAc               >  2:165322/1‑151 (MQ=2)                            cAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTAAAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCGAAAc               >  2:241001/1‑151 (MQ=2)                              aGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAAcc              <  2:560871/151‑1 (MQ=31)                                  gTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGt          <  2:636902/151‑1 (MQ=31)                                   tGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGt          <  2:614075/150‑1 (MQ=31)                                    gggAGCATGCTTAGGCGTGTGACTGCGTTCCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTGa        <  2:455429/151‑1 (MQ=14)                                      gAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTGATc      <  1:732666/151‑1 (MQ=31)                                          aTGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAGGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTGATCTAGc  >  2:752243/1‑151 (MQ=14)                                                                                                                                                  |                                               GAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTGATCTAGC  >  NC_000913/226258‑226449

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