Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,937,768	A→T	33.3%	L429L (CTT→CTA)	PP_2570 ←	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,937,768	0	A	T	33.3%	6.0 / 4.6	12	L429L (CTT→CTA)	PP_2570	transposase
Reads supporting (aligned to +/- strand):  ref base A (6/2);  new base T (2/2);  total (8/4)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGGGATATCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAACTCCTC  >  NC_002947/2937621‑2937815                                                                                                                                                    |                                                aTCGGGATATCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGaa                                                <  6:297150/149‑1 (MQ=255)          tcactcacACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGcc                                        <  2:232846/148‑1 (MQ=255)              tcacACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATTCAACCCATTGTTTTCTAAGCACTCAAGCGGACACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGAc                                   >  5:257334/1‑149 (MQ=255)              tcacACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGAACGAACGCGCCCACAAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGAAAAGCCCGCAc                                   >  3:256589/1‑149 (MQ=17)                acacCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGcc                                        >  1:20166/1‑142 (MQ=255)                acacCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGcc                                        <  2:20166/142‑1 (MQ=255)                                  tGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTtggcttctg                <  3:240806/148‑1 (MQ=39)                                   gCCTGTGCCGCCCCCAGCGCGGGTGAACCCGCTCTCATCAAACAAAAGCCAACCCATTTTTTTCTACGCACTCAAGCGGCCACCAGTTTATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTtggcttctgg               <  6:70670/148‑1 (MQ=11)                                           cGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTACGCACTCAAGCGGCCACCAGATGATACCCCAGTTTCCGGCGTTTGCGTAGCAACGCCCGGACCTGCCGGGCTTGGCTTCTGGCTTCGCAct      >  5:256071/1‑149 (MQ=11)                                           cGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCtt                        >  5:128091/1‑131 (MQ=33)                                           cGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAAct      >  3:275480/1‑149 (MQ=32)                                            ggCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAActc     >  3:205323/1‑149 (MQ=32)                                                ccTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAActcctc  >  7:97303/1‑148 (MQ=31)                                                                                                                                                    |                                                ATCGGGATATCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAACTCCTC  >  NC_002947/2937621‑2937815

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: ‑