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–258675	258683	9–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.910)	34/466	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

AACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCA  >  NC_000913/257700‑258131                                                                                                                                                                                                                |                                                                                                                                                                                                                                 aaCCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtat                                                                                                                                                                                        <  2:56616‑M1/250‑43 (MQ=255) aaCCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGAGCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatg                                                                                                                                                                                       <  1:215747‑M1/251‑44 (MQ=255) aaCCCGGTGTAGCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCGTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgagtctttttcgattgtctggctgtat                                                                                                                                                                                        <  2:203301‑M1/250‑43 (MQ=255)       gtgtATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGCTAGAAAAACGCTTGTGTTACAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCCCAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGAAGATTGATCAAAAAATTTACCGAACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtca                                                                                                                                                                                  >  1:33277‑M1/1‑202 (MQ=255)           aTCGGAATTTTGACAAAGGTCTTCAGGAAATCTCCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgt                                                                                                                                                                              <  2:342878‑M1/250‑53 (MQ=255)                    ttGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAGGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggc                                                                                                                                                                     <  2:430263‑M1/250‑62 (MQ=255)                                  aagGAAAACTCCGACTGGATCGAAAACCGCTTGTGTGAAAAAGTTGCTAAATTTTGCCAAGTTTGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGAGCAAAAAATTTACCGCACTAGGCCCGTAGATTCGTGAAGGTAagtgcaaagataagcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaacggcaccggaagtgcgtga                                                                                                                                                      <  2:141442‑M1/250‑77 (MQ=255)                                                    aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggat                                                                                                                                    >  1:160436‑M1/1‑157 (MQ=255)                                                    aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggat                                                                                                                                    >  1:329543‑M1/1‑157 (MQ=255)                                                               cTTGTGTTAAAAATTTGCTAAATTTTGCCACTTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaacaggcaccggaagtgcgtgaattctggggctggtggatggagctgga                                                                                                                           >  2:53979‑M1/1‑146 (MQ=255)                                                                 tgtgTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcgttttcgattgtctggctgtaggcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagc                                                                                                                        <  2:1216‑M1/250‑107 (MQ=255)                                                                          aaTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcc                                                                                                              >  2:136090‑M1/1‑135 (MQ=255)                                                                                cTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaatgataatcgtttctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgtggctggtggatggagctggaagcgcaggactcccgtttt                                                                                                        >  2:13092‑M1/1‑129 (MQ=255)                                                                                cTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtttt                                                                                                        >  1:355159‑M1/1‑129 (MQ=255)                                                                                    aTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacct                                                                                                    <  2:219329‑M1/251‑127 (MQ=255)                                                                                         gCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatctattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagt                                                                                               <  1:269543‑M1/251‑132 (MQ=255)                                                                                          ccAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagtt                                                                                              >  1:381389‑M1/1‑119 (MQ=255)                                                                                             aTTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttacc                                                                                           >  2:95633‑M1/1‑116 (MQ=255)                                                                                               ttGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttacc                                                                                           >  1:233980‑M1/1‑114 (MQ=255)                                                                                                ttgtaAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATGCGTGAAGGGAagtgcaaagataatcgattcgttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagcgggaagcgcaggaatccctttttacctccagttaccag                                                                                         <  2:280510‑M1/248‑138 (MQ=255)                                                                                                      aaCAGTTGCATCACAACAGGAGAAAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACAGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgagaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtc                                                                                  >  2:182031‑M1/1‑107 (MQ=255)                                                                                                                   caacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctaca                                                                                                 >  2:315281‑M1/1‑94 (MQ=255)                                                                                                                   caacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctaca                                                                                                 <  1:315281‑M1/223‑130 (MQ=255)                                                                                                                       aGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagaaggc                                                                 >  2:252356‑M1/1‑90 (MQ=255)                                                                                                                                   tGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgacaccggcaccggaagtgcgtgaattctggggctggtggatggagctgg                                                                                                                            >  1:192898‑M1/1‑78 (MQ=255)                                                                                                                                   tGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGGATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctgg                                                                                                                            <  2:192898‑M1/180‑103 (MQ=255)                                                                                                                                       gTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggatctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttc                                                 <  1:402345‑M1/251‑178 (MQ=255)                                                                                                                                       gTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggcagagtgttc                                                 >  1:102921‑M1/1‑74 (MQ=255)                                                                                                                                         ttccgagtgggcttccgaagagcagattgaTCAAAAAATTTACCGCACGAGGCCCGTATATGCGTGAAGGTAagtgcaaagataatcgactctttttcgattgtctggctgtatgcgtcaacgtcaaaccggcaccggaagtgcttgaactctggggcgggtggatggagcttgaagcgcatgaatcccgtttcacctacagttaccagtttgtgctgttagataaagcaggcgactggaagagttttcc                                                <  2:261087‑M1/236‑179 (MQ=255)                                                                                                                                           ccGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTTCCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagttgggtctgttcgataaagcaggcgactggaagagtgttccggt                                             <  2:261443‑M1/251‑182 (MQ=255)                                                                                                                                              aGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaa                                           >  2:324120‑M1/1‑67 (MQ=255)                                                                                                                                                       ccGAAGAGCAGATTGATCAAAAACTTTACCGCACTAGGCCCGTATATTCGTGAAGGCCagtgcaacgataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatgaagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactg                                  >  2:143115‑M1/1‑58 (MQ=255)                                                                                                                                                              gCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttagattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgggaattctggggctgggggatggagctggaagcgcaggaatcccgttttacatacagttacaagtttggtctgttcgataaagcaggcgactggaagagtgttacggtaaaagacactgaggtggtt                          >  1:166876‑M1/1‑51 (MQ=255)                                                                                                                                                                      aTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgact                  <  1:104964‑M1/251‑209 (MQ=255)                                                                                                                                                                                      cACTAGGACCGTAACTTCGTGAAGGTAagtgcaaagataatcgcttctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgt  >  1:337010‑M1/1‑27 (MQ=255)                                                                                                                                                                                            gCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggcagcgcaggaatcccgtt                                                                                                          >  2:158427‑M1/1‑21 (MQ=255)                                                                                                                                                                                            gCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtt                                                                                                          <  1:158427‑M1/141‑121 (MQ=255)                                                                                                                                                                                                                |                                                                                                                                                                                                                                 AACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCA  >  NC_000913/257700‑258131

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