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	1979278	776	20 [0]	[0] 20	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)	20 (0.310)	20/472	1.1	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

CATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT‑TGCGCACAGTCAAAATAACTCACC‑GTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAA  >  NC_000913/1979033‑1979503                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                   cgccgggaagccccggtaaaaaataattagcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                                                                                                                                                 >  2:118209‑M2/247‑250 (MQ=255)   ccgggaagccccggtaaaaaataattagcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGt                                                                                                                                                                                                                              >  1:98171‑M2/245‑251 (MQ=255)     gggaagccccggtaaaaaataattagcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTc                                                                                                                                                                                                                             >  1:87488‑M2/243‑250 (MQ=255)             ccggtaaaaaataattagcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                                                                                                                                                    <  2:87488‑M2/17‑1 (MQ=255)             ccggtaaaaaataattagcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                                                                                                                                                    >  1:411994‑M2/235‑251 (MQ=255)                   aaaaataattagcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaagtacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                                                                                                              <  1:227255‑M2/23‑1 (MQ=255)                                             cgataaactgaaataagcagaaccctctttttgttttaatatgtccttacaaatagaaatgggtccttacacttatctaatatttttccgaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggagggcttcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaattGATAAGTTATACTTGGGTTTTGGGTGGATTTCAATCTGTTAAATAGttt                                                                                                                                                                                    <  2:564848‑M2/49‑1 (MQ=255)                                                           aagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactagtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGGTTTGGGTGTAGTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                                                                                                                                                      <  2:159992‑M2/63‑1 (MQ=255)                                                                gaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagca                                                                                                                                                                 <  1:501405‑M2/68‑1 (MQ=255)                                                                              gtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                                                                                                                                                    >  2:90188‑M2/170‑250 (MQ=255)                                                                                                                            agatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtggtgcagatcacacaacacactcaattacttaacataaatgTATAAGGCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGATACGCTAGCAAGCAAACATGAAACAGGAATCATCGAACTGGGCTGT‑TGCGCACAGTCAAAATAACTCACC‑GTa                                                                                                    >  1:82499‑M2/124‑250 (MQ=255)                                                                                                                                         atcgacgcaactgtactcgtacctacacgcacatacacgggtggggggatgcgattttcactaatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTGATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAACAGGAAACAGGAATAATCGAAATGGGA‑GTGTGCGCACAGTCAAAATACCT‑ACCTGAAAATAATCATCTGCt                                                                                      >  2:551861‑M2/111‑251 (MQ=255)                                                                                                                                                      tactcgtcactacacgcacatacaacggaggggggctgcgattttccataatgcgtgatgcagatcacaaaaaacactcaattacttaacataactgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAACAATGAAACAGGAATAATCGAAATGGGAt                                                                                                                                   >  1:253099‑M2/98‑195 (MQ=255)                                                                                                                                                      tactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTGTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt                                                                                                                                   <  2:253099‑M2/98‑1 (MQ=255)                                                                                                                                                       actcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT‑TGCGCACAGTCAAAATAACTCACC‑GTAAATAATCATCTGCTATAAATAATCACt                                                                         <  1:90188‑M2/154‑1 (MQ=255)                                                                                                                                                                                ggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT‑TGCGCACAGTCAAAATAACTCACC‑GTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTAttt                                               <  1:118209‑M2/180‑1 (MQ=255)                                                                                                                                                                                                        atgcgtgatgcagatcactcaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT‑TGCGCACAGTCAAAATAACTCACC‑GTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCacaac                       <  2:158935‑M2/204‑1 (MQ=255)                                                                                                                                                                                                                  cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATTCTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT‑TGCGCACAGTCAAAATAACTCACC‑GTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGTCACGTACCGCACAACAAACTAATTa             <  2:98171‑M2/214‑1 (MQ=255)                                                                                                                                                                                                                       cacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT‑TGCGCACAGTCAAAATAACTCACC‑GTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATAc        >  1:158158‑M2/33‑251 (MQ=255)                                                                                                                                                                                                                             aaacactcaattacttaacataaatgTATAAGTCAGACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATTAAACAGGAATAATCGAAATGGGATGT‑TGCGCACAGTCAAAATAACTCACC‑TTAAATAATCATCTGCTATAAATTATCTCTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACaaagaa  <  1:43024‑M2/225‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                   CATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT‑TGCGCACAGTCAAAATAACTCACC‑GTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAA  >  NC_000913/1979033‑1979503

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