Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	2,069,187	0	A→C	100.0%	1.1	8	V224G (GTG→GGG)	insD1	IS2 transposase TnpB
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (1/2):  total (6/2)

CAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAG‑TCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCACCGCCGTGTTCTTCGGTTCAAGTCCCAACATC  >  NC_000913/2068901‑2069219                                                                                                                                                                                                                                                                                                |                                 cAACGCTTTGGAATGCCGTAAACAAACCCGCATTCGTGTCTGGATTTGCCCATAGTTATCCAGACATCCGTTATCACTAAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGACGAGCGATAACCCCGCGCACTATGCTGATGCCAGTCGATATAATGCTCGAACGCCTCTGCAATGTTCTTTGTTGCCGTTAACTCGTCTGGTTTGGGCATGATACTGATGTAG‑TCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCACCGCCGTGttcttc                    <  2:183986/301‑1 (MQ=11)   aCTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGCTTTGCCCCTATATTTCCAGACCTCTGCCATCACTTCACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGCTAACCCACCGCCCTATTCGGATGCCATTCGTCATACTGCTCGAACGCCCTCGCACCGTTCTTTGCTACCGTTAACCCGTCCGGTTTGGGCCTGATACTGATGTAC‑TCACGCTTTATCGTCTTCACGTCGCTCCCTGCTCTTCCGTTCCTCTCCGGCGTCCGCCCCGCCGCGttcttgtg                  >  1:241029/1‑298 (MQ=14)   aCTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATACCCCAGCGCACTCTGCGGATGCCATTCGTTATATTGATCGAACGCCGCGGCAAGGTTCGTTGCTGCCGTTAACCCGGCTGGTTTGGGCCTGATAACGATGTTGTTCAC‑CTGTATGGTTTTCACAGAGCTATTGGATTTCACGTTCCTTTCCGGTCTCCGCCccaccgcttttttacg                  >  1:696350/1‑291 (MQ=14)    cTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGCTGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCGGTCTGGTTTGGGCATGATACTGATGTAG‑TCAACATTTATCGTTTTCACGAAGCTCTCTTCTATTCCGTTGATTCCCGGCCTCCGCCCCGCCGTGTTCTTCgtt                 >  2:804424/1‑299 (MQ=14)      ccTTGGCAGGCCGTAACCAAATCCGCATTCGTGGCTGAATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACCAGCCCGCTGCCGCAGATATTCCCGTTGCGAGCGCTAACCCCGCGCAATCTGCAACTGCCATTCGTTATCATCCTCGACCGCCTCTGCCAGGTTCTTTGCTGCCGTTAACCCTTCTGGTTTGGGCCTGCTCCTGCTTTAG‑TCACGCTTTAGCGTTTTCCCGCACCGCTATTCATATCTGTTACTCTCCGGACTCCTCATCGCCTTTTTTTTCCGttt               >  1:1507780/1‑300 (MQ=14)      ccTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTTGCGAGCGATAACCCAGCGCACTATGCGGCTGCCATTCGTTATAATGCTCGAACGCCTCTGCCAGGTTCTTTGCTGCCGTTAACACGTCTGGTTTGGCCATGCTACTGATGTCT‑TCACGCTTTATCGTTTTCACTATACTCTCTTCTATTCCGTTACTCTCCAGACCCTCCCCCACCGCTTTCTTCGTtta               >  1:1051484/1‑300 (MQ=14)               gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCCAGCGATAACCCAGCCCACTATGCTGCTGCCATTCGTTCTAATGCTCGAACGCCTCTGCCAGGTTCTTTGCTGCCGTTAACCCGTCTTGTTTCGGCATGATACTGATTTAG‑TCACGCTTTAACGTTTTCACGACGCTCTCTGCTATTCCGTTACTCTCCGGACCTCGCCCCCTCgtgtcttcacgttccgttccgc       >  1:612447/1‑282 (MQ=14)                   tAAAAAAATCCGCATTAGTGGCTGGATTTGCCCCTAGATTGCCAGGCCTCTGTTATCACTTAACCCATTACAAGCCCTCTGCCGTAGATATTGCGGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTAGAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGTCATGATACTGATGTAG‑TCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCACCGCCGTGTTCTTCGGTTCAAGTCCCAACATc  <  2:1442599/301‑1 (MQ=11)                                                                                                                                                                                                                                                                                                |                                 CAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAG‑TCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCACCGCCGTGTTCTTCGGTTCAAGTCCCAACATC  >  NC_000913/2068901‑2069219

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.