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

GTGCCAGCCCAGGCCCTGCATCTCCAGGCCTGCCGCTTCCGGCTCGGGGCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAGGCCCATCTGCACGGCGGCCAGCGGGTTTTCCAGGTTACGCTCGGCCGACAGGTCGCGATTCTGC  >  NC_002947/4168989‑4169277                                                                                                                                                     |                                                                                                                                             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)                          tggaaaacccgctggccgccgtgcagatgggcctgatctatgtgaacccggaaggcccggaaggtaaccccgaccccgtggcctcgggcaaggacatccgtgaaacctttggccatgcggccaAAGGTTTCACGGATGTCCTTGCCCGAg                                                                                                                    <  1:175905‑M2/27‑1 (MQ=255)                              cTGCCGCTTCCGGCTCGGGGCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTc                                                                                                                                       >  1:714589/1‑127 (MQ=255)                              cTGCCGCTTCCGGCTCGGGGCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTc                                                                                                                                       <  2:714589/127‑1 (MQ=255)                                                ggCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCAcggggtcgg                                                                                                       <  1:676289/141‑1 (MQ=255)                                                ggCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCAcggggtcgg                                                                                                       >  2:676289/1‑141 (MQ=255)                                                    gacaacgTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTAccttccgg                                                                                          <  2:35666/144‑1 (MQ=255)                                                             ggcggGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGGCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTAccttccgggccttccgg                                                                                 <  1:814752/150‑1 (MQ=255)                                                                   gCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTcac                                                                           <  2:748402/150‑1 (MQ=255)                                                                       gCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGt                                                                               <  1:603261/142‑1 (MQ=255)                                                                       gCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGt                                                                               >  2:603261/1‑142 (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                                                                                                        <  2:463101/58‑1 (MQ=255)                                                                                                                                  gTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCAcggggtcg                                                                                                        >  1:463101/1‑58 (MQ=255)                                                                                                                                  acctttggccgcatggccaAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTAcct                                                                                               <  2:230726‑M2/48‑1 (MQ=255)                                                                                                                                  acctttggccgcatggccaAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTAcct                                                                                               >  1:230726‑M2/20‑67 (MQ=255)                                                                                                                                            aTGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCAcggggtcgggg                                                                                                     >  2:376332/1‑51 (MQ=255)                                                                                                                                            aTGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCAcggggtcgggg                                                                                                     <  1:376332/51‑1 (MQ=255)                                                                                                                                            aTGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAGGCCCATCTGCACGGCGGCCAGCGGGTTTTCCAGGTTACGCTCGGCCGACAGGTCGCGATTCTGc  >  1:621267/1‑150 (MQ=255)                                                                                                                                                   caAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAGGCCCATc                                                           <  2:407472‑M2/84‑1 (MQ=255)                                                                                                                                                   caAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAGGCCCATc                                                           >  1:407472‑M2/3‑86 (MQ=255)                                                                                                                                                     |                                                                                                                                             GTGCCAGCCCAGGCCCTGCATCTCCAGGCCTGCCGCTTCCGGCTCGGGGCCGACGTTGTCGGCGGGGCCGGCACCATGGGTCTTGCCGAACGCGTGGCCACCGGCGATCAAGGCCACGGTTTCTTCATCGTTCATGGCCATGCGGCCAAAGGTTTCACGGATGTCCTTGCCCGAGGCCACGGGGTCGGGGTTACCTTCCGGGCCTTCCGGGTTCACATAGATCAGGCCCATCTGCACGGCGGCCAGCGGGTTTTCCAGGTTACGCTCGGCCGACAGGTCGCGATTCTGC  >  NC_002947/4168989‑4169277

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.