Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	699,722	0	A	G	42.9%	1.4 / 12.7	14	intergenic (+115/‑129)	PP_t06/PP_23SD	tRNA‑Ala/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base G (2/4);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAA‑ATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGA  >  NC_002947/699576‑699787                                                                                                                                                    |                                                                  gggcAGCGTTTCGATCCCGCTTGGCTCCACCACTTTCGCCGGACGCAGTTACTTGTCAGAACTTAGAA‑ATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGa                                                                   <  5:127928/144‑1 (MQ=255)          ttCGCTCCCGCTTGGCTCCACCACTTTCGACGTACGCAGTTACTTGTCAGAACTTAGAA‑ATGAACATTCAGCTCTAATGTTGATTTTTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACAcc                                                          <  6:253891/147‑1 (MQ=11)             gATCCCGCTTGGCTCCACCACTTTCTCCGTACTCAGTTACTTGTCAGAACTTAGAA‑ATGAACATTCATCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAg                                                        <  7:27162/146‑1 (MQ=255)                cccGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAA‑ATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTc                                                    <  1:242496/147‑1 (MQ=37)                     ttGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAA‑ATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATAGTTCTTTAAAAATTCGGATATGTGATATATATAGACTGAACACCAGTTTCActgc                                               >  2:87528/1‑147 (MQ=255)                           ccaccaCTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAA‑ATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtgg                                         >  3:194984/1‑147 (MQ=37)                           ccaccaCTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAA‑ATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtgg                                         >  7:141717/1‑147 (MQ=37)                                        cGTACGCAGTTACTTGTCAGAACTTAGAAGATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg                              >  1:162817/1‑146 (MQ=255)                                             gCAGTTACTTGTCAGAACTTAGAA‑ATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTaaaa                        <  6:14133/146‑1 (MQ=37)                                              cAGTTACTTGTCAGAACTTAGAA‑ATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAAtt                      <  1:177011/147‑1 (MQ=37)                                               aGTTACTTGTCAGAACTTAGAA‑ATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAAttt                     >  5:106745/1‑147 (MQ=37)                                               aGTTACTTGTCAGAACTTAGAA‑ATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAAttt                     <  8:141717/147‑1 (MQ=37)                                                               ttAGAA‑ATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGGCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaa     <  6:39939/147‑1 (MQ=21)                                                                tAGAA‑ATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCg       >  1:120316/1‑144 (MQ=34)                                                                   aa‑aTGAACATTCAGCTCGAATGTTGATTTCTTGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATATATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAaga  <  7:135469/146‑1 (MQ=21)                                                                                                                                                    |                                                                  GGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAA‑ATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGA  >  NC_002947/699576‑699787

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