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	19 [0]	[1] 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)	19 (0.310)	18/470	1.3	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

CCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCA  >  NC_000913/1978268‑1978719                                                                                                                                                                                                                                           |                                                                                                                                                                                                                          cccGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                                                                                                           <  1:186019‑M1/251‑17 (MQ=255)                                   gATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                                                                                                         >  2:489557‑M1/1‑201 (MQ=255)                                    aTAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                                                                                                        >  1:282662‑M1/1‑200 (MQ=255)                                              aTAAGCAGAACCACCTTTTTGGTTTAATCTGTCCTTACAAATAGAAATGGGCCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGCACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCACTAATGCGTGAGGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcattctgttaaaaagtttttcgctacgct                                                                                                                                                               >  2:37173‑M1/1‑190 (MQ=255)                                                   cAGAACCACCTTTTTGGTTTAATAAGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACCGTACTCGTCCCTACACGCACATTCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATTCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag                                                                                                                                                         <  1:467453‑M1/251‑67 (MQ=255)                                                                  ggTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcagactgttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                                                                                                                                          <  2:282662‑M1/251‑82 (MQ=255)                                                                      tAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACTCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGTGCTGCGATTTTCAATAATGCGTTATTCAGATCACACAAAACACTCAATTGCTTAACATAAATGtataagtcatactgttgttttgggtgtattgcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat                                                                                                                                       <  1:498055‑M1/250‑85 (MQ=255)                                                                        atatGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                                                                                                    <  1:432344‑M1/251‑88 (MQ=255)                                                                                              gggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcgacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                                                                                               <  2:552401‑M1/250‑109 (MQ=255)                                                                                                            aTCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtttcgcacagtcaaaataactca                                                                                                 <  1:199747‑M1/250‑123 (MQ=255)                                                                                                                  gATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa                                                                                          <  2:366126‑M1/251‑130 (MQ=255)                                                                                                                                                            acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaatacc                                                    >  1:6831‑M1/1‑80 (MQ=255)                                                                                                                                                                 cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTATCATAAATGtataagtcatacttttgttttgggtgtatttcaagctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgctcacagtcaaaataactcaccgtaaataatcttctgctataaataatcactttcatgcaataccagataagct                                           <  1:569818‑M1/251‑177 (MQ=255)                                                                                                                                                                                          ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat                                                                                                                         >  2:441685‑M1/1‑50 (MQ=255)                                                                                                                                                                                          ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat                                                                                                                         <  1:441685‑M1/148‑99 (MQ=255)                                                                                                                                                                                               gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcattaccagataagctatttttaaacagacacttaccgcacaaca              >  2:436921‑M1/1‑45 (MQ=255)                                                                                                                                                                                               gCGTGATGAAGATCACCCAAAACACTCAATTACTTACCATAAATGtataagtcatacttttgttttgggtgtatatcactctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaccggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctctttttaaacagacacttaccgcacaaca              >  2:362139‑M1/1‑45 (MQ=255)                                                                                                                                                                                                          ataacacaaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgataaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttgtaaacagacacttaccgcacaacaaactaattaaat  <  2:261181‑M1/248‑218 (MQ=255)                                                                                                                                                                                                           tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaat  >  2:439434‑M1/1‑33 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                          CCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCA  >  NC_000913/1978268‑1978719

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