Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,942,761	0	A	C	37.5%	9.6 / 12.3	16	L429R (CTT→CGT) ‡	PP_2570	transposase
Reads supporting (aligned to +/- strand):  ref base A (5/5);  new base C (2/4);  total (7/9)
Fisher's exact test for biased strand distribution p-value = 6.33e-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.

GATCGGGATATCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTC  >  NC_002947/2942612‑2942798                                                                                                                                                      |                                      gATCGGGATATCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGaa                                       <  2:85908/150‑1 (MQ=255) gATCGGAATATACCTAGCCCCTGAAGGTAGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGaa                                       <  2:719319/150‑1 (MQ=17)          atCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGAACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGccc                              <  2:1100012/150‑1 (MQ=18)                    ccTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTGTCCGGGCTTTGCTTCGCAAAGCCCGGaaactggg                    <  2:1144604/150‑8 (MQ=17)                     cTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCg                             <  1:376328/140‑1 (MQ=255)                     cTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCg                             >  2:376328/1‑140 (MQ=255)                       gAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTc                 >  1:957031/1‑150 (MQ=255)                         aGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCtt               <  1:754939/150‑1 (MQ=255)                         aGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCtt               <  1:446585/150‑1 (MQ=255)                         aGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCtt               <  2:957031/150‑1 (MQ=255)                          ggTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCg                             <  1:661756/135‑1 (MQ=255)                          ggTCGAGTTTGCCTGTGACGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCg                             >  2:661756/1‑135 (MQ=255)                                gTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGcc                               <  1:821814/127‑1 (MQ=255)                                gTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGcc                               >  2:821814/1‑127 (MQ=255)                                gTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTtggcttct        >  2:396514/1‑150 (MQ=255)                                gTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTtg              >  1:646147/1‑144 (MQ=255)                                gTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTtg              <  2:646147/144‑1 (MQ=255)                                      cggtgCCTGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGTGCTTTGCGAAGCAAAGCCCTGACCTGCCGGTCTtggcttctggcttt  >  1:185245/3‑149 (MQ=255)                                                                                                                                                      |                                      GATCGGGATATCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTC  >  NC_002947/2942612‑2942798

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