Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	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	1299499	1300697–1299500	2–1199	40 [0]	[2] 43	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	= 1299498	0 (0.000)	40 (1.050)	38/276	0.1	98.8%	intergenic (+253/‑1684)	ychE/oppA	UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein
	?	NC_000913	1300698 =	1 (0.030)					intergenic (+1453/‑484)	ychE/oppA	UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

GATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGT  >  NC_000913/1300556‑1300840                                                                                                                                               |                                                                                                                                               cattccaggtagttatgaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaa                                                                                                                                                <  2:450011‑M2/1‑1 (MQ=255)   ttacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaata                                                                                                                                              >  2:233366‑M2/141‑143 (MQ=255)      cGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCACACATGAAGAAATGAAATGACAGAGCCAGCCGAGAAGAATTTACCCGCTTAATCGCACCTTCCCTAaataatc                                                                                                                                           >  2:356151/1‑142 (MQ=11)             ttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAAtt                                                                                                                                    <  2:228127‑M2/13‑1 (MQ=255)               ataaggatttccttaactgcttctcctcaccatcatgtgattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGt                                                                                                                                  >  2:350291‑M2/129‑143 (MQ=255)                taaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaa                                                                                                                                                >  1:242301‑M2/128‑128 (MQ=255)                taaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaa                                                                                                                                                <  2:242301‑M2/1‑1 (MQ=255)                        tccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAattatt                                                                                                                         <  1:658417‑M2/24‑1 (MQ=255)                             aactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttaatgaagaataattgaaataatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCACTTGTTAAATTATTGTGCa                                                                                                                    >  1:451054‑M2/115‑143 (MQ=255)                                tgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATGATTGTGCAttt                                                                                                                 >  1:157346‑M2/112‑143 (MQ=255)                                  cttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCa                                                                                                               <  2:309625‑M2/34‑1 (MQ=255)                                      tcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCactac                                                                                                           >  2:462532‑M2/106‑143 (MQ=255)                                                   gttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTGCTGGAACTGTAATc                                                                                              >  1:273947‑M2/93‑143 (MQ=255)                                                         ttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGaaaa                                                                                        <  2:333634‑M2/57‑1 (MQ=255)                                                          tcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTAATGGAACTGTAATCAGAagag                                                                                       >  2:423576‑M2/86‑143 (MQ=255)                                                           cgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGa                                                                                      >  1:645889‑M2/85‑143 (MQ=255)                                                                   cataatcctgcgctttctgaagaataattgaaatgatattattaatgccactgactttggtagaggaaagagctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATg                                                                              <  1:141174‑M2/67‑1 (MQ=255)                                                                       atcctgggcttgctgaagaagaattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTa                                                                          <  1:272360‑M2/71‑1 (MQ=255)                                                                          ctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGcc                                                                       >  2:381272‑M2/70‑143 (MQ=255)                                                                           tgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCa                                                                      >  1:596788‑M2/69‑143 (MQ=255)                                                                           tgggcttgctgaagaataattgaaacgatattattaattccactgtctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCa                                                                      <  1:467031‑M2/75‑1 (MQ=255)                                                                            gggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCaa                                                                     <  1:344683‑M2/76‑1 (MQ=255)                                                                            gggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCaa                                                                     >  1:670521‑M2/68‑143 (MQ=255)                                                                             ggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGTAACTGTAATCAGAAAAGATAGACATGCTTAGCCAAt                                                                    <  1:462532‑M2/77‑1 (MQ=255)                                                                                     gaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTAttt                                                            >  1:576746‑M2/59‑143 (MQ=255)                                                                                        gaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATttgat                                                         <  1:9701‑M2/88‑1 (MQ=255)                                                                                        gaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATttgat                                                         >  1:475333‑M2/56‑143 (MQ=255)                                                                                                  aatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGaaa                                               <  1:713812‑M2/98‑1 (MQ=255)                                                                                                         ttattaattccactggttttggtatatgaaagtgctaaATAATAATAAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgttt                                        <  2:25224‑M2/105‑1 (MQ=255)                                                                                                          tattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtttg                                       >  1:172250‑M2/38‑143 (MQ=255)                                                                                                            ttaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAggtttgtt                                     >  2:548284‑M2/36‑143 (MQ=255)                                                                                                            ttaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTCCACTACTGGAACTGTATTCAGACAAGATAGACATCCTTAGCCAGTCTCTATTTGATAGAACAGACAGAagtttttt                                     >  2:242562‑M2/36‑143 (MQ=255)                                                                                                                   cactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCAt                              >  2:69465‑M2/29‑143 (MQ=255)                                                                                                                    actgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATg                             >  2:183224‑M2/28‑143 (MQ=255)                                                                                                                    actgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAAACAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATg                             <  1:366250‑M2/116‑1 (MQ=255)                                                                                                                      tgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCTCTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGa                           >  2:356821‑M2/26‑143 (MQ=255)                                                                                                                       gcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGAt                          <  1:350291‑M2/119‑1 (MQ=255)                                                                                                                                agatgaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAATATAGACATGCTTAGCCAATCTCTATTCGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCtat                 <  1:235174‑M2/128‑1 (MQ=255)                                                                                                                                    gaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATCGAt             >  2:514539‑M2/12‑143 (MQ=255)                                                                                                                                     aaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGATAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGAtt            <  2:576746‑M2/133‑1 (MQ=255)                                                                                                                                       agtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACAACCGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCt          >  1:303523‑M2/9‑143 (MQ=255)                                                                                                                                              aataaccatCAATGGTTAACGTATGGTGCATGTCACTACGGGAGCCGAAATCAGAAAAGATAGACATCTAGAGCCAAGCTCTATAGGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACgg   <  2:209183/143‑1 (MQ=255)                                                                                                                                               ataataatCAATTGTCAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGt  <  2:157346/143‑1 (MQ=255)                                                                                                                                               |                                                                                                                                               GATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGT  >  NC_000913/1300556‑1300840

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