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–258676	258683	8–776	37 [0]	[1] 38	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)	37 (0.680)	35/472	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

CATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCA  >  NC_000913/257662‑258140                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          catcatCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGAGTGATTAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtg                                                                                                                                                                                                                                       <  1:446976‑M1/250‑5 (MQ=255)         gctgACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagata                                                                                                                                                                                                                               <  1:373040‑M1/250‑13 (MQ=255)         gctgACGGTAAATTATTAGAGAACCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCTTATATTCGTGAAGGTAagtgcaaagat                                                                                                                                                                                                                                <  1:363612‑M1/249‑12 (MQ=255)           tgACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatc                                                                                                                                                                                                                            >  1:471903‑M1/1‑236 (MQ=255)                gTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGCATTTTGACAAAGGTATTCAGGCAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGAAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGCTTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattc                                                                                                                                                                                                                       >  1:451513‑M1/1‑231 (MQ=255)                gTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAACAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCAAAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCACAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGGAagtgcaaagctaatcgattt                                                                                                                                                                                                                       >  2:383951‑M1/1‑231 (MQ=255)                  aaaTTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattct                                                                                                                                                                                                                      <  2:413759‑M1/250‑22 (MQ=255)                            agaTCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                                                                                                                                                                                                           <  1:302837‑M1/251‑33 (MQ=255)                             gaTCCCATTTAACACGGTGTAACGGGATTTTGACAATGTTCTTCCGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTCAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACACTAGGAGAGAGCAATGACGTTACCGAGCGGACACCCGCAGAGCAGATTGATAAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAcgtgaaaagatacacggttctttttcgattgtc                                                                                                                                                                                                          >  2:332745‑M1/1‑218 (MQ=255)                                  cATTGAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGGTAAAAAGTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGGGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATGCTTGAAAGTAagtgcaaagataatcgtgtcggttgcgattgtctggct                                                                                                                                                                                                     <  2:451513‑M1/251‑39 (MQ=255)                                            ggTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAAAGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagctaatcgattctttttcgattgtctggctgtatgcgtc                                                                                                                                                                                            >  2:48883‑M1/1‑203 (MQ=255)                                                 aTCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatctattctttttcgattgtctggctgtatgcgtcaacgtg                                                                                                                                                                                      <  1:242601‑M1/251‑54 (MQ=255)                                                   cGGAATTTTGACCAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtga                                                                                                                                                                                     >  1:147552‑M1/1‑196 (MQ=255)                                                     gAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaac                                                                                                                                                                                  <  1:342506‑M1/251‑58 (MQ=255)                                                        ttttGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccgg                                                                                                                                                                               >  1:82431‑M1/1‑191 (MQ=255)                                                                                                aaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagct                                                                                                                                       >  2:126912‑M1/1‑151 (MQ=255)                                                                                                                aaTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcc                                                                                                                       >  2:298124‑M1/1‑135 (MQ=255)                                                                                                                      cTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtttt                                                                                                                 >  1:321963‑M1/1‑129 (MQ=255)                                                                                                                      cTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtttt                                                                                                                 >  2:323801‑M1/1‑129 (MQ=255)                                                                                                                                  aaTTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggagctggtggatggagctggaagcgcaggaatcccgttttacctacagttac                                                                                                     >  1:348714‑M1/1‑117 (MQ=255)                                                                                                                                      tGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagt                                                                                                 >  1:162002‑M1/1‑113 (MQ=255)                                                                                                                                       ggTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtt                                                                                                <  2:430908‑M1/251‑140 (MQ=255)                                                                                                                                           aaaCAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagggcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggt                                                                                            <  2:452525‑M1/251‑144 (MQ=255)                                                                                                                                                    cATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAACAGATTGATCAAAAAATTTACCGCACTAGTCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctgtttcgattgtctggctgtatgcgtcatcgtgaaaccggctccggaagtgcgtgaagtctggggctggtggatggagctggaagcgtaggaatcccgttttacccacagttttcgttttcgtctgttcgat                                                                                   <  2:296239‑M1/251‑153 (MQ=255)                                                                                                                                                                       aaTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctgggggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgcctgga                                                                   >  2:315034‑M1/1‑80 (MQ=255)                                                                                                                                                                          gACGTTACCGAGTGGACTCCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagt                                                              <  1:257832‑M1/250‑174 (MQ=255)                                                                                                                                                                                 ccGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccgtaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtcggttcgataaagcaggcgactggaagagtgttccggt                                                      <  2:471903‑M1/251‑182 (MQ=255)                                                                                                                                                                                              cGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTACATTAGTGAAGGTAagtgaaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaacccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagggttcaggtaaaagacactga                                          >  1:44074‑M1/1‑57 (MQ=255)                                                                                                                                                                                               gAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactga                                          >  2:363170‑M1/1‑56 (MQ=255)                                                                                                                                                                                                 agagCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtg                                      >  1:22526‑M1/1‑54 (MQ=255)                                                                                                                                                                                                        aTTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaac                               <  2:168992‑M1/251‑205 (MQ=255)                                                                                                                                                                                                              cAAAAAATTTACCGCACGAGGCCCGTATATTCGGTAAGGTAagtgcaaagataatcgatgcgttttcgagtgtctggcggtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgggggtggtggatggagctggaagggcaggaatcccgttttacctacagtttccagtttggtcgtttcgataaaggaggctactggaagagtgttccggtaaaagacactgaagtggttgaacgactgg                         <  1:280224‑M1/251‑211 (MQ=255)                                                                                                                                                                                                                aaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactgg                         >  2:193741‑M1/1‑39 (MQ=255)                                                                                                                                                                                                                 aaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagc                      <  1:107710‑M1/251‑214 (MQ=255)                                                                                                                                                                                                                        aCCGCACTAGGCCCGGGTATTCGTGAAGGTAagtgcaaagataatcgcttctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggcagtgcgtgaattctggggctggtggatgtagatggatgagcaggaatcccgttttaactacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttacggtaaaagacactgcaggggttgaacgactgaggaccaccct               >  2:268733‑M1/1‑31 (MQ=255)                                                                                                                                                                                                                                  gCCCGTATATTCGTGAAGGTAagtgcaaagataagcgattcttttgcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgggagttt     <  2:48526‑M1/251‑231 (MQ=255)                                                                                                                                                                                                                                     cGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctggtcgataaagcaggcgactggaagagtgttccggtaaaagacgctgaagtggttgaacgactggagcacaccctgcgtgagtttcac  <  2:50546‑M1/251‑234 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          CATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCA  >  NC_000913/257662‑258140

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: ‑