Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,434,139	IS2 (+) +5 bp	coding (450‑454/510 nt)	def →	peptide deformylase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2996361 =	NA (NA)	32 (1.210)	32/490	0.1	96.9%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	3434139 =	1 (0.040)					coding (450/510 nt)	def	peptide deformylase
*	?	NC_000913	3186096 =	NA (NA)	15 (0.560)	15/490	0.6	93.6%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	= 3434143	1 (0.040)					coding (454/510 nt)	def	peptide deformylase

ATAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCC  >  NC_000913/2996139‑2996604                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    ataGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATc                                                                                                                                                                                                                          <  1:203020/250‑1 (MQ=255)                    aaaCCACATGCGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCAAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCAtt                                                                                                                                                                                                       >  2:84567/1‑250 (MQ=255)                             gtgATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                                                                                                                                              <  2:91852/250‑1 (MQ=255)                                                            gTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAAc                                                                                                                                                              <  2:105681/251‑1 (MQ=255)                                                                                    tttCCAGCGGTATCGCTGATAAATAATCCATTAGTCCATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGCCATCTGTTATCACTTAACCCATTACACGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAtt                                                                                                                                      >  2:172467/1‑251 (MQ=255)                                                                                      tCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCAATGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCg                                                                                                                                    >  2:167442/1‑251 (MQ=255)                                                                                             gTATAGCTGATAAATAATCCCTTAGTACAGCAGGAACCTTATAAAAAAACCCGCGGTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGATTAAACAGATCAATGGAAGTAAATTCCAATTGTTGATTGGATTTGCCCCTATGTTTCCAGACATCTGTTAGCACTTAACCCATTACAAGCCCGCTGCCGCATATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGAGGCCATTCGTTATAAt                                                                                                                             <  2:91690/251‑1 (MQ=25)                                                                                              tataGCTGATAAATAATCCCTCAGTACATCAAGCACCTTCTAAAAAAACCCGCGTTCTGCAAACCCTTCTCGCTAAAAGTGACAACCACTGAGTAAACCGCTCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGCGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATg                                                                                                                            >  2:16659/1‑251 (MQ=21)                                                                                                                   tAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGt                                                                                                        <  1:29832/250‑1 (MQ=255)                                                                                                                                               ccgcgtcctgctaccctctcaagaaaaaaGTGAAAATAACTGAATAAACAGATCAATCGACGTAAATTCCATTTGCTTATTGGATTTGCTCTTATATTTCAAGACATCTGTTATCACTTATACCATTAAAAGTCCGCTGCAGCAGATATCCACGTGGAGAGAGATAACCCTACCCAATATGCGGATGCCCTTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCTGTTTACCCGCCTGGTTTg                                                                           >  2:153847/26‑251 (MQ=9)                                                                                                                                                   gTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGc                                                                        >  2:240895/1‑250 (MQ=255)                                                                                                                                                   gTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGc                                                                        >  1:2241/1‑250 (MQ=255)                                                                                                                                                                    cGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGt                                                        <  2:53822/249‑1 (MQ=255)                                                                                                                                                                              gACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATc                                             <  2:2241/250‑1 (MQ=255)                                                                                                                                                                                      cTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGCTAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACg                                    >  1:231335/1‑251 (MQ=255)                                                                                                                                                                                                           attAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGATGCCGCAGATATTCCCGTGTCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATACTGCTCGCACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGtt                >  2:164924/3‑250 (MQ=255)                                                                                                                                                                                                                        tGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTcc  >  2:1891/1‑251 (MQ=31)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    ATAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCC  >  NC_000913/2996139‑2996604

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ 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.