Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,493,846	0	A	.	62.5%	‑5.8 / 20.8	16	coding (1439/1533 nt)	PP_3985	transposase
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base . (0/10);  total (2/14)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTTCAGCGCTTCGTAGGCCCCCTTCTGACCCACGATGGGGGTAGCACTTAACTATGTATTTTAGGCTTGCAATTCCACGGCAGCAGCGCTTCGTAGTCTTCCACGGAGGCCGCCTGCGGCAGGCGCTCTAGTGCGTGGC‑GCAGCCACGCATAGGGCTCTTGGCCGTTGGCTTTGGCCGTCTCGACCAGGCTGTAAAGTTGCGCACTGGCCGTCGCCCCTCTGGGCGTGTCACTGAACAACCAGTTCTTTCTCCCGAT  >  NC_002947/4493705‑4493961                                                                                                                                               |                                                                                                                    cTTCAGCGCTTCGTAGGCCCCCTTCTGACCCACGATGGGGGTAGCACTTAACTATGTATTTTAGGCTTGCAATTCCACGGCAGCAGCGCTTCGTAGTCTTCCACGGAGGCCGCCTGCGGCAGGCGCTCTAGTGCGTGGCTGC‑GCCACGCa                                                                                                             <  1:3125660/150‑1 (MQ=0) cTTCAGCGCTTCGTAGGCCCCCTTCTGACCCACGATGGGGGTAGCACTTAACTATGTATTTTAGGCTTGCAATTCCACGGCAGCAGCGCTTCGTAGTCTTCCACGGAGGCCGCCTGCGGCAGGCGCTCTAGTGCGTGGCTGC‑GCCACGCa                                                                                                             <  1:6415168/150‑1 (MQ=0)  ttCAGCGCTTCGTAGGCCCCCTTCTGACCCACGATGGGGGTAGCACTTAACTATGTATTTTAGGCTTGCCAACCCACGGCAGCAGCGCTTCGGAGTCTTCCACGGAGGCCGCCTGCGGCAGGCGCTCTAGTGCGTGGCTGC‑GCCACGCac                                                                                                            <  1:880343/150‑2 (MQ=1)  ttCAGCGCTTCGTAGGCCCCCTTCTGACCCACGATGGGGGTAGCACTTAACTATGTATTTTAGGCTTGCAATTCCACGGCAGCAGCGCTTCGTAGTCTTCCACGGAGGCCGCCTGCGGCAGGCGCTCTAGTGCGTGGCTGC‑GCCACGCac                                                                                                            <  1:5039574/150‑2 (MQ=0)  ttCAGCGCTTCGTAGGCCCCCTTCTGACCCACGATGGGGGTAGCACTTAACTATGTATTTTAGGCTTGCAATTCCACGGCAGCAGCGCTTCGTAGTCTTCCACGGAGGCCGCCTGCGGCAGGCGCTCTAATGCGTGGCTGC‑GCCACGCac                                                                                                            <  2:4230361/150‑2 (MQ=1)    cAGCGCTTCGTAGGCCCCCTTCTGACCCACGATGGGGGTAGCACTTAACTATGTATTTTAGGCTTGCAATTCCACGGCAGCAGCGCTTCGTAGTCTTCCACGGAGGCCGCCTGCGGCAGGCGCTCTAGTGCGTGGCTGC‑GCCACGCacta                                                                                                          <  2:5348230/150‑4 (MQ=0)    cAGCGCTTCGTAGGCCCCCTTCTGACCCACGATGGGGGTAGCACTTAACTATGTATTTTAGGCTTGCAATTCCACGGCAGCAGCGCTTCGTAGTCTTCCACGGAGGCCGCCTGCGGCAGGCGCTCTAATGCGTGGCTGC‑GCCACGCacta                                                                                                          <  1:6463873/150‑4 (MQ=1)    cAGCGATTCGTAGGCCCCCTTCTGACCCACGATGGGGGTAGCACTTAACTATGTATTTTAGGCTTGCAATTCCACGGCAGCAGCGCTTCGTAGTCTTCCACGGAGGCCGCCTGCGGCAGGCGCTCTAATGCGTGGCTGC‑GCCACGCacta                                                                                                          <  1:1916604/150‑4 (MQ=1)      gcgcTTCGTAGGCCCCCTTCTGACCCACGATGGGGGTAGCACTTAACTATGTATTTTAGGCTTGCACTTCCACGGCAGTAGCGCTTCGATTGCATGGTCTGAGGCCGCCTGCGGCAGGCGCTCTAGTGCGTGGCTGC‑GCCACGCactaga                                                                                                        <  1:6245566/150‑6 (MQ=1)      gcgcTTCGTAGGCCCCCTTCTGACCCACGATGGGGGTAGCACTTAACTATGTATTTTAGGCTTGCAATTCCACGGCAGCAGCGCTTCGTAGTCTTCCACGGAGGCCGCCTGCGGCAGGCGCTCTAGTGCGTGGCTGC‑GCCACGCactaga                                                                                                        <  1:5231107/150‑6 (MQ=1)                               cgcgttGGGGGTAGCACTTAATTATGTATTTTGGGCTTGCAATTCCGCGGCAGCAGCGCTTCGTAGTCTTCCACGGAGGCCGCCTGCGGCAGGCGCTCTAGTGCGTGGC‑GCAGCCACGCATAGGGCTCTTGGCCGTTGGCTTTGGCCGtc                                                                               <  2:2764944/145‑1 (MQ=0)                                       ggttAGCACTTAACTATGTATTTTAGGCTTGCAATTCCACGGCAGCAGCGCTTCGTAGTCTTCCACGGAGGCCGCCTGCGGCAGGCGCTCTAGTGCGTGGC‑GCAGCCACGCATAGGGCTCTTGGCCGTTGGCTTTGGCCGTCTCGACCAg                                                                       >  1:3778974/4‑150 (MQ=1)                                                                                  aggagcGCTTCGTAGTTTCGCACGGAGGCCGCCTGCGGCAGGCGCTCTAGTGCGTGGC‑GCAGCCACGCATAGGGCTCTTGGCCGTTGGCTTTGGCCGTCTCGCCCAGGCTGTAAAGTTGCGCACTGGCCGTCGCCCCTCTGGGCGTGTCa                            <  2:1631123/147‑1 (MQ=1)                                                                                   tcatcgcgTCGTACTCTTCCCCGGCCGCCGCGTGCCGCCTGCTCTCTGCCGGTTGGC‑GCAGCCACTCATAGGGCGCTTGGCCGTTGGCTTTGGCCGTCTCGACCAGGCTGTAAAGTTGCGCACAGGCCGTCGCCCCTCTGGGCGTGTCAc                           <  2:3605944/146‑1 (MQ=0)                                                                                                          tagngcgccTGCCGCAGGCGCTCTAGTGCGTGGC‑GCAGCCACGCATAGGGCTCTTGGCCGTTGGCTTTGGCCGTCTCGACCAGGCTGTAAAGTTGCGCACTGGCCGTCGCCCCTCTGGGCGTGTCACTGAACAACCAGTTCTTTCTCCCg    >  2:4470225/6‑150 (MQ=0)                                                                                                            gggcgccTGCGGCAGGCGCTCTGGCGCGTGGC‑GCAGCCACGGATAGGGCTCTTGGCCGTTGGCTTTGGCCGTCTCGAGCAGGCTGTAAAGGTGCGCACTGGCCGTCGCCCCTCTGGGCGTGTCACTGAACAACCAGTTCTTTCTCCCGAt  <  1:1799924/147‑1 (MQ=1)                                                                                                                                               |                                                                                                                    CTTCAGCGCTTCGTAGGCCCCCTTCTGACCCACGATGGGGGTAGCACTTAACTATGTATTTTAGGCTTGCAATTCCACGGCAGCAGCGCTTCGTAGTCTTCCACGGAGGCCGCCTGCGGCAGGCGCTCTAGTGCGTGGC‑GCAGCCACGCATAGGGCTCTTGGCCGTTGGCTTTGGCCGTCTCGACCAGGCTGTAAAGTTGCGCACTGGCCGTCGCCCCTCTGGGCGTGTCACTGAACAACCAGTTCTTTCTCCCGAT  >  NC_002947/4493705‑4493961

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑