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	25 [0]	[3] 28	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)	25 (0.760)	24/470	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

CATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCAC  >  NC_000913/257665‑258144                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              catcTGCTGACGGTAAATTCTTAGAGATCCCATTTAACCCGGTGTATCGGACTTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTCAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACACAAGGAGATAGCAATGACGTTACCGAGAGGACACCCGAAGAGCAGATTGATCAAAAAATTTACTGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaa                                                                                                                                                                                                                                        >  2:349918‑M1/1‑243 (MQ=255)      gctgACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGTTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTTCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTAAATTCGTGAAGGTAagtgcaaagat                                                                                                                                                                                                                                    <  1:223601‑M1/249‑12 (MQ=255)           cGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgat                                                                                                                                                                                                                             <  2:140142‑M1/251‑19 (MQ=255)            ggTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGACAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAACAATTTACCGCACTAGGCCCGTATATTCATGAAGGTTagtgcaaagataatcgatt                                                                                                                                                                                                                            >  1:86666‑M1/1‑232 (MQ=255)                       agagATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgatt                                                                                                                                                                                                                 <  1:144148‑M1/251‑31 (MQ=255)                           aTCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGCGTGGACCCCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtct                                                                                                                                                                                                             >  1:307960‑M1/1‑217 (MQ=255)                                    aaCCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtat                                                                                                                                                                                                     <  1:252832‑M1/250‑43 (MQ=255)                                      cccGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtcaggctgtatgcg                                                                                                                                                                                                  <  1:162602‑M1/251‑46 (MQ=255)                                                 ggAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaa                                                                                                                                                                                        >  1:231390‑M1/1‑195 (MQ=255)                                                 ggAATTTTGACAAAGGTCTTCAGGAAATAACCGCCTGGATAGAAAAAAGCTTGTGTTAAACATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagagcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaa                                                                                                                                                                                        >  2:207346‑M1/1‑195 (MQ=255)                                                               ggTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAACATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccgg                                                                                                                                                                             >  1:310286‑M1/1‑181 (MQ=255)                                                               ggTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagagaatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccgg                                                                                                                                                                             <  2:310286‑M1/247‑67 (MQ=255)                                                                   ttCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgt                                                                                                                                                                     >  1:723‑M1/1‑177 (MQ=255)                                                                                                cGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAACACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGCGACCCGAAGAGCAGAATGATAAACAAATTTACCGCACCAGGCCCGTATATTCGTGAAGGTAcgcgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattatggggctggtggatggagctgg                                                                                                                                         >  1:168420‑M1/1‑148 (MQ=255)                                                                                                                           tGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATCGCAATGACGTTACCGAGTGGACACCCGCAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctac                                                                                                               >  2:64101‑M1/1‑121 (MQ=255)                                                                                                                                atatggtAAAACAGTTGCATCTCAACAGGAGATAGCAATGACGTTACCGAGTGGACGCCCCAAGAGCAGATTGATCAAAAAATTAACCGCACTATGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttctattgtctggctgtatgcgtcaacgtaaaaccggcaccggaagtgcgtgaattctggggctggtgg                                                                                                                                                   <  1:234867‑M1/205‑93 (MQ=255)                                                                                                                                aTTTGGAAAAACAGTTGCATCACAACAGGAGCTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattatttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtgg                                                                                                                                                   >  2:234867‑M1/1‑116 (MQ=255)                                                                                                                                  ttGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttacc                                                                                                        >  1:112690‑M1/1‑114 (MQ=255)                                                                                                                                       aaaaCAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttgg                                                                                                 >  2:222970‑M1/1‑109 (MQ=255)                                                                                                                                               tGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtaggcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatctcgtttgacctacagttaccagtttggtctgttc                                                                                          <  2:55960‑M1/250‑150 (MQ=255)                                                                                                                                                                          gTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttc                                                              <  1:87983‑M1/251‑178 (MQ=255)                                                                                                                                                                                           cGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactga                                              >  1:331589‑M1/1‑57 (MQ=255)                                                                                                                                                                                                                aaTTTACCGCACTCGGCTCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggcagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcac                        >  2:145505‑M1/1‑36 (MQ=255)                                                                                                                                                                                                                    tACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccc                    >  1:5169‑M1/1‑32 (MQ=255)                                                                                                                                                                                                                                       atTCGTGAAGGTAaatgcaacaataatcgattctttttcgattgtctggctgtatgcgtccacgtgcaaccggaaccggaagagcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagaggttgaacgactggagcacaccctgcgtgagtttcacgaga  >  2:335107‑M1/1‑13 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              CATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCAC  >  NC_000913/257665‑258144

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