New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	1846333 =	7 (0.400)	4 (0.250)	4/286	NT	33.9%	coding (190/642 nt)	pncA	nicotinamidase/pyrazinamidase
	?	CP009273	1846630 =	9 (0.550)					coding (487/642 nt)	pncA	nicotinamidase/pyrazinamidase

GGATGATTAAGGAGACTGTAATGCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAG  >  CP009273/1846124‑1846505                                                                                                                                                                                                                  |                                                                                                                                                                             ggATGATTAAGGATGCTTTAATGCCCCCGCGCGCCCTGTTACTGGTCGATTTACAAAATGATTGCGGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGa                                                                                                              <  1:173947/274‑1 (MQ=255)                                                                                         agagacaGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATc                     <  1:125444/269‑1 (MQ=255)                                                                                            gggaGGTGTCAGTACGGTGGATGTCGCTAATCGCCTGATTGACTGGTGCCTGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACACCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCtt                  <  1:105714/272‑1 (MQ=255)                                                                                                    aCAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCTCGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTGGAGCCTTATACGCCAGGCCAACTCGACAGTTTGCCCCAAACCTTCTGGCCAGACCACTGTGTGCAGAACAGTCAAGGCCCACAATTACATCCGTTACTTCACCAAAAAGCGACCGCAGCGGTGTTCCATAAAGGCGAAAATACATTAGTTTAc          >  1:37901/1‑275 (MQ=255)                                                                                                            gTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACATTCTGGCACGAGCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCAGTTACTGACACAAAAAGCGAACGCCGCGGAGGTCCATAAAGGCGAAAATCCATTAATTGNcctgtgccg  >  2:125442/1‑267 (MQ=255)                                                                                                                atctccataaacgcgtgcgcactgtcctggggctggatattcacgccacgacaaccatcggtaatcacgtttaccttataacctaactgtaacgcgtcCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACtgtg                                                                                                      <  2:95140‑M1/73‑1 (MQ=255)                                                                                                                    ccataaacgcgtgcgcactgtcctggggctggatattcacgccacgacaaccatcggtaatcacgtttaccttataacctaactgtaacgcgtcCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCa                                                                                                  >  1:95142‑M1/95‑171 (MQ=255)                                                                                                                                                          acgccacgacaaccatcggtaatcacgtttaccttataacctaactgtaacgcgtcCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGAt                                                                                                             <  1:90760‑M1/66‑1 (MQ=255)                                                                                                                                                              cacgacaaccatcggtaatcacgtttaccttataacctaactgtaacgcgtcCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACt                                                                                                         >  2:90758‑M1/53‑122 (MQ=255)                                                                                                                                                                  aTCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAgg                                                                                      >  1:172719/1‑137 (MQ=255)                                                                                                                                                                      ccAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGc                                                                                                                                                                           <  2:21898/48‑1 (MQ=255)                                                                                                                                                                      ccAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGca                                                                                  <  2:172715/137‑1 (MQ=255)                                                                                                                                                                                                                  |                                                                                                                                                                             GGATGATTAAGGAGACTGTAATGCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAG  >  CP009273/1846124‑1846505

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