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–258674	258683	10–776	27 [0]	[0] 27	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)	27 (0.880)	26/476	0.3	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

ATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTT‑GGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCC  >  NC_000913/258435‑258902                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                            cctcatcatctgctgacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAg                                                                                                                                                                                                                            <  1:76429‑M2/2‑1 (MQ=255)               gacggtaaattattagagatcccatttaacccggtgtatcggaagtttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATc                                                                                                                                                                                                               <  2:233289‑M2/15‑1 (MQ=255)                         tattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCttttt                                                                                                                                                                                                     <  2:264979‑M2/25‑1 (MQ=255)                                           acccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGc                                                                                                                                                                                  <  2:250240‑M2/44‑1 (MQ=255)                                            cccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCg                                                                                                                                                                                 >  2:60526‑M2/207‑251 (MQ=255)                                            cccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCg                                                                                                                                                                                 <  1:318783‑M2/45‑1 (MQ=255)                                                             ttgacaaaggtcttcaagaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaatgttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcgctaggcccgtatattcgtgaaggtaAGTTCAAAGATTAGCGATTCTGGGTCGATTGTCTGGCTGTATGCGTCAACGTGAGACCGGc                                                                                                                                                                 <  1:205507‑M2/61‑1 (MQ=255)                                                                             ggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAAt                                                                                                                                                >  1:200821‑M2/174‑251 (MQ=255)                                                                                     ccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGc                                                                                                                                        >  2:133953‑M2/166‑251 (MQ=255)                                                                                                gaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTTCGTGAATTCTGGGGCTCGtggatgga                                                                                                                             <  2:43977‑M2/97‑1 (MQ=255)                                                                                                          tgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGa                                                                                                                                                                                                                    <  2:290763‑M2/10‑1 (MQ=255)                                                                                                          tgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGa                                                                                                                                                                                                                    >  1:290763‑M2/145‑154 (MQ=255)                                                                                                              ttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAgg                                                                                                               >  1:105636‑M2/141‑251 (MQ=255)                                                                                                                                           gtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacccccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGttt                                                                                  >  1:213107‑M2/112‑251 (MQ=255)                                                                                                                                               aacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAAGCCCGTTTTACCTACAGTTACCAGTTT‑GGTc                                                                             <  2:200821‑M2/144‑1 (MQ=255)                                                                                                                                                          cacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTT‑GGTCTGTTCGATaa                                                                   <  2:71645‑M2/154‑1 (MQ=255)                                                                                                                                                                 ggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatataccgcactaggcccgtatagtcgtgaatgtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAATCTCAGGAATCCCGTTTTACCTACAGTTTCCAGTTT‑GGTCTGTTCGATAAAGCAGGCg                                                           <  1:104261‑M2/162‑1 (MQ=255)                                                                                                                                                                    gatagcaaggacggtaccgaggggaataccgaaaagcagattgatcaaaaaagttaccgcactaggcccgtatattcgggaaggtaAGTGCAAAGATAATCGATTCTTTGTCGATTGTCTGGCGGTATGCGTCAACGTGAACCCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTTTATGTAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTAGGT‑TGTCCGATAAAGCAGGCGACt                                                        <  2:57121‑M2/165‑1 (MQ=255)                                                                                                                                                                        gcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTT‑GGTCTGTTCGATAAAGCAGGCGACTGgat                                                    <  2:52166‑M2/169‑2 (MQ=255)                                                                                                                                                                          aatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgcatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTT‑GGTCTGTTCGATAAAGCAGGCGACTGGAaga                                                  <  2:93540‑M2/171‑1 (MQ=255)                                                                                                                                                                             gacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTT‑GGTCTGTTCGATAAAGCAGGCGACTGGAAGAgtg                                               <  1:267773‑M2/174‑1 (MQ=255)                                                                                                                                                                             gacgtcaccgagtggcccctcgaagagcagattgatccaaaaatttaccgctctaggcccgtatattagtgacggtgAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGCGAAACCGGGAACGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGCTACCAGTTT‑GGTCTGTTCGATAAAGCAGGCGACTGGCAGAgtg                                               >  2:157414‑M2/78‑251 (MQ=255)                                                                                                                                                                                   accgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTT‑GGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCg                                          <  2:278407‑M2/179‑1 (MQ=255)                                                                                                                                                                                                 cgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTT‑GGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGa                            >  2:69095‑M2/58‑250 (MQ=255)                                                                                                                                                                                                      agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTT‑GGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGt                      >  2:33966‑M2/53‑251 (MQ=255)                                                                                                                                                                                                                  aaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTT‑GGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTgg           >  2:90067‑M2/41‑250 (MQ=255)                                                                                                                                                                                                                          taccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTT‑GGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACAccc  >  1:280163‑M2/33‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                            ATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTT‑GGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCC  >  NC_000913/258435‑258902

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