Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,923,718	0	G	C	20.7%	64.8 / 13.1	29	intergenic (‑38/+26)	rsmG/mnmG	16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein/5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34
Reads supporting (aligned to +/- strand):  ref base G (12/11);  new base C (2/4);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTT  >  NC_000913/3923575‑3923863                                                                                                                                                |                                                                                                                                                  tGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGc                                                                                                                                               <  6:49066/148‑1 (MQ=255) tGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGc                                                                                                                                               <  6:46713/148‑1 (MQ=255)      cATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCACCCGCTTACCa                                                                                                                                         <  5:130942/149‑1 (MQ=255)       aTTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCACCCGCTTACCAg                                                                                                                                        >  8:112334/1‑149 (MQ=255)       aTTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCACCCGCTTACCAg                                                                                                                                        <  7:112334/149‑1 (MQ=255)       aTTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGAGCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCACCCGCTTACCAg                                                                                                                                        <  5:147024/149‑1 (MQ=255)                       cACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGt                                                                                                                         <  6:115180/148‑1 (MQ=255)                                             ttCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATAcc                                                                                                  >  1:154897/1‑149 (MQ=255)                                                tGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTg                                                                                               >  2:102001/1‑149 (MQ=255)                                                          tAAGCGAAATACCTGCGTCTTCCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAAt                                                                                                                             >  1:39605/1‑109 (MQ=255)                                                                  aTACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACcag                                                                             <  1:102001/149‑1 (MQ=255)                                                                        gCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAAt                                                                       >  4:88902/1‑149 (MQ=255)                                                                                   cagTAAGGAGAGTTTGTTGAGCACGGTGATTACCGGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTc                                                                                        <  2:99071/121‑1 (MQ=255)                                                                                                  gttgAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGc                                                >  1:11586/1‑146 (MQ=255)                                                                                                       gCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGtt                                                                                             >  5:42174/1‑96 (MQ=255)                                                                                                       gCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGtt                                                                                             <  6:42174/96‑1 (MQ=255)                                                                                                             tGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGTAGATGGCCGCAGTCGTGACGACAGAAATAcgcg                                    >  8:138273/1‑147 (MQ=255)                                                                                                                 tACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGc                                                       >  4:70507/1‑124 (MQ=255)                                                                                                                 tACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGc                                                       <  3:70507/124‑1 (MQ=255)                                                                                                                      gTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGcc                          >  5:129758/1‑148 (MQ=255)                                                                                                                          ttGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGtttttt                                                                                         <  2:162858/81‑1 (MQ=255)                                                                                                                          ttGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGtttttt                                                                                         >  1:162858/1‑81 (MQ=255)                                                                                                                                  ttGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGt             >  7:142242/1‑149 (MQ=255)                                                                                                                                      cTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTtgt         <  8:142242/149‑1 (MQ=255)                                                                                                                                      cTGGTAAGCACCCGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGga                                                                    <  7:136203/90‑1 (MQ=39)                                                                                                                                      cTGGTAAGCACCCGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGga                                                                    >  8:136203/1‑90 (MQ=39)                                                                                                                                       tGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCtt                              >  1:166432/1‑127 (MQ=255)                                                                                                                                       tGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCtt                              <  2:166432/127‑1 (MQ=255)                                                                                                                                             gCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGtt  <  2:11586/149‑1 (MQ=255)                                                                                                                                               gggTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATaga                    <  7:93989/129‑1 (MQ=255)                                                                                                                                               gggTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATaga                    >  8:93989/1‑129 (MQ=255)                                                                                                                                                |                                                                                                                                                  TGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTT  >  NC_000913/3923575‑3923863

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