New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4067814 =	26 (1.120)	3 (0.130)	3/244	NT	10.0%	coding (1294/3108 nt)	mdtC	multidrug transporter membrane protein
	?	NC_002947	4067854 =	29 (1.300)					coding (1254/3108 nt)	mdtC	multidrug transporter membrane protein

GCGACACCACCAGCGAGACGATGATCGCCGCCGCCAGGGTGATGGAAAACTCCTGGAACAGGTTGCGCACGATGCCGCCCATGAACAGGATGGAGACGAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATGGCATCGTCCACCACCAGCCCGGTGGCAAGGATCAGCGCCATCAGCGACAGGTTGTTGAGCGA  >  NC_002947/4067673‑4067961                                                                                                                                              |                                                                                                                                                    gCGACACCACCAGCGAGACGATGATCGCCGCCGCCAGGGTGATGGAAAACTCCTGGAACAGGTTGCGCACGATGCCGCCCATGAACAGGATGGAGACGAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCAcctcct                                                                                                                                             >  2:92958/1‑150 (MQ=255) gCGACACCACCAGCGAGACGATGATCGCCGCCGCCAGGGCGATGGAAAACTCCTGGAACAGGTTGCGCACGATGCCGCCCATGAACAGGATGGAGACGAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCAccgcct                                                                                                                                             >  2:11040/1‑150 (MQ=255)     caccacCAGCGAGACGATGATCGCCGCCGCCAGGGTGATGGAAAACTCCTGGAACAGGTTGCGCACGATGCCGCCCATGAACAGGATGGAGACGAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGgc                                                                                                                                         <  1:92958/150‑1 (MQ=255)        caccaGCGAGACGATGATCGCCGCCGCCAGGGTGATGGAAAACTCCTGGAACAGGTTGCGCACGATGCCGCCCATGAACAGGATGGAGACGAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACGCACCTCCTTGGCGcc                                                                                                                                      >  1:113514/1‑150 (MQ=255)                    gatgatCGCCGCCGCCAGGGTGATGGAAAACTCCTGGAACAGGTTGCGCACGATGCCGCCCATGAACAGGATGGAGACGAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGgccgc                                                                                                                          >  1:17846/1‑150 (MQ=255)                              gccgccAGGGTGATGGAAAACTCCTGGAACAGGTTGCGCACGATGCCGCCCATGAACAGGATGGAGACGAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGt                                                                                                                >  1:145213/1‑150 (MQ=255)                                  tcgctcaacaccctgtcgctgatggcgctgctccttgccaccgggctggtggtggacgatgccatcgtggtgctggaaaacatctccaggcatatcgagaacggccagCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGtggct                                                                                                            >  2:198517‑M1/109‑150 (MQ=255)                                           tGGAAAACTCCTGGAACAGGTTGCGCACGATGCCGCCCATGAACAGGATGGAGACGAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTtct                                                                                                   <  2:474713/150‑1 (MQ=255)                                                        gAACAGGTTCCGCACGATGCCGCCCATGAACAGGATGGAGACGAACACCGCCACCAGAGACACGTTCATCGACAGCAAGGTGAAACCAACCTCCGTCGGGCCCAGGTAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGga                                                                                      <  2:54808/150‑1 (MQ=255)                                                                 gcgcACGATGCCGCCCATGAACAGGATTGATACGAACACCTCCACCAGCGCCACGTTCATCGACAAGCGTTTTAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTc                                                                             <  1:11040/150‑1 (MQ=255)                                                                      cGATGCCGCCCATGAACAGGATGGAGACGAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGca                                                                        <  2:73740/150‑1 (MQ=255)                                                                       gATGCCGCCCATGAACAGGATGGAGACGAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGga                                                                                      >  1:477113/1‑135 (MQ=255)                                                                       gATGCCGCCCATGAACAGGATGGAGACGAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGga                                                                                      <  2:477113/135‑1 (MQ=255)                                                                                  gtggtggacgatgccatcgtggtgctggaaaacatctccaggcatatcgagaacggccagCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGc                                                                                           >  1:178286‑M1/61‑119 (MQ=255)                                                                                  gtggtggacgatgccatcgtggtgctggaaaacatctccaggcatatcgagaacggccagCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGc                                                                                           <  2:178286‑M1/59‑1 (MQ=255)                                                                                         gATGGAGACGAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATGGCATCGTccac                                                     >  2:373204/1‑150 (MQ=255)                                                                                            ggAGACGAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATGGCATCGTccaccac                                                  <  1:589753/150‑1 (MQ=255)                                                                                            ggAGACGAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATGGCATCGTccaccac                                                  <  2:17846/150‑1 (MQ=255)                                                                                                  gAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATGGCATCGt                                                         <  1:477642/137‑1 (MQ=255)                                                                                                  gAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATGGCATCGt                                                         >  2:477642/1‑137 (MQ=255)                                                                                                   aaCACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATGGCATCGTCCACCACCAGCCCg                                           <  1:373204/150‑1 (MQ=255)                                                                                                     cacCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATg                                                                >  1:141552/1‑127 (MQ=255)                                                                                                     cacCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATg                                                                <  2:141552/127‑1 (MQ=255)                                                                                                                    cacGTTCATCGACAGCAAGGTGAACCCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATGGCATCGTCCACCACCAGCCCGGTGGCAAGGATCAgcgc                          <  2:412607/150‑1 (MQ=255)                                                                                                                       gTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATGGCATCGTCCACCACCAGCCCGGTGGCaa                                    <  1:578720/137‑1 (MQ=255)                                                                                                                       gTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATGGCATCGTCCACCACCAGCCCGGTGGCaa                                    >  2:578720/1‑137 (MQ=255)                                                                                                                           aTCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTACTCGATATGCCTGGAGATGTTTTCCAGCACCACGATGGCATCGTCCACCACCAGCCCGGTGGCAAGGATCAGCGCCATCAGc                   <  1:46656/150‑1 (MQ=255)                                                                                                                                  gCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATGGCATCGTCCACCACCAGCCCGGTGGCAAGGATCAGCGCCATCAGCGACAGgt            >  1:349885/1‑150 (MQ=255)                                                                                                                                      ggTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATGGCATCGTCCACCACCAGCCCGGTGGCAAGGATCAGCGCCATCAGCGACAGgttgtt        <  2:349885/150‑1 (MQ=255)                                                                                                                                            aCCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATGGCATCGTCCACCACCAGCCCGGTGGCAAGGATCAGCGCCATCAGCGACAGGTTGTTGAGCga  <  2:113514/150‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                    GCGACACCACCAGCGAGACGATGATCGCCGCCGCCAGGGTGATGGAAAACTCCTGGAACAGGTTGCGCACGATGCCGCCCATGAACAGGATGGAGACGAACACCGCCACCAGCGACACGTTCATCGACAGCAAGGTGAAACCCACCTCCTTGGCGCCCAGGAAGGCCGCCTTCATCGGTGGCTGGCCGTTCTCGATATGCCTGGAGATGTTTTCCAGCACCACGATGGCATCGTCCACCACCAGCCCGGTGGCAAGGATCAGCGCCATCAGCGACAGGTTGTTGAGCGA  >  NC_002947/4067673‑4067961

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.