New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4169123 =	28 (0.880)	3 (0.100)	3/238	NT	9.9%	coding (814/2256 nt)	katG	catalase‑peroxidase
	?	NC_002947	4169137 =	28 (0.920)					coding (800/2256 nt)	katG	catalase‑peroxidase

AAACGTGTTGTGCCAGCCCAGGCCCTGCATCTCCAGGCCTGCCGCTTCCGGCTCGGGGCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAGGCCCATCTGCACGGCGGCCAGCGGGTTTTCCAGGTTACGCTCGGCCGACA  >  NC_002947/4168980‑4169263                                                                                                                                                |                                                                                                                                             aaaCGTGTTGTGCCAGCCCAGGCCCTGCATCTCCAGGCCTGCCGCTTCCGGCTCGGGGCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGccgc                                                                                                                                        >  1:53933/1‑148 (MQ=255)     gtgtTGTGCGAGCGCAGGCCGTGCATTTCCAGGCCTGCCGCTTCCGGCTCGGGGCCGACGTTGTCGGCGGGGCCGGCACCATGGATCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCgg                                                                                                                                    <  2:432893/150‑1 (MQ=255)          gtgCCAGCCCAGGCCCTGCATCTCCAGGCCTGCCGCTTCCGGCTCGGGGCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCaaa                                                                                                                               <  1:244750/150‑1 (MQ=255)               aGCCCAGGCCCTGCATCTCCAGGCCTGCCGCTTCCGGCTCGGGGCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGttt                                                                                                                          <  2:164680/150‑1 (MQ=255)                  ccAGGCCCTGCATCTCCAGGCCTGCCGCTTCCGGCTCGGGGCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCAAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCAc                                                                                                                       <  2:483507/150‑1 (MQ=255)                   cAGGCCCTGCATCTCCAGGCCTGCCGCTTCCGGCTCGGGGCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGTTCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACg                                                                                                                      >  2:842418/1‑150 (MQ=255)                                       cTGCCGCTTCCGGCTCGGGGCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTc                                                                                                                         <  2:714589/127‑1 (MQ=255)                                       cTGCCGCTTCCGGCTCGGGGCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTc                                                                                                                         >  1:714589/1‑127 (MQ=255)                                                         ggCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCAcggggtcgg                                                                                         >  2:676289/1‑141 (MQ=255)                                                         ggCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCAcggggtcgg                                                                                         <  1:676289/141‑1 (MQ=255)                                                            gatgggcctgatctatgtgaacccggaaggcccggaaggtaaccccgaccccgtggcctcgggcaaggacatccgtgaaaccttTg                                                                                                                                             >  1:407472‑M1/85‑86 (MQ=255)                                                            gatgggcctgatctatgtgaacccggaaggcccggaaggtaaccccgaccccgtggcctcgggcaaggacatccgtgaaaccttTg                                                                                                                                             <  2:407472‑M1/2‑1 (MQ=255)                                                             gacaacgTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTAccttccgg                                                                            <  2:35666/144‑1 (MQ=255)                                                                      ggcggGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGGCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTAccttccgggccttccgg                                                                   <  1:814752/150‑1 (MQ=255)                                                                            gCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTcac                                                             <  2:748402/150‑1 (MQ=255)                                                                                gCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGt                                                                 >  2:603261/1‑142 (MQ=255)                                                                                gCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGt                                                                 <  1:603261/142‑1 (MQ=255)                                                                                    cATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCa                                                     >  2:352731/1‑150 (MQ=255)                                                                                      tGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAgg                                                   >  1:457443/1‑150 (MQ=255)                                                                                          tCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAGGCCCa                                               >  1:517757/1‑150 (MQ=255)                                                                                                aggtaaccccgaccccgtggcctcgggcaaggacatccgtgaaaccttTGGCCATGCGGCCAAAGGt                                                                                                                            <  2:230726‑M1/19‑1 (MQ=255)                                                                                                aggtaaccccgaccccgtggcctcgggcaaggacatccgtgaaaccttTGGCCATGCGGCCAAAGGt                                                                                                                            >  1:230726‑M1/49‑67 (MQ=255)                                                                                                 gAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGtt                                                                                     >  1:837164/1‑105 (MQ=255)                                                                                                 gAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGtt                                                                                     <  2:837164/105‑1 (MQ=255)                                                                                                       gTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTAccttccgggccttccggg                                                                  <  1:600859/118‑1 (MQ=255)                                                                                                       gTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTAccttccgggccttccggg                                                                  >  2:600859/1‑118 (MQ=255)                                                                                                             aCCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAGGCCCATCTGCACGGCGGCCAGCgg                            <  1:352731/150‑1 (MQ=255)                                                                                                                 gCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAGGCCCATCTGCACGGCGGCCAGCGGGttt                        <  2:752012/150‑1 (MQ=255)                                                                                                                     ctcgggcaaggacatccgtgaaaccttTGGCCGCATGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAGGCCCATCTGCACGGCGGCCAGCGGGTTTTCCa                    <  1:175905‑M1/123‑1 (MQ=255)                                                                                                                       aaGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAGGCCCATCTGCACGGCGGCCAGCGGGTTTTCCAgg                  <  1:48280/150‑1 (MQ=255)                                                                                                                               ggTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAGGCCCATCTGCACGGCGGCCAGCGGGTTTTCCAGGTTACGCTc          >  1:272655/1‑150 (MQ=255)                                                                                                                                     ttcATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTGGCCCGAGGCCACGGGGTCGGGGTTACCTTCCCCGGAAACGGGCTTGACATAGAGCAGGGCCATCTTAATGCCGGCCAGCGCGGTTTCCTGGTTACGCTCGGCCGa    >  1:232630/1‑150 (MQ=255)                                                                                                                                       cATCGTTCATGGCCATGCGGCCAGACATAAGGCCTATGGCCGTGCGCGGTGCCACGTGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAGGCCCATCTGCACGGCGGCCAGCGGGTTTTCCAGGTTACGCTCGGCCGACa  <  2:520954/150‑1 (MQ=255)                                                                                                                                       cATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAGGCCCATCTGCACGGCGGCCAGCGGGTTTTCCAGGTTACGCTCGGCCGACa  <  2:457443/150‑1 (MQ=255)                                                                                                                                           gTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCAcggggtcg                                                                                          >  1:463101/1‑58 (MQ=255)                                                                                                                                           gTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCAcggggtcg                                                                                          <  2:463101/58‑1 (MQ=255)                                                                                                                                                |                                                                                                                                             AAACGTGTTGTGCCAGCCCAGGCCCTGCATCTCCAGGCCTGCCGCTTCCGGCTCGGGGCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAGGCCCATCTGCACGGCGGCCAGCGGGTTTTCCAGGTTACGCTCGGCCGACA  >  NC_002947/4168980‑4169263

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ 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.