Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,498,182	Δ1,336 bp	IS2‑mediated	insC6–insD6	insC6, insD6

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4498182–4499511	4499517	7–1336	27 [0]	[0] 27	insC6–insD6	insC6,insD6

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4498181	0 (0.000)	27 (0.580)	26/474	0.8	100%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	4499518 =	0 (0.000)					intergenic (+18/‑1742)	insD6/yjgZ	IS2 insertion element protein InsB/uncharacterized protein YjgZ

TATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGC  >  NC_000913/4497977‑4498419                                                                                                                                                                                                             |                                                                                                                                                                                                                                               tattTTCTTTGTTAAAAAAGACAAACTGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtaca                                                                                                                                                                                                   <  2:348040‑M1/250‑46 (MQ=255)      tCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCACTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTTCTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGTGATTTGATTGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcga                                                                                                                                                                                              <  1:499831‑M1/250‑51 (MQ=255)         ttGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgt                                                                                                                                                                                          >  1:83299‑M1/1‑197 (MQ=255)                                              actatttacCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcagg                                                                                                                                                      <  2:459007‑M1/250‑91 (MQ=255)                                               ctatttacCGCTCACTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcagga                                                                                                                                                     >  1:155596‑M1/1‑159 (MQ=255)                                                    tacCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGCTTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaat                                                                                                                                                >  1:494462‑M1/1‑154 (MQ=255)                                                     acCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaata                                                                                                                                               >  1:258169‑M1/1‑153 (MQ=255)                                                     acCGCTCCCTGCAGGGGGGATTATAAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCAAAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagagcctgtatgacagtctgcaatatatgattacctgtgatcaggaaaaata                                                                                                                                               >  2:357391‑M1/1‑153 (MQ=255)                                                                 aGGGTGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCGGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtt                                                                                                                                  <  2:346‑M1/251‑111 (MQ=255)                                                                                 aTCCGCGGTTATTCATTATTTATATACATTCAATTAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAAATTACGTAAGGAGTGCTGCCGGGATTTGATGGTCCCTTGTgctgatatgaatacgggtattttatacaacatacttttcagtacaaacgactgtattaaagtctgccatatatgattacctttgatcaggaagaataaggggctggggtttcaaactttccgccgg                                                                                                                  >  2:210388‑M1/1‑125 (MQ=255)                                                                                    cgcgGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtat                                                                                                               >  1:406801‑M1/1‑122 (MQ=255)                                                                                        gTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccaaatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattc                                                                                                           >  2:21742‑M1/1‑118 (MQ=255)                                                                                                  atttatATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctccct                                                                                                 >  2:366612‑M1/1‑108 (MQ=255)                                                                                                              tCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccata                                                                                      <  2:149726‑M1/250‑155 (MQ=255)                                                                                                                    aaaaaGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgactacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggtt                                                                               >  1:195338‑M1/1‑90 (MQ=255)                                                                                                                     aaaaGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacgggtattttatacagaatacttttcagtacaagcgactgtattacatactgccatatatgattacctgtgatcaggaagaataagtgactgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggtt                                                                               >  1:488457‑M1/1‑89 (MQ=255)                                                                                                                            accacataaAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatct                                                                       <  2:494462‑M1/247‑170 (MQ=255)                                                                                                                                 aTAAAAAATCCTGTACTCACAGGGCGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcacttaaaaaccataacaggttagttatcttcatc                                                                  >  1:165564‑M1/1‑77 (MQ=255)                                                                                                                                  tAAAAAATCCTGTACTCACAGGGAGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttagtct                                                                 >  2:415714‑M1/1‑76 (MQ=255)                                                                                                                                   aaaaaaTCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctg                                                                <  1:415714‑M1/251‑177 (MQ=255)                                                                                                                                                   acaGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatgg                                                <  1:473295‑M1/251‑193 (MQ=255)                                                                                                                                                           aTGCGGGATAACTTACGTAAGGAGTGCTGCAGGTATTTGATGGTCCCTTGTgctgctatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaaatgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcgga                                            >  2:178006‑M1/1‑51 (MQ=255)                                                                                                                                                               gggATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcggtatattttc                                    >  1:277211‑M1/1‑47 (MQ=255)                                                                                                                                                                         aCGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcggtatattttcagtacaaata                          <  1:210388‑M1/251‑215 (MQ=255)                                                                                                                                                                             aaGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagagtaagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcggtatattttcagtacaaataccgg                      <  1:158464‑M1/251‑219 (MQ=255)                                                                                                                                                                                                 gATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccaaatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcggtatattttcagtacaaataccggtcatccactgaaaaatacc   <  1:21742‑M1/250‑238 (MQ=255)                                                                                                                                                                                                 gATGGTCCCTTGTgctgatatgaatacggatatgttattcagaatacttttcagtacaagcgactttattacagtctgccatatatgattacctgtgatcaggaataataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcggtatattttcagtacaaataccggtcatccactgaaaaatacct  <  2:112832‑M1/251‑239 (MQ=255)                                                                                                                                                                                                             |                                                                                                                                                                                                                                               TATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGC  >  NC_000913/4497977‑4498419

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