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–258672	258683	12–776	29 [0]	[2] 31	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)	29 (0.770)	27/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

CTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCAC  >  NC_000913/257660‑258141                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           cTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGCTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                                                                                                                                                                                                         >  1:47206‑M1/1‑248 (MQ=255)   catcatCTGCTGAGGGTAAATTATTAGAGATCCCATTTAACCCGGTGTCTAGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgc                                                                                                                                                                                                                                       >  2:370799‑M1/1‑246 (MQ=255)           gctgACGGTAAAGTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataa                                                                                                                                                                                                                               <  2:137111‑M1/251‑14 (MQ=255)                     aaTTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGAATGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcttt                                                                                                                                                                                                                     >  1:367152‑M1/1‑228 (MQ=255)                         attaGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAACTCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGGTAAATTTTGCCAATTTGGTAAAACAGATGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCACAGTGATCAAAAGATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagatattcgcttctttttc                                                                                                                                                                                                                  <  2:281074‑M1/250‑27 (MQ=255)                         attaGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttc                                                                                                                                                                                                                  <  1:147874‑M1/250‑27 (MQ=255)                           taGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcga                                                                                                                                                                                                                <  1:196501‑M1/250‑29 (MQ=255)                                aTCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCTCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtct                                                                                                                                                                                                          <  2:148701‑M1/251‑35 (MQ=255)                                 tCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctg                                                                                                                                                                                                         >  1:96580‑M1/1‑216 (MQ=255)                                  cccATTTACCCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTTGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACTGAGTGGACACCCGAAGAGCAGATTAATCAAAAAATGTACCGCACTAGGCCCGGATATTCGTGAAGGTAagtgcaaagataatcgatacgttttcgattgtctgg                                                                                                                                                                                                        <  2:201443‑M1/251‑37 (MQ=255)                                  cccATTTAACCAGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATAACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGCTAGCAATGACGTTACCGAGTGGACACCCGACGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctgg                                                                                                                                                                                                        >  2:228894‑M1/1‑215 (MQ=255)                                                 gtATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaac                                                                                                                                                                                          <  1:393254‑M1/250‑51 (MQ=255)                                                         aTTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccg                                                                                                                                                                                 <  1:167909‑M1/251‑60 (MQ=255)                                                                      tCTTCGGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATAACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAACAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaacccggcaccggaagtgc                                                                                                                                                                    >  2:151480‑M1/1‑179 (MQ=255)                                                                              aaaTCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattc                                                                                                                                                            >  1:98066‑M1/1‑171 (MQ=255)                                                                                         tGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGGCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAACGACGTTACCGAGTGGCCACCCGAAGAGAAGATTGACCAAAAAATTTACCGCACTAGCCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgagaattctggggctggtg                                                                                                                                                 >  2:160832‑M1/1‑160 (MQ=255)                                                                                                                                   cAATGTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCTCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttctattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagtta                                                                                                       <  1:238091‑M1/251‑134 (MQ=255)                                                                                                                                         ggTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagt                                                                                                  <  1:94711‑M1/250‑139 (MQ=255)                                                                                                                                            aaaaCAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAAGTTACCGCACGAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttg                                                                                               <  1:228894‑M1/250‑142 (MQ=255)                                                                                                                                                 aGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAAGTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgt                                                                                         <  2:307803‑M1/251‑148 (MQ=255)                                                                                                                                                                   gaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgg                                                                                                                                                         <  2:12971‑M1/169‑84 (MQ=255)                                                                                                                                                                   gaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgg                                                                                                                                                         >  1:12971‑M1/1‑86 (MQ=255)                                                                                                                                                                            gACGTTACCGAGTGGACACCCGAAGAGCAGATTGATAAAAAAATTTACCGCACTAGGACAGTATATTCGTGAAGGTAagtgcaatgataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatccagttttacctacagctaccagtttggtctgttcgataaagcaggcgactggaagagtg                                                              >  2:99631‑M1/1‑77 (MQ=255)                                                                                                                                                                              cGTTACCGAGTGGACACCCGAAGAGCAGCTTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagcgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcagccgactggaagagtgat                                                            >  2:168028‑M1/1‑75 (MQ=255)                                                                                                                                                                                 tACCGAGTGAACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtccacgtgaaaccggctccggaagtgcgtgaattctggggctggcggatggagctggacgcgcatgaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggagattggaagagtgttccg                                                         >  2:219251‑M1/1‑72 (MQ=255)                                                                                                                                                                                                  aaGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagt                                        >  1:337826‑M1/1‑55 (MQ=255)                                                                                                                                                                                                               tCAAAACATTCACCGCACTAGGCCCGTATATTCGTGACGGTAagtgcaaagataatcgattctttttcgattgtctggatgtatgcgtaaacgtgaaaccggcaccggaagtgcgtgaattgtggggctggtggatggagctggaagcgcaggaataccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagactgttccggtaacagccactcaagtggtggaacgtgtg                           >  2:259856‑M1/1‑42 (MQ=255)                                                                                                                                                                                                                                        gtttattcgTGAAGGAAagagcaaagataatcgattattttttgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtggattctggggttggtggatggagctggaagcgcaggaataccgttttacctacagttaccagtttggtctgttcgataacgcaggcgactggaagagtgttccggtaaaagatactgaagtggttgaacgactggagcaaaacctgcgtgagtttcacg  >  2:128694‑M1/4‑17 (MQ=255)                                                                                                                                                                                                                                        gTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtacaagacactgaagtggttgaacgactggagcacaccctgcgtgagtttcacg  >  1:251517‑M1/1‑17 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           CTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCAC  >  NC_000913/257660‑258141

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