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

TAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCT  >  NC_000913/4499313‑4499755                                                                                                                                                                                                              |                                                                                                                                                                                                                                              tattttctttgttaaaaaagacaaactgtattaactgatgtatttactatttaccgctccctgcaggggggattatcagaatccgcggttattcattatttatatacattcaataaaaaagtaacctcataaaaaatcctgtactcacagggtgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACa                                                                                                                                                                                                   <  2:348040‑M2/45‑1 (MQ=255)      tctttgttaaaaaagacaaacggtattaactgatgtatttactatttaccgctcactgcaggggggattatcagaatccgcggttattcattatttatatacattcaataaaaaagtaacctcataaaaaatcctgttctcacagggtgatgcgggataactgacgtaaggagtgctgcagtgatttgattgtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                                                                                                                                              <  1:499831‑M2/50‑1 (MQ=255)         ttgttaaaaaagacaaacggtattaactgatgtatttactatttaccgctccctgcaggggggattatcagaatccgcggttattcattatttatatacattcaataaaaaagtaacctcataaaaaatcctgtactcacagggtgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGt                                                                                                                                                                                          >  1:83299‑M2/198‑251 (MQ=255)                                              actatttaccgctccctgcaggggggattatcagaatccgcggttattcattatttatatacattcaataaaaaagtaacctcataaaaaatcctgtactcacagggtgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAgg                                                                                                                                                      <  2:459007‑M2/90‑1 (MQ=255)                                               ctatttaccgctcactgcaggggggattatcagaatccgcggttattcattatttatatacattcaataaaaaagtaacctcataaaaaatcctgtactcacagggtgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGga                                                                                                                                                     >  1:155596‑M2/160‑250 (MQ=255)                                                    taccgctccctgcaggggggattatcagaatccgcggttattcattatttatatacattcaataaaaaagtaacctcataaaaaatcctgtactcacagggtgatgcgggataactgacgtaaggagtgctgcagggctttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAAt                                                                                                                                                >  1:494462‑M2/155‑250 (MQ=255)                                                     accgctccctgcaggggggattatcagaatccgcggttattcattatttatatacattcaataaaaaagtaacctcataaaaaatcctgtactcacagggtgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATa                                                                                                                                               >  1:258169‑M2/154‑250 (MQ=255)                                                     accgctccctgcaggggggattataagaatccgcggttattcattatttatatacattcaataaaaaagtaacctcataaaaaatcctgtactcaaagggtgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGAGCCTGTATGACAGTCTGCAATATATGATTACCTGTGATCAGGAAAAATa                                                                                                                                               >  2:357391‑M2/154‑250 (MQ=255)                                                                 agggtggattatcagaatccgcggttattcattatttatatacattcaataaaaaagtaacctcataaaaaatcctgtactcacagggtgatgcgggataactgacgtaaggagtgcggcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGtt                                                                                                                                  <  2:346‑M2/110‑1 (MQ=255)                                                                                 atccgcggttattcattatttatatacattcaattaaaaagtaacctcataaaaaatcctgtactcacagggtgatgcgggataaattacgtaaggagtgctgccgggatttgatggtcccttgtGCTGATATGAATACGGGTATTTTATACAACATACTTTTCAGTACAAACGACTGTATTAAAGTCTGCCATATATGATTACCTTTGATCAGGAAGAATAAGGGGCTGGGGTTTCAAACTTTCCGCCgg                                                                                                                  >  2:210388‑M2/126‑251 (MQ=255)                                                                                    cgcggttattcattatttatatacattcaataaaaaagtaacctcataaaaaatcctgtactcacagggtgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGtat                                                                                                               >  1:406801‑M2/123‑251 (MQ=255)                                                                                        gttattcattatttatatacattcaataaaaaagtaacctcataaaaaatcctgtactcacagggtgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCAAATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTc                                                                                                           >  2:21742‑M2/119‑251 (MQ=255)                                                                                                  atttatatacattcaataaaaaagtaacctcataaaaaatcctgtactcacagggtgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCt                                                                                                 >  2:366612‑M2/109‑251 (MQ=255)                                                                                                              tcaataaaaaagtaacctcataaaaaatcctgtactcacagggtgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATa                                                                                      <  2:149726‑M2/154‑1 (MQ=255)                                                                                                                    aaaaagtaacctcataaaaaatcctgtactcacagggtgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGACTACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGgtt                                                                               >  1:195338‑M2/91‑251 (MQ=255)                                                                                                                     aaaagtaacctcataaaaaatcctgtactcacagggtgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGGTATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACATACTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGACTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGgtt                                                                               >  1:488457‑M2/90‑250 (MQ=255)                                                                                                                            accacataaaaaatcctgtactcacagggtgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCt                                                                       <  2:494462‑M2/169‑1 (MQ=255)                                                                                                                                 ataaaaaatcctgtactcacagggcgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCACTTAAAAACCATAACAGGTTAGTTATcttcatc                                                                  >  1:165564‑M2/78‑248 (MQ=255)                                                                                                                                  taaaaaatcctgtactcacagggagatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTAGTCt                                                                 >  2:415714‑M2/77‑251 (MQ=255)                                                                                                                                   aaaaaatcctgtactcacagggtgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTg                                                                <  1:415714‑M2/176‑1 (MQ=255)                                                                                                                                                   acagggtgatgcgggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATgg                                                <  1:473295‑M2/192‑1 (MQ=255)                                                                                                                                                           atgcgggataacttacgtaaggagtgctgcaggtatttgatggtcccttgtGCTGCTATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAAATGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGcgga                                            >  2:178006‑M2/52‑246 (MQ=255)                                                                                                                                                               gggataactgacgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTc                                    >  1:277211‑M2/48‑251 (MQ=255)                                                                                                                                                                         acgtaaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATa                          <  1:210388‑M2/214‑1 (MQ=255)                                                                                                                                                                             aaggagtgctgcagggatttgatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAGTAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCgg                      <  1:158464‑M2/218‑1 (MQ=255)                                                                                                                                                                                                 gatggtcccttgtGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCAAATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATAcc   <  1:21742‑M2/237‑1 (MQ=255)                                                                                                                                                                                                 gatggtcccttgtGCTGATATGAATACGGATATGTTATTCAGAATACTTTTCAGTACAAGCGACTTTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAATAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCt  <  2:112832‑M2/238‑1 (MQ=255)                                                                                                                                                                                                              |                                                                                                                                                                                                                                              TAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCT  >  NC_000913/4499313‑4499755

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