Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	257,908	Δ776 bp		insB9–[crl]	insB9, insA9, [crl]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	257908	258683	776	22 [0]	[1] 23	insB9–[crl]	insB9,insA9,[crl]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 257907	0 (0.000)	22 (0.750)	22/474	0.4	100%	intergenic (+8/+16)	crl/insB9	RNA polymerase holoenzyme assembly factor Crl/IS1 family transposase B
	?	NC_000913	258684 =	0 (0.000)					pseudogene (9/331 nt)	crl	RNA polymerase holoenzyme assembly factor Crl

TATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGC  >  NC_000913/258446‑258926                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    gctgacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGa                                                                                                                                                                                                                                            <  1:311576‑M2/10‑1 (MQ=255)  ctgacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAAt                                                                                                                                                                                                                                        >  1:67676‑M2/238‑251 (MQ=255)      cggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaagttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTTCAAAGATAATCGAt                                                                                                                                                                                                                                    <  1:56053‑M2/18‑1 (MQ=255)          aaattattagatatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccggatattcgtgaaggtaAGTGCAAAGATAATCTATTCt                                                                                                                                                                                                                                 <  2:204514‑M2/21‑1 (MQ=255)                         ccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGc                                                                                                                                                                                                                >  2:68516‑M2/215‑251 (MQ=255)                               aacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggaacgaaaaacgcttgtgttaaaaatatgctaaattttgacaatttggtaaaacagttgcatcacaacaggagatagccctgaagttaccgagtggacacccgaagcgcagattgatcaaaaaatttaccgcactagacccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATg                                                                                                                                                                                                          >  1:25717‑M2/209‑251 (MQ=255)                                  ccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggctcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggcgatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGt                                                                                                                                                                                                       >  2:114805‑M2/206‑251 (MQ=255)                                      tgtatcggaattttgccaaaggtcttcaggaaatcacagaatggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagaagattgatcaaaaaatttaccgcactaggcccgaatattcgtgaagggaATTGCAAAGATACTCGCTTCTTTTTCACTTG‑CTGGGCTGTTTGCGTCAAc                                                                                                                                                                                                   >  2:107702‑M2/202‑251 (MQ=255)                                                 tttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatatgcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGTCAACGTGAAACCgg                                                                                                                                                                                         <  1:114805‑M2/60‑1 (MQ=255)                                                      caaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGTCAACGTGAAACCGGCACCgg                                                                                                                                                                                   >  2:150172‑M2/186‑251 (MQ=255)                                                                            gactggatcgaaaaacgcttgcgttcaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgcgtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGTCAACGTGAAACAGGCACCGGACGTGCGTGAATTCTGGGGCtg                                                                                                                                                             >  2:221722‑M2/164‑251 (MQ=255)                                                                                                       aaatttgctaaaatttgccaattaggtaaaacagttgcatcacaacaggagatagcaatgacgttacagaggggacacccgaagagcagattgatcaaaaaatttaccgcactaggcacgaatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCAGAATGCGTCAACGAGAAACCGGAAAAGGAAGTGCGTGAATTCTGGGGCGGGTGGATGGAGATGGAAGCGCAggaaac                                                                                                                                  >  2:126795‑M2/137‑249 (MQ=255)                                                                                                              ctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGTCACCGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGtttt                                                                                                                           >  1:172405‑M2/130‑251 (MQ=255)                                                                                                                     ttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTTGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTAc                                                                                                                     >  1:149506‑M2/123‑250 (MQ=255)                                                                                                                          aatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTAc                                                                                                               >  1:85287‑M2/118‑251 (MQ=255)                                                                                                                                        gttgcataacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGtt                                                                                                 >  1:78162‑M2/104‑251 (MQ=255)                                                                                                                                                                            agtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTaa                                                              >  1:127377‑M2/68‑250 (MQ=255)                                                                                                                                                                                 ccacccgaagagcatattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGaca                                                        <  1:126795‑M2/189‑1 (MQ=255)                                                                                                                                                                                 acacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGac                                                         >  1:211479‑M2/63‑250 (MQ=255)                                                                                                                                                                                                    atcaaaaaatttaccgcactaggcccgtatttgcgtgaagggaAGTGCAAAGATAATCGATTCTTTTTCGATGGTCT‑GGCTGTATGCGTCAACGTGAAACCGGCACCGGAATTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAAGCCCGTTTTACCTACATTTACCAGTTTGGTCTGTTCGATACAGCAGGCGACTGGAATAGTGTTCCGGTAAAAGACACTGAAGTTGTTGAACGACt                                     <  2:108497‑M2/208‑1 (MQ=255)                                                                                                                                                                                                        aaaaatttaccgcactaggcccggatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAg                                 <  1:9452‑M2/212‑1 (MQ=255)                                                                                                                                                                                                                          gcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGttt               <  1:11955‑M2/230‑1 (MQ=255)                                                                                                                                                                                                                                        gAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGc  <  1:68516/250‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    TATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT‑GGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGC  >  NC_000913/258446‑258926

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