Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,704,350	Δ390 bp		blr	blr

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1704350	1704739	390	28 [0]	[0] 27	blr	blr

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1704349	0 (0.000)	27 (0.780)	27/486	0.4	100%	intergenic (‑41/‑202)	ydgJ/blr	putative oxidoreductase YdgJ/beta‑lactam resistance protein
	?	NC_000913	1704740 =	0 (0.000)					intergenic (+64/‑209)	blr/cnu	beta‑lactam resistance protein/H‑NS‑ and StpA‑binding protein

ACCCCCTGAAGTTAATGGATCAATCAAGAGAGATGTGGGCTGTAATGAATCGTCTTATTGAATTAACAGGTTGGATCGTTCTTGTCGTTTCAGTCATTCTTCTTGGCGTGGCGAGTCACATTGACAACTATCAGCCACCTGAACAGAGTGCTTCGGTACAACACAAGTAAGCTCTGCACTTGTGGAGCGACATGCTGCCCGTCCGGGTGCATGTTTTCACTTGTCGGATATTAAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCATTCTTATTATTGTAATAACACAAGCGTTCACTTTTGGTGAAATTTTATCATCACTAATAATACCTACCCCTTTGGTATTAATTATTTTATGGACCTTTTATGACTGTTCAGGACTACTTATTAAAATTTCGCAA  >  NC_000913/1704507‑1704983                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    gttgggatcgttaaacagatgcttcggctcagagacaaccgtaaccgttggccagtcggcttttacttttgtttcatcactgctggagattaccgccagttcctgcccgggcgtgcccgcaatcaggggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATc                                                                                                                                                                                                                                    <  1:176887‑M2/18‑1 (MQ=255)  ttgggatcgttaaacagacgcttcggctcagagacaaccgtaaccgttggccagtcggcttttacttttgtttcatcactgctggagattaccgccagttcctgcccgggcgtgcccgcaatcaggggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCt                                                                                                                                                                                                                                   >  2:54182‑M2/233‑251 (MQ=255)      gatcgttaaacagatgcttcggctcagagacaaccgtaaccgttggccagtcggcttttacttttgtttcatcactgctggagattaccgccagttcctgcccgggcgtgcccgcaatcaggggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGtt                                                                                                                                                                                                                               >  1:129479‑M2/229‑251 (MQ=255)         cgttaaacagatgcttcggctcagagacaaccgtaaccgttggccagtcggtttttacttttgtttcatcactgctggagattaccgccagttcctgcccgggcgtgcccgcaatcaggggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTAGCTTGTTCg                                                                                                                                                                                                                             <  1:289477‑M2/25‑1 (MQ=255)                 agatgcttcggctcagagacaaccgtaaccgttggccagtcggcttttacttttgtttcatcactgctggagattaccgccagttcctgcccgggcgtgcccgcaatcaggggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGttgttg                                                                                                                                                                                                                    >  1:671‑M2/218‑251 (MQ=255)                   atgctgcgtctcagagacaaccgtaaccgttggccagtcggcttttacttttgtttcatcacggctgtagactaccgccagttcctgcccgggcgcgcccgcaatcaggggcgcatggaaggttttgctcgcataaccatacccaatcaacctaacacggatgttgtcgctcatgtcattgcctctcactgaagtatgatggctatttgaccctaAACCAGGAATTTAGTATCTTGTTCGATGTTGTTGGt                                                                                                                                                                                                                  <  1:274787‑M2/36‑1 (MQ=255)                         gggctcagagacacccgtaaccgttggccagtcggcttttacttttgtttcatcactgctggagattaccgccagttcctgcccgggcgtgcccgcaatcaggggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGt                                                                                                                                                                                                            <  1:33871‑M2/42‑1 (MQ=255)                               agagacaaccgtaaccgttggacagtcggcttttccttttgtttcatcactgctggagattaccgccagttcctgcccgggcgtacccgcaatcaggggcgcatggcaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGAGTGTCAggg                                                                                                                                                                                                       >  2:161384‑M2/204‑250 (MQ=255)                                               gttggccagtcggcttttacttttgtttcatcactgctggagattaccgccagttcctgcccgggcgtgcccgcaatcaggggcgcatggaaggttttgctcgcataaccatacccaatcaacctaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTAc                                                                                                                                                                                      >  2:248882‑M2/188‑251 (MQ=255)                                                         cggcttttgctttgttttcatcactgctggtgatttccgccagttcctgcccgggcgtgcccgcaagcaggggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGGAAAGCGCGTTGCTCTTTTCaga                                                                                                                                                                            <  2:43536‑M2/74‑1 (MQ=255)                                                           gcttttacttttgtttcctcatcgtgggagatttccgccagttcctgcccgggcgtgcccgcaatcagtggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctagttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGTTTGTCAGGGATAGTAAAGCGGGTTACTCTTTTCAGAGc                                                                                                                                                                          <  2:239976‑M2/76‑1 (MQ=255)                                                                         ttcatcactgctggagattaccgccagttcctgcccgggcgcgcccgcaatcaggggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATccc                                                                                                                                                             <  2:117226‑M2/89‑1 (MQ=255)                                                                                      gagattaccgccagttcctgcccgggcgtgcccgcaatcaggggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACGCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCAt                                                                                                                                               <  1:54475‑M2/103‑1 (MQ=255)                                                                                             ccgccagttcctgcccgggcgtgcccgcaatcaggggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcagtgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCa                                                                                                                                        <  2:297545‑M2/110‑1 (MQ=255)                                                                                                  agttcctgcccgggcgtgcccgcaatcaggggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCATTCtt                                                                                                                                   <  1:325788‑M2/115‑1 (MQ=255)                                                                                                          cccgggcgtgcccgcaatcaggggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCATTCTTATTATTGt                                                                                                                           <  1:20558‑M2/123‑1 (MQ=255)                                                                                                                             aggggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCAt                                                                                                                                       <  1:353530‑M2/111‑1 (MQ=255)                                                                                                                             aggggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCAt                                                                                                                                       >  2:353530‑M2/110‑220 (MQ=255)                                                                                                                               gggcgcatggaaggttttgctcgcataaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCATTCTTATTATTGTAATAACACAAGCGTTCACttt                                                                                                      <  1:4816‑M2/144‑1 (MQ=255)                                                                                                                                                         aaccatacccaatcaacccaacacggatgttgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGCATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCTGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCCTTCTTATTATTGTAATAACACAATCGTTCACTTTTGGTGAAATTTTATCATCACtaata                                                                             >  2:242058‑M2/82‑250 (MQ=255)                                                                                                                                                                        acccaacacggatgttgtcgctcatgacattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCATTCTTATTATTGTAATAACACAAGCGTTCACTTTTGGTGAAATTTTATCATCACTAATAATACCTACCCCTTTgg                                                             >  2:311289‑M2/67‑251 (MQ=255)                                                                                                                                                                                       tgtcgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCATTCTTATTATTGTAATAACACAAGCGTTCACTTTTGGTGAAATTTTATCATCACTAATAATACCTACCCCTTTGGTATTAATTATTTTAt                                              >  2:302704‑M2/52‑251 (MQ=255)                                                                                                                                                                                          cgctcatgtcattgcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCATAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCATTCTTATTATTGTAATAACACAAGCGTTCACTTTTGGTGAAATTTTATCATCACTAATAATACTTACCCCGTTGGTATTAATGATTTTATgg                                            <  1:113890‑M2/202‑1 (MQ=255)                                                                                                                                                                                                       gcctctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCATTCTTATTATTGTAATAACACAAGCGTTCACTTTTGGTGAAATTTTATCATCACTAATAATACCTACCCCTTTGGTATTAATTATTTTATGGACCTTTTATGACTg                              >  1:39499‑M2/36‑251 (MQ=255)                                                                                                                                                                                                          tctcattgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCATTCTTATTATTGTAATAACACAAGCGTTCACTTTTGGTGAAATTTTATCATCACTAATAATACCTACCCCTTTGGTATTAATTATTTTATGGACCTTTTATGACTGTTc                           >  2:308639‑M2/33‑251 (MQ=255)                                                                                                                                                                                                                tgaagtatgatggctatttgacactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCATTCTTATTATTGTAATAACACAAGCGTTCACTTTTGGTGAAATTTTATCATCACTAATAATACCTACCCCTTTGGTATTAATTATTTTATGGACCTTTTATGACTGTTCAGGact                     >  1:125789‑M2/27‑251 (MQ=255)                                                                                                                                                                                                                                    acactaAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCATTCTTATTATTGTAATAACACAAGCGTTCACTTTTGGTGAAATTTTATCATCACTAATAATACCTACCCCTTTGGTATTAATTATTTTATGGACCTTTTATGACTGTTCAGGACTACTTATTAAAATTTCGCaa  >  1:179832‑M2/7‑250 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    ACCCCCTGAAGTTAATGGATCAATCAAGAGAGATGTGGGCTGTAATGAATCGTCTTATTGAATTAACAGGTTGGATCGTTCTTGTCGTTTCAGTCATTCTTCTTGGCGTGGCGAGTCACATTGACAACTATCAGCCACCTGAACAGAGTGCTTCGGTACAACACAAGTAAGCTCTGCACTTGTGGAGCGACATGCTGCCCGTCCGGGTGCATGTTTTCACTTGTCGGATATTAAACCAGGAATTTATTATCTTGTTCGATGTTGTTGGTGATTGTCAGGGATAGTAAAGCGCGTTACTCTTTTCAGAGCAGTTTTTCATCCCTGAATATTTGCCATCAGCGCATTCTTATTATTGTAATAACACAAGCGTTCACTTTTGGTGAAATTTTATCATCACTAATAATACCTACCCCTTTGGTATTAATTATTTTATGGACCTTTTATGACTGTTCAGGACTACTTATTAAAATTTCGCAA  >  NC_000913/1704507‑1704983

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