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)	45 (1.100)	42/490	0.1	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	3434139 =	0 (0.000)					coding (450/510 nt)	def	peptide deformylase
*	?	NC_000913	3186096 =	NA (NA)	46 (1.130)	43/490	0.1	100%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	= 3434143	0 (0.000)					coding (454/510 nt)	def	peptide deformylase

TACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAATA  >  NC_000913/3185854‑3186344                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                         aaCCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGc                                                                                                                                                                                                                                                  <  1:412435/250‑1 (MQ=255)              cAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAAGCGCAGTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGAGCAGATGTATAAAAAATTAAAATTAACCACAAGATGCGAAGTGATGAAAAATATTTATTGGTCAGTATCCGTTATTAGACTGGCCCCCTGAAtctc                                                                                                                                                                                                                                      <  2:59157/250‑1 (MQ=255)                            cAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTACTCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATc                                                                                                                                                                                                                       >  2:379436/1‑251 (MQ=255)                                gCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACtt                                                                                                                                                                                                                   <  1:226992/251‑1 (MQ=255)                                         gAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCCCCGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCAGGGATCGACCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAACATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTCTTAGACTGGCCCCCTGAATCTCCAGACAACCACTATCACTTAAATAAGTg                                                                                                                                                                                                          >  1:27327/1‑251 (MQ=255)                                                  aGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCTCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTa                                                                                                                                                                                                 <  1:65255/251‑1 (MQ=255)                                                    cAAAGATGACAACCCGGGCAACGTTTACGCTCCGGCTGACTAAAAGATTGCTTACGCGTAATCGCAATTGGCCGGGTTCGCCCGGCAACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAACAAATTAAAATTAACCACACTAAGCGAATTGATTACAAAAATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTACATACGCGATAGTCTTAAt                                                                                                                                                                                               >  2:269213/1‑251 (MQ=255)                                                     aaaGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAATTGATAGTCTTAATa                                                                                                                                                                                              <  1:375585/251‑1 (MQ=255)                                                                                 cTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAc                                                                                                                                                                  >  1:173968/1‑251 (MQ=255)                                                                                             aCAAGATTGCTTACGAGTAATAGCATTTGGCCGGGCTCGCCCGGGCACATCATTAACGGATTAATGAGAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTGATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGAAGAACAAATGATTGATg                                                                                                                                                      >  1:361593/1‑251 (MQ=37)                                                                                               aaGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTc                                                                                                                                                    >  1:436801/1‑251 (MQ=255)                                                                                                aGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCt                                                                                                                                                   <  1:313926/251‑1 (MQ=255)                                                                                                                              ggTTCGCCCGGCAACAACATTAACGGATTAATGATAAGTGGATCAGATGGATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAGACTAGTTGTTAGACTAGTCATTGGAGAACAGATGATTGAGGTCTTAGGGCCGGAGAAACGCAGACGGCGTACca                                                                                                                     <  2:173968/251‑1 (MQ=255)                                                                                                                               gTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACcac                                                                                                                    >  1:163354/1‑251 (MQ=255)                                                                                                                                               catTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGc                                                                                                     <  1:379436/250‑1 (MQ=255)                                                                                                                                                           taTGATAAGCGGCTCAGATGTATAAAAAATTAAAATTAACCACCATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTAAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGTCGGCGTACCACACAGGAAAAGATCGCAATTGTTcagc                                                                                          >  2:263453/2‑249 (MQ=255)                                                                                                                                                            aTGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTcagc                                                                                          >  2:296254/1‑248 (MQ=255)                                                                                                                                                             tGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGACATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGTCGTACCACACAGAAAAAGATCGCAATTGTTcagcag                                                                                        >  1:238821/1‑249 (MQ=255)                                                                                                                                                                     ggATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAAc                                                                              <  1:340904/251‑1 (MQ=255)                                                                                                                                                                        tCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTACAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCgg                                                                           >  2:185578/1‑251 (MQ=255)                                                                                                                                                                        tCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCgg                                                                           >  1:207059/1‑251 (MQ=255)                                                                                                                                                                         cAGATGAATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTACAGACAACCAATATCACTTCAATACGAGATAGTCTTAATCCTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGTGCCGGAGAAACGCAGGCGGCGTAGCACACAGGAAAAGATCACAATTGTTCAGCAGAGCTTTGGACCggg                                                                          >  1:397119/1‑251 (MQ=18)                                                                                                                                                                          aGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCgggg                                                                         >  1:234563/1‑251 (MQ=255)                                                                                                                                                                                                aaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATTCTAGTTGTTAGACTAGTCATTGGAGAACAGATGATTTATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGcc                                                    <  1:397336/250‑1 (MQ=255)                                                                                                                                                                                                     caATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAACAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCa                                               >  1:376415/1‑250 (MQ=255)                                                                                                                                                                                                       aTAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAAc                                             >  1:102639/1‑250 (MQ=255)                                                                                                                                                                                                                      aaaaaTATTTATTGTTCATTATCCGTTATTAGACTGGACCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTACCCGGCAACATGGTGTAGcagcc                               >  1:391743/1‑249 (MQ=255)                                                                                                                                                                                                                            atttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGtt                       >  2:264766/1‑251 (MQ=39)                                                                                                                                                                                                                                attGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTAttt                   >  1:371662/1‑251 (MQ=38)                                                                                                                                                                                                                                   gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTtctc                >  1:103456/1‑251 (MQ=35)                                                                                                                                                                                                                                      cATTATCCTTTATTATACTGGCCCCCTGAAGCACCAGTCAACCAATATCACGTAAATAAGTGATAGTCTTAAGACTAGTTTTTAGACTAGTCATTGGAGAACAGATTATTGATCTCTTCGCGCCGGAGAAACTCCGACGGCGTACCCCACAGGAAAAGATCTGAAGTGTTCAGCAGTGCTTTGAGGCTGGGAAGACGTTCCCCCGCGTTCCCCTGCAACATGGTGGAGCAGcccgccagctgcgtctctg              <  2:416851/250‑12 (MQ=2)                                                                                                                                                                                                                                                 attaGACTGGCCCCCTGAATCGCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATTCTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGGCTGAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAATa  <  2:130735/251‑1 (MQ=2)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                         TACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAATA  >  NC_000913/3185854‑3186344

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.