Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,884,742	0	G	A	66.7%	9.5 / 11.2	18	D84N (GAT→AAT)	rnpA	RNase P protein component
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base A (8/4);  total (12/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACGGAATTGAGGCATTGCGC  >  NC_000913/3884575‑3884953                                                                                                                                                                        |                                                                                                                                                                                                                    ntGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGAATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGcgccgc                                                                                                                                  >  1:64601/2‑251 (MQ=255)   ggCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAACGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCgccgcc                                                                                                                                 >  1:283048/1‑250 (MQ=255)    gCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCGCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCACTCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCCGCGCCAACATGAACTCCCGGCTATGAATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCt                                                                                                                                                                                >  1:275218/1‑202 (MQ=255)           ncAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATaa                                                                                                                                                                          >  1:197225/2‑201 (MQ=255)                                   gAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCCAGAAAAACGATCGACGAGCCCATGAACGAAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCAAACATGAACTCCCGGCTATGAATTTCgcggtggt                                                                                                                                                                                                        >  1:168542/1‑147 (MQ=255)                                         gCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGAATTTCGTGGTGGTGGCGAAAAAAGGGGTTGC‑AACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCCTTCgg                                                                                         >  1:151927/1‑251 (MQ=255)                                                                  nnaGNNNNNAAGNAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATaa                                                                                                                                                                          <  2:197225/144‑1 (MQ=255)                                                                         cnAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGAATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGt                                                                                       >  1:253618/1‑222 (MQ=255)                                                                                    gTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAACTTCCGTCTGCGCCCACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGCAAAAAGGGGTTGCCGCCCTCGATa                                                                                                                                                                           >  1:301448/1‑127 (MQ=255)                                                                                                        cAATCGGATTAAACTTCTTACGCGTGAAAGCTTCCGTCTTCGCGAACATCAACCCCCTGCTATGAAGTTTGTGGTTGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCg                                                                                                                                                       <  1:315093/127‑1 (MQ=255)                                                                                                                                 naaaGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGAATTTCGTGGTGGTGGCGAAAAAAg                                                                                                                                                                                             >  1:200741/2‑64 (MQ=255)                                                                                                                                 gAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGAATTTCGTGGTGGTGGCGAAAAAAg                                                                                                                                                                                             <  2:200741/64‑1 (MQ=255)                                                                                                                                 gAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGAATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACGGAATTGAGGCATTgcgc  <  1:258647/251‑1 (MQ=255)                                                                                                                                      cTTCCGTCTGCGCCAACATGAACTCCCGGCTATGAATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGc                                                                                                                                                        <  1:188304/96‑1 (MQ=255)                                                                                                                                      cTTCCGTCTGCGCCAACATGAACTCCCGGCTATGAATTTAGTGGTGGAGGAGAAAAAAGGGGTTGCAGACCTAGATAACAGTGTTCTCTCGGCAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTcc                                                                                                          >  1:232741/1‑142 (MQ=255)                                                                                                                                                    aaCATGAACTCCCGGCTATGAATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCACGGGTCCTGATAACCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACg                   >  1:269402/1‑215 (MQ=255)                                                                                                                                                        tGAACTCCCGGCTGTGGATTTCGTGGGGGTGGCGAAAAAAGGGGTTGCCGACCCCGATAACCGTGCTCTCTCGGAAGCGGTGGAAAAATTATGGCGCCGCCACTGTCGCCGGGTTCGCGGGTCCTGATAGCCCTCATTTGGGTCTATCa                                                                                 <  1:130995/149‑1 (MQ=255)                                                                                                                                                        tGAACNCCCGGNTNNGGATTTCGTGGTGGTGGCGNAAAAAGGGGTTGCCGNCCTCNATAACNgancn                                                                                                                                                                   >  2:130995/1‑63 (MQ=255)                                                                                                                                                                        |                                                                                                                                                                                                                    CTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACGGAATTGAGGCATTGCGC  >  NC_000913/3884575‑3884953

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