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	40 [0]	[0] 40	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.800)	28/260	0.5	100%	intergenic (+8/+16)	crl/insB9	RNA polymerase holoenzyme assembly factor Crl/IS1 transposase B
	?	NC_000913	258684 =	0 (0.000)					pseudogene (9/331 nt)	crl	RNA polymerase holoenzyme assembly factor Crl

GCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTT  >  NC_000913/258547‑258822                                                                                                                                          |                                                                                                                                           aaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                                                                                                          <  2:319717‑M2/3‑1 (MQ=255)             attttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATc                                                                                                                              <  1:604519‑M2/15‑1 (MQ=255)                    caatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCtt                                                                                                                       <  1:350997‑M2/22‑1 (MQ=255)                              aaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                                                                                             <  1:668424‑M2/32‑1 (MQ=255)                                       catcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTa                                                                                                    <  2:136304‑M2/41‑1 (MQ=255)                                            caacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGt                                                                                               >  2:344524‑M2/95‑140 (MQ=255)                                               caggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCaa                                                                                            >  1:140261‑M2/92‑140 (MQ=255)                                                      tagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGaaa                                                                                     >  1:691856‑M2/85‑140 (MQ=255)                                                       agcaatgacgttgccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAAc                                                                                    >  2:785332‑M2/84‑140 (MQ=255)                                                       agcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAAc                                                                                    >  2:461897‑M2/84‑140 (MQ=255)                                                       agcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAAc                                                                                    >  1:397164‑M2/84‑140 (MQ=255)                                                       agcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAAc                                                                                    >  1:805461‑M2/84‑140 (MQ=255)                                                        gcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAAcc                                                                                   <  2:632750‑M2/58‑1 (MQ=255)                                                        gcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAAcc                                                                                   <  2:140261‑M2/58‑1 (MQ=255)                                                              acgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCAcc                                                                             >  1:240881‑M2/77‑140 (MQ=255)                                                              acgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCAcc                                                                             >  1:281847‑M2/77‑140 (MQ=255)                                                                 ttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGa                                                                          >  1:74600‑M2/74‑140 (MQ=255)                                                                         tggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTg                                                                  <  2:358488‑M2/75‑1 (MQ=255)                                                                              acccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTc                                                             >  1:706276‑M2/61‑140 (MQ=255)                                                                              acccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTc                                                             >  2:348813‑M2/61‑140 (MQ=255)                                                                               cccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCt                                                            >  1:309634‑M2/60‑140 (MQ=255)                                                                                  gaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTggg                                                         >  1:563871‑M2/57‑140 (MQ=255)                                                                                            ttgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGtggat                                               <  2:563871‑M2/94‑1 (MQ=255)                                                                                                  aaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGtgga                                                <  2:783769‑M2/93‑1 (MQ=255)                                                                                                  aaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGtgga                                                >  1:783769‑M2/41‑133 (MQ=255)                                                                                                    aaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTgg                                       <  2:74600‑M2/102‑1 (MQ=255)                                                                                                       atttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAg                                    <  2:240881‑M2/105‑1 (MQ=255)                                                                                                           accgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCa                                <  2:691856‑M2/109‑1 (MQ=255)                                                                                                                 ctaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGaa                            >  2:535242‑M2/26‑138 (MQ=255)                                                                                                                 ctaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGaa                            <  1:535242‑M2/113‑1 (MQ=255)                                                                                                                    ggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCg                       >  1:331064‑M2/23‑140 (MQ=255)                                                                                                                     gcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGt                      <  2:706276‑M2/119‑1 (MQ=255)                                                                                                                          tatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTAc                 >  2:125255‑M2/17‑140 (MQ=255)                                                                                                                          tatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTAc                 >  1:343725‑M2/17‑140 (MQ=255)                                                                                                                          tatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTAc                 >  2:668359‑M2/17‑140 (MQ=255)                                                                                                                          tatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTAc                 >  2:698238‑M2/17‑140 (MQ=255)                                                                                                                            tattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCt               <  2:343725‑M2/126‑1 (MQ=255)                                                                                                                                   gaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTa        <  1:461897‑M2/133‑1 (MQ=255)                                                                                                                                    aaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTAc       >  1:742905‑M2/7‑140 (MQ=255)                                                                                                                                         aAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGtt  >  1:823238‑M2/2‑140 (MQ=255)                                                                                                                                          |                                                                                                                                           GCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTT  >  NC_000913/258547‑258822

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