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	25 [0]	[0] 24	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)	24 (0.770)	23/472	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

GACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATAT  >  NC_000913/257673‑258138                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                        gACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGAAATTTTGACAAAGGTCTTCAGGAAAGCACCGACTGGAGCGAAAAACGCTTGTTTTAAAAAGTTGCTAAATTTTGCCAATGTGGTAAAACAGTTACATCACAACAGGAGAAAACAATGACGTGAGCGAGTGGACACCCTATGAGCAGATTGATCAAAAAAAATACCTCACTAGGCCCTTATATTCGTGAaagtatgaaaaaatataatcg                                                                                                                                                                                                                         <  2:161633/251‑17 (MQ=255) gACGGTAAATTATTAGAGATCCCATGTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGATTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTGGCATCACAACAGGATATAGCAATGACGTTACCTAGTGTACACCCGAAGATCAGATTGATCAAAAAATTTGCCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatc                                                                                                                                                                                                                          <  1:139007‑M1/250‑16 (MQ=255)  acgctaaatTATTAGAGATCCCATTTAACCCTGTGTATCGGAATTTTTACAAATGTCTTCAGGAAATCACCGACTGGATCGAAAAACGATTGTGTTAAAAATTTGCTAAATGTTGCCAATTTGGTAAAACAGTTGCATCACAACAGTAGATAGCAATTACGTTACCGCTTGGACTCCCGAAGATCAGATTGATCAAAAAATTTACCGCACTAGGCCCTTATATTCGTGAAGGTAagtgcaaagatatcgga                                                                                                                                                                                                                        <  1:186570‑M1/247‑18 (MQ=255)        aaTTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGCTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGCAGAGCAGATTGCTCACAACATTTACCGCACCAGGCCCGTCTCTTCGTGACGGAAagcgcacagatactcccttctct                                                                                                                                                                                                                  >  2:68209‑M1/1‑228 (MQ=255)           tatgagagatccCATTTAACCCGGTGTATCGGAATTTTGAAAAAGGGCTTCAGGAAATCACCTATTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAGGACGTTACCGAGTGGACCTCCGAAGAGCAGATTGAGCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcttttt                                                                                                                                                                                                                <  1:4790‑M1/246‑26 (MQ=255)                    tCCCATTTAACCCGGTGTATCGGAATTTTGACAACGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaacgataatcgattctttttcgattgtctg                                                                                                                                                                                                      >  2:73499‑M1/1‑216 (MQ=255)                        aTTTAACCCGGTGTATCGGCATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctg                                                                                                                                                                                                  >  2:111143‑M1/1‑212 (MQ=255)                         tttAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgt                                                                                                                                                                                                 >  2:85115‑M1/1‑211 (MQ=255)                                                                gAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATGTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGGGAAGGTAagtgcaaagataatcgattctggttcgagtgtcgggctgtatgcgtcaacgtgaaaccgggatctgaagggcgtgaatt                                                                                                                                                          <  2:310973‑M1/251‑80 (MQ=255)                                                                 aaaTCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattc                                                                                                                                                         >  2:94601‑M1/1‑171 (MQ=255)                                                                       ccGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggc                                                                                                                                                   <  1:307437‑M1/251‑87 (MQ=255)                                                                                                     aaTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcc                                                                                                                     >  1:114905‑M1/1‑135 (MQ=255)                                                                                                        ttGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgt                                                                                                                  >  1:22684‑M1/1‑132 (MQ=255)                                                                                                            tAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtttt                                                                                                               >  2:166667‑M1/1‑128 (MQ=255)                                                                                                                      cAATTTGGGAAAAAAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagtt                                                                                                     >  2:272198‑M1/1‑118 (MQ=255)                                                                                                                           tGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagagcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggcagtgcgcgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagt                                                                                               >  2:130314‑M1/1‑113 (MQ=255)                                                                                                                               aaaaCAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttgg                                                                                           >  1:244374‑M1/1‑109 (MQ=255)                                                                                                                                            cacaACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaa                                                                              <  2:324854‑M1/251‑156 (MQ=255)                                                                                                                                              caacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagc                                                                            >  1:325219‑M1/1‑94 (MQ=255)                                                                                                                                                                  gTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAAGTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttc                                                        <  1:174280‑M1/251‑178 (MQ=255)                                                                                                                                                                      ccGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggt                                                    >  1:132209‑M1/1‑70 (MQ=255)                                                                                                                                                                                 cccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactg                                         <  2:190642‑M1/251‑193 (MQ=255)                                                                                                                                                                                                                         ccGTATATTCGTGAAGGGAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccgacaccggaagtgcgtgaattctggggctgggggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagtttc  >  2:196866‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                           gTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaag                                               <  2:19224‑M1/203‑187 (MQ=255)                                                                                                                                                                                                                           gTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaag                                               >  1:19224‑M1/1‑17 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                        GACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATAT  >  NC_000913/257673‑258138

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