Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,978,503	Δ776 bp		insB5–insA5	insB5, insA5

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979270	1979278	9–776	25 [0]	[0] 25	insB5–insA5	insB5,insA5

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	25 (0.450)	24/466	0.8	100%	intergenic (‑305/+16)	flhD/insB5	DNA‑binding transcriptional dual regulator FlhD/IS1 family protein InsB
	?	NC_000913	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA5/uspC	IS1 family protein InsA/universal stress protein C

GGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGG  >  NC_000913/1979036‑1979506                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                    cgggaagccccggtaaaaaataattagcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTc                                                                                                                                                                                                                              <  2:100871‑M2/8‑1 (MQ=255)      agccccggtaaaaaataattagcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgattcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                                                                                                                                                                                         <  2:544280‑M2/13‑1 (MQ=255)          ccggtaaaaaataattagcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaagagaaatgggtctttacacttatctaagatttttcctaaaacgacgcaactgtactctgcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                                                                                                                                                     <  1:594922‑M2/17‑1 (MQ=255)                          agcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                                                                                                                                                                                      <  1:230347‑M2/32‑1 (MQ=255)                          agcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                                                                                                                                                                                      >  1:431369‑M2/219‑250 (MQ=255)                               tagaatagttgcgataagctgccatacgcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttctctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                                                                                                                                                                >  2:380289‑M2/214‑251 (MQ=255)                                agaatagttgcgataagctgcaataagcataaccaccttttttgtttaatatgtccttacaaatagaaatgggtctttacacttagctaagattgttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcagataaaacactcaattacttaacataaatgTATAAGGCAGACTTTTGTTTTGGGTGTATTTCAATCTGt                                                                                                                                                                                               <  2:564768‑M2/39‑1 (MQ=255)                                         gcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                                                                                                                                                                      <  1:505931‑M2/48‑1 (MQ=255)                                          cgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTGTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                                                                                                                     <  2:80134‑M2/49‑1 (MQ=255)                                          cgataagctgcaaaatgcagaaccacctatttggttgtatatgtccgtgcaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggagggctgcgattttcaataatgcgttatgcagatcacacaaaacaggcaattacttaacagaaatgTAGAAGTCATACTTTTTTTTTGTGTGTAGTTCAAGCAGTTAAAAAGttt                                                                                                                                                                                     <  2:260973‑M2/49‑1 (MQ=255)                                                              aaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacaagcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                                                                                                                      >  1:597188‑M2/183‑198 (MQ=255)                                                              aaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatcgaagattttgccgaaatcgacgcaactgtactcgtcactacactcacatactacggaggggtgctgcgattttcaataatgcgtgatgcagatcacacaaaacacgcaacttcttaacataaatgTATAAGTCATTCtttt                                                                                                                                                                                                                      <  2:597188‑M2/16‑1 (MQ=255)                                                                         tggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTGTGGGTGTATTTCAATCTTTTAAAAAGTTGTTCGCTACTCTAGCAAGCAAAAATGAAACa                                                                                                                                                       <  2:532340‑M2/79‑1 (MQ=255)                                                                                  atgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTGTCGTTACGCTAGCAAGCAAAAATGAAACAGGAATAATc                                                                                                                                              <  1:380289‑M2/88‑1 (MQ=255)                                                                                                                                        cgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCtat                                                                                       >  1:353567‑M2/109‑251 (MQ=255)                                                                                                                                            gcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTGAAAAAGTTTTTCGCTACGCGAGCAAGCAAAAATGAAACAGTAATAATCGAAATGGGATGTTGCGCACATTCAAAATAACTCACCGTAaataa                                                                                                  <  2:299354‑M2/132‑1 (MQ=255)                                                                                                                                            gcaactgtactcgtcactacacgcacatacaacggaggggggctgagattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaataa                                                                                                  >  1:299354‑M2/105‑236 (MQ=255)                                                                                                                                                                                     gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAg                                          <  1:113117‑M2/188‑1 (MQ=255)                                                                                                                                                                                      ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAAc                                            >  1:317423‑M2/63‑248 (MQ=255)                                                                                                                                                                                                   taatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGcac                             >  1:280071‑M2/50‑250 (MQ=255)                                                                                                                                                                                                           gatgcagatcacacaaaacactaaattacttaacataaattTATTAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCATAGTCAAAATAACTCACCGTAAATAATTATCTGCTATAAATAATCTCTTTCATGCAATACCAGATAAGCTAGTGTTAAACAGACTCGTTCAGTACAACAAACTa                    <  2:353567‑M2/210‑1 (MQ=255)                                                                                                                                                                                                                agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTa                >  1:445018‑M2/37‑250 (MQ=255)                                                                                                                                                                                                                  atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATTCCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAAt             <  1:51567‑M2/217‑1 (MQ=255)                                                                                                                                                                                                                      cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACa          >  1:39949‑M2/31‑250 (MQ=255)                                                                                                                                                                                                                             cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAgg  >  1:517918‑M2/24‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                    GGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGG  >  NC_000913/1979036‑1979506

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑