Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,937,768	A→T	22.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	22.3%	17.6 / 5.2	18	L429L (CTT→CTA) ‡	PP_2570	transposase
Reads supporting (aligned to +/- strand):  ref base A (10/4);  new base T (2/2);  total (12/6)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

GGATATCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAACTCCTCA  >  NC_002947/2937625‑2937816                                                                                                                                                |                                                 ggATATCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCa                                              >  1:58806/1‑148 (MQ=255)  gATATCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCaaa                                            <  5:81166/149‑1 (MQ=255)     atCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAg                                                <  5:245382/142‑1 (MQ=255)     atCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAg                                                >  6:245382/1‑142 (MQ=255)          tcacACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCCCCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGAc                                    >  8:97998/1‑149 (MQ=255)          ccacaCCTCAAGGGCGAGTTTGCCTGTGCGGCCCCCTTGGCGGGTGAACCCGCTCTCATCAAAAAAATTGCAACCCTTTTTTTTCTACCCACTCACGGGGCCACCTGTTTATACCCCTGTTTCCGGGCTTTGCTAAGCAAAGCCCGGAc                                    <  8:137468/148‑1 (MQ=14)            acacCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCCGGCTTTGCGAAGCAAAGCCCGGAc                                    >  7:137468/1‑147 (MQ=255)              acCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGc                                >  2:76482/1‑149 (MQ=255)              acCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGc                                >  2:98332/1‑149 (MQ=255)                      gTCGTGTTGGCCTTGGCCGCCCCCTCCGCGGGTCACCCCGCGCTCATAAAAAAAAATCCAACCCATTTTTTTCTACCCACACACGCGCCCACCAGTTCAAACCCCCGTTCCCGGGCTTGGCTAACAAAAGCCCGGACCTGCCGGTCTtg                        <  4:104316/149‑1 (MQ=2)                      gTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCAGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTtg                        >  5:255971/1‑149 (MQ=255)                                 cTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGccc                                        <  7:152746/122‑1 (MQ=255)                                 cTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGccc                                        >  8:152746/1‑122 (MQ=255)                                      ccGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGccc                                        <  5:232500/117‑1 (MQ=255)                                      ccGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGccc                                        >  6:232500/1‑117 (MQ=255)                                          ccccTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAActcc     >  1:214829/1‑148 (MQ=32)                                           cccTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAActcct    >  8:241627/1‑148 (MQ=32)                                            ccTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAActcgtag  >  3:182433/1‑145 (MQ=14)                                                                                                                                                |                                                 GGATATCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAACTCCTCA  >  NC_002947/2937625‑2937816

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