Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913_3_bme_pgi	1,299,499	Δ1,199 bp	intergenic (+254/‑485)	ychE → / → oppA	UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_bme_pgi	1299499–1300697	1300697	1–1199	38 [0]	[0] 38	ychE/oppA	UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 1299498	0 (0.000)	35 (0.800)	33/278	0.4	100%	intergenic (+253/‑1684)	ychE/oppA	UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein
	?	NC_000913_3_bme_pgi	1300698 =	0 (0.000)					intergenic (+1453/‑484)	ychE/oppA	UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

TAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATA  >  NC_000913_3_bme_pgi/1300558‑1300845                                                                                                                                             |                                                                                                                                                    ttacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATc                                                                                                                                             <  1:466280‑M2/9‑1 (MQ=255)   acagttagtgataaggagttccgtaactgcttgtcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcttgtggtagaggaaagtgctaaATAATAATCaa                                                                                                                                           <  2:170617‑M2/11‑1 (MQ=255)     agttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAAtt                                                                                                                                         <  2:287430‑M2/13‑1 (MQ=255)           ttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTaaa                                                                                                                                   <  1:114877‑M2/19‑1 (MQ=255)                aggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAattatt                                                                                                                              <  2:157899‑M2/24‑1 (MQ=255)                aggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcatgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAattatt                                                                                                                              <  2:359570‑M2/24‑1 (MQ=255)                          taactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCa                                                                                                                    <  1:22505‑M2/34‑1 (MQ=255)                             ctgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATAATTGTGCATTTCacta                                                                                                                 <  1:684475‑M2/37‑1 (MQ=255)                               gcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCactact                                                                                                               >  2:625725‑M2/111‑149 (MQ=255)                                cttctcctcaccatcctgttattttcgccacctcacaatcatgggcttgctgaagaataattgaaatgatattattaattccactgcttttgggagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTg                                                                                                              <  1:657799‑M2/40‑1 (MQ=255)                                 ttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattcaactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCTTTTCACTACTgg                                                                                                             >  2:75268‑M2/109‑149 (MQ=255)                                   ctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGaa                                                                                                           <  1:328346‑M2/43‑1 (MQ=255)                                   ctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcatttggtagaggaaagtgctaaATAATAATCAATTTTTAAATTATTGTACATTTCACGACTTGaa                                                                                                           >  2:235064‑M2/107‑149 (MQ=255)                                      ctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTg                                                                                                        >  1:513707‑M2/104‑149 (MQ=255)                                       tcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGTAACTGt                                                                                                       >  2:173821‑M2/103‑149 (MQ=255)                                         accatcatgttattttcgccacataatcattttgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTaa                                                                                                     >  2:370989‑M2/101‑149 (MQ=255)                                          ccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAAt                                                                                                    <  1:677697‑M2/50‑1 (MQ=255)                                                          gccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACAt                                                                                    >  1:313174‑M2/84‑149 (MQ=255)                                                          gccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACAt                                                                                    >  2:517509‑M2/84‑149 (MQ=255)                                                           ccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATg                                                                                   >  1:746001‑M2/83‑149 (MQ=255)                                                            cacatcataatcctgggcttgctgaagaataattgagatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTACTCAGAAAAGATAGACATGc                                                                                  >  2:396586‑M2/82‑149 (MQ=255)                                                            cacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGc                                                                                  <  1:625725‑M2/68‑1 (MQ=255)                                                                     atcctgggcttgctgaagaataatggaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAAt                                                                         <  2:617315‑M2/77‑1 (MQ=255)                                                                           ggcttgctgaagaataattgaaattatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTAt                                                                   <  1:397512‑M2/83‑1 (MQ=255)                                                                            gcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTAtt                                                                  <  1:387547‑M2/84‑1 (MQ=255)                                                                            gcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTAtt                                                                  >  1:30570‑M2/66‑149 (MQ=255)                                                                             cttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTAttt                                                                 <  2:417775‑M2/85‑1 (MQ=255)                                                                                  tgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATttgattg                                                            >  1:241678‑M2/60‑149 (MQ=255)                                                                                          aattgaaatgatattattaattccactgcctttggtagaggaaagaaatccAAAGTAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGaaa                                                    <  1:23337‑M2/98‑1 (MQ=255)                                                                                            ttgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGa                                                  >  2:333454‑M2/50‑149 (MQ=255)                                                                                              gaaatgatattattaattccacagcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAAGTATTGTGCATTTCACTACTGGATCTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtg                                                <  1:173821‑M2/102‑1 (MQ=255)                                                                                                             attccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATgg                                 <  1:346509‑M2/117‑1 (MQ=255)                                                                                                                      cctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCGAATGTTTTGGTTGTTGTATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCt                        >  1:61916‑M2/24‑149 (MQ=255)                                                                                                                                 ggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGa             <  1:370989‑M2/137‑1 (MQ=255)                                                                                                                                       tgctaaATAATAATCAATTGGTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGt       <  2:311007‑M2/143‑1 (MQ=255)                                                                                                                                          taaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCCCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCaa    >  1:495769‑M2/4‑149 (MQ=255)                                                                                                                                            aATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATAGGTTCAATa  >  2:274653‑M2/2‑149 (MQ=255)                                                                                                                                            aATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCGATAGATTCTGATACGgtcacca  >  1:144462‑M2/2‑146 (MQ=255)                                                                                                                                             |                                                                                                                                                    TAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATA  >  NC_000913_3_bme_pgi/1300558‑1300845

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