New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 1792950	11 (0.880)	3 (0.250)	3/246	1.3	20.4%	coding (609/750 nt)	btuD	vitamin B12 ABC transporter ATPase
	?	NC_000913_3_pae_tpiA	= 1792973	13 (1.100)					coding (586/750 nt)	btuD	vitamin B12 ABC transporter ATPase

CAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTG  >  NC_000913_3_pae_tpiA/1792824‑1793122                                                                                                                                                      |                                                                                                                                                      cAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAAtgt                                                                                                                                                       <  3:15232/150‑1 (MQ=255)             aCCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATc                                                                                                                                           >  3:23842/1‑150 (MQ=255)                 tCGATAGTCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTTCGTCCACTGGCAAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGa                                                                                                                                       <  4:90661/150‑1 (MQ=255)                              ggCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACaa                                                                                                                          <  4:23842/150‑1 (MQ=255)                                                   ccTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGAcaca                                                                                                     >  1:145571/1‑150 (MQ=255)                                                     tGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAgc                                                                                                   >  4:79203/1‑150 (MQ=255)                                                                                                    cAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACt                                                    >  2:32248/1‑150 (MQ=255)                                                                                                             tCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCgg                                           <  2:132252/150‑1 (MQ=255)                                                                                                                     ttAGCAACCACGCCCGATGCGCATGACGCAATGTgtggttgctaaaaggtggaaaaatgctggccagtggacgcagggaagaggtgctcacgccgccaaatctggcgcaggcctatgggatgaattttcgccgtctggatatcgaaggtc                                   <  4:80639‑M2/150‑117 (MQ=255)                                                                                                                       aGCAACCACGCCCGATGCGCATGACGCAATGTgtggttgctaaaaggtggaaaaatgctggccagtggacgcagggaagaggtgctcacgccgccaaatctggcgcaggcctatgggatg                                                               <  2:109084‑M2/120‑89 (MQ=255)                                                                                                                       aGCAACCACGCCCGATGCGCATGACGCAATGTgtggttgctaaaaggtggaaaaatgctggccagtggacgcagggaagaggtgctcacgccgccaaatctggcgcaggcctatgggatg                                                               >  1:109084‑M2/1‑32 (MQ=255)                                                                                                                          aaCCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAagc                              <  3:79203/150‑1 (MQ=255)                                                                                                                                   cGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCa                                                               >  3:24794/1‑108 (MQ=255)                                                                                                                                   cGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCa                                                               <  4:24794/108‑1 (MQ=255)                                                                                                                                    gATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGc                    >  4:20340/1‑150 (MQ=255)                                                                                                                                                      tgtgGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTtgtg  <  2:109799/150‑1 (MQ=255)                                                                                                                                                      tgtgGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTtgtg  >  2:101954/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTG  >  NC_000913_3_pae_tpiA/1792824‑1793122

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.