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	10 [0]	[0] 10	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)	10 (0.200)	8/468	2.2	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

AGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAAC  >  NC_000913/1979084‑1979510                                                                                                                                                                                                    |                                                                                                                                                                                                                                        ctgcaataagcagaaccacctttttggtttaatatgaccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgccgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                                                                                                                                                                   >  1:465915‑M2/196‑250 (MQ=255)                                          ccaaatagaattgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaaggggactgcgattttcaataatgcgtgatgcatatcacacaaaacactccattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAAGCTAGCTAGCAAAAATGCAACAGGCATAATCGAAAtgggg                                                                                                                                         >  2:17032‑M2/155‑250 (MQ=255)                                              atagaaatgggtctttacacttatgtcagatttttcctaaatcgacgcaaatgtactcgtcactacacgcacatacaaaggaggggggatgcgattatcaataatgcgtgatgcagatcacacaaaacactaaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGAATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGAAAAAATGAAACAGAAATAATCGAAATGGGATGt                                                                                                                                      >  2:509004‑M2/151‑250 (MQ=255)                                              atagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGtt                                                                                                                                     >  1:16044‑M2/151‑251 (MQ=255)                                                                               ttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaataat                                                                                                     <  1:562205‑M2/133‑1 (MQ=255)                                                                                 cctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaatactgcgtgatgcagatcacacaaaacactcaattacttaacataactgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGAAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat                                                                                                  >  1:23254‑M2/116‑251 (MQ=255)                                                                                                                acgcacatacaacggaggggtgctgcgagtttcaataattcttgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTATTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCAAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATTTTTCGCACAGTCAAAATAACTCACCGTAAATAATCATCCGCTATAAATAATCCCTTTCATGCAATAcc                                                                   <  2:465915‑M2/167‑1 (MQ=255)                                                                                                                                        gcgatattcaataatgcgggaggcagatcacaaaaaacactcaattactcaacataaatgGATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCCCAGTCAAAATAACTCACCGTAAATAATCATCTGCGATAAATAATCACTTGCATGCAATACCAGATAATCTATTATCAAACAGa                                             <  2:23254‑M2/189‑1 (MQ=255)                                                                                                                                                                             cactcaattacttaacataaatgTTTAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATAGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACTACAAACTAATTAAATACAAAGAAAgg      >  2:407023‑M2/24‑251 (MQ=255)                                                                                                                                                                                  aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAggcacc  >  1:105126‑M2/19‑248 (MQ=255)                                                                                                                                                                                                    |                                                                                                                                                                                                                                        AGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAAC  >  NC_000913/1979084‑1979510

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