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	38 [0]	[3] 41	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)	38 (0.800)	34/468	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

CCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCA  >  NC_000913/257659‑258143                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             ccgaatcatctgctGACGGTAAATTAGTAGAGATCCCATGTAACCCGGTGGATCGGAATTTTGACAAAGGTCGGCAGGAAAGCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCAACACAACAGGAGATAGCAATGACGTTACCGAGGGGACACCCGAAGAGCAGATTGATTAAAAAATGTACCTCACTAGGCCCGTATAGTCGTGAAGGTAag                                                                                                                                                                                                                                            <  2:206836‑M1/247‑3 (MQ=255)                    tAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGCATTTTGAAACAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCACAAAATTTACCGCACTAGGCCTGTATATTCGTGAAGGTAagtgcaacgataatcgattct                                                                                                                                                                                                                         >  1:443198‑M1/1‑230 (MQ=255)                                           aCCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgc                                                                                                                                                                                                  >  1:361352‑M1/1‑207 (MQ=255)                                               ggTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGCTTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtc                                                                                                                                                                                               >  2:154311‑M1/1‑203 (MQ=255)                                                        gAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGCTCGCAAAACGCTAGTGTTAAAAATTTGCTAAATTTTGCCAATTTCGTAAAACAGTTGCATCACAACAGGAGCTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaac                                                                                                                                                                                     >  1:208027‑M1/1‑194 (MQ=255)                                                        gAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGAATCGAAAAACGCTTGTGTTAAAAATTTGCTACATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaac                                                                                                                                                                                     >  2:205946‑M1/1‑194 (MQ=255)                                                         aaTTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtattcgtcagcgtgaaacc                                                                                                                                                                                    <  1:136302‑M1/251‑59 (MQ=255)                                                             ttGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggc                                                                                                                                                                                 <  2:349573‑M1/250‑62 (MQ=255)                                                                        cTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcg                                                                                                                                                                     >  1:179848‑M1/1‑178 (MQ=255)                                                                          tCAGGAAATCACCGACTGGATCGAAAAACGCGTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGAGAGCAATGACGGTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctgtttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgc                                                                                                                                                                      <  1:299219‑M1/248‑73 (MQ=255)                                                                           cAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtg                                                                                                                                                                   >  1:157694‑M1/1‑175 (MQ=255)                                                                                 aTCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGAGAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctg                                                                                                                                                            <  2:356292‑M1/251‑83 (MQ=255)                                                                                       gACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctg                                                                                                                                                      <  1:252891‑M1/251‑89 (MQ=255)                                                                                                                  aaagttgctaAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATGCGTGAAGGTAagtgcaaagataatctattcgttttcgattgtctgtctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggagggagctggaagcgcaggaatc                                                                                                                           <  2:41036‑M1/247‑116 (MQ=255)                                                                                                                                              aaaCAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggt                                                                                               >  2:360680‑M1/1‑108 (MQ=255)                                                                                                                                                    ttGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAACATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagt                                                                                                    >  1:301640‑M1/1‑102 (MQ=255)                                                                                                                                                    ttGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagt                                                                                                    <  2:301640‑M1/240‑139 (MQ=255)                                                                                                                                                         tCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCAGAAGAGCAGATTGATCAAAAAATTTACCGCCCCACAACAGTATATTCGTGAAGGCAagagcaaagatcatccattctttttcgattatctggctgtatgcgtcaacgtgaacccgacaccggaagtgcgtgaattctggggctgggggatggagctggaagcgcaggactcccgttttacctacagttaacagtttggtctgttcgata                                                                                     >  1:274467‑M1/1‑97 (MQ=255)                                                                                                                                                           acaacaGGAGATAGCAATGACGTTACCGCGTCGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctccagttaccagtttggtctgttcgataa                                                                                    >  2:166111‑M1/1‑95 (MQ=255)                                                                                                                                                                   agaTAGCAATGACGTTACCGAGCGGACACCCGAAGAACAGCTTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGAAagtgcaaagataatcgattctttttcgattgtctggctgtatgagtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctggtcgataaagcaggcgac                                                                          >  1:394101‑M1/1‑87 (MQ=255)                                                                                                                                                                    gattgcaatGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgaaggaatcccgttttacctacagttaccaggttggtctgttcgataaagcaggcgact                                                                         >  1:17826‑M1/5‑86 (MQ=255)                                                                                                                                                                    gaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgact                                                                         >  1:17841‑M1/1‑86 (MQ=255)                                                                                                                                                                        gCAATGACGTTACCGAGTCGACACCCGAAGAGCAGATTGATCGAAAAAGTTACCGCACTAGGCCCGTATATTCGTTAAGGTAagtgcaaagataatcgtttcgttttcgattgtctggctgtaggcgtcatcgttaaaccggcaccggaagtgggtgaattctggggctgggggatggagctggatgcgttggaatcccgttttacctacagtgaccagtttggtctgttcgataaagcaggcgactggaa                                                                     <  1:460916‑M1/251‑170 (MQ=255)                                                                                                                                                                                   aCCGAGTGGCCACCCGAAGAGCAGCTTGCTCAAAAAATTTACCGAACTAGGACCGTATATTCGTGAAGGAAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaacaggcaccggaagtgcgtgaattctggggctgggggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgaaaaagaaggcgactggaagagtgttccgg                                                          >  2:442863‑M1/1‑71 (MQ=255)                                                                                                                                                                                           gACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGCGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagac                                                  >  1:68901‑M1/1‑63 (MQ=255)                                                                                                                                                                                                 cGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagac                                                  <  2:194397‑M1/245‑189 (MQ=255)                                                                                                                                                                                                 cGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagac                                                  >  1:194397‑M1/1‑57 (MQ=255)                                                                                                                                                                                                        cAGATTGATCAAAAAATTGACCGCAC‑AGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtcgggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggatctggaagcgcaggaagcccgtgttacctacagttaccaggttggtctgttcgaggaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttga                                    <  1:214171‑M1/251‑203 (MQ=255)                                                                                                                                                                                                         aGATTGATCAAAAAATTTACCGCACTCGGCCCGCCTATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggccccggaagtgcgtgaattctggggctggtggctggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagcgttccggtaaaagacactgaagtggttg                                     >  2:317040‑M1/1‑49 (MQ=255)                                                                                                                                                                                                                  aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactgg                            >  1:355635‑M1/1‑40 (MQ=255)                                                                                                                                                                                                                      aaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcac                       <  1:361080‑M1/251‑216 (MQ=255)                                                                                                                                                                                                                         ttACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacac                     <  1:133979‑M1/250‑218 (MQ=255)                                                                                                                                                                                                                                cctAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcgttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactgtagcacaccctgcgtg             <  2:57108‑M1/250‑226 (MQ=255)                                                                                                                                                                                                                                    ggCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagtt         <  2:17841‑M1/251‑230 (MQ=255)                                                                                                                                                                                                                                    ggCCCGTATAGTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagtt         <  2:17826‑M1/251‑230 (MQ=255)                                                                                                                                                                                                                                    ggCCCATATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtcgggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggcttgtggatggagctggaagcgcaggaatcccgttttgccgacagttaccagtttggtctggtcgataaagcatgcgactggaatagtgttccggtaaaagacagtgaaggggtggaacgactggagcacacccggcgtgcatt         <  2:409616‑M1/251‑230 (MQ=255)                                                                                                                                                                                                                                     gCCCGTATATTCGTGAAGGTCagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagttt        >  2:209142‑M1/1‑21 (MQ=255)                                                                                                                                                                                                                                             atTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggatgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgacttctggggctggtggatggagctggaagcgcaggaatcccgttttacctacaggtaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaggacactgaagtggttgaacgactggagcacaccctgagtgagtttcacgag  >  2:87759‑M1/1‑13 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             CCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCA  >  NC_000913/257659‑258143

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