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

TGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGC  >  NC_000913_3_pae_tpiA/1792804‑1793094                                                                                                                                                   |                                                                                                                                                 tGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGcc                                                                                                                                               <  1:21083/150‑1 (MQ=255)         gATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGCAATTTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGTACCTCTTCCCTGCGTCCACTGGCCAGCATCTTTCCACCTTTTTGCAACCACGCCCGATgcgc                                                                                                                                       >  3:90661/1‑150 (MQ=255)                     cAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAAtgt                                                                                                                           <  3:15232/150‑1 (MQ=255)                                gACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACacattgcgtcatgcgcatcgggcgtggttgctaa                                                                                                                <  4:80639‑M1/150‑35 (MQ=255)                                 aCCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATc                                                                                                               >  3:23842/1‑150 (MQ=255)                                     tCGATAGTCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTTCGTCCACTGGCAAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGa                                                                                                           <  4:90661/150‑1 (MQ=255)                                                  ggCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACaa                                                                                              <  4:23842/150‑1 (MQ=255)                                                            cATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACacattgcgtcatgcgcatcgggcgtggttgct                                                                                                                  <  2:109084‑M1/120‑33 (MQ=255)                                                            cATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACacattgcgtcatgcgcatcgggcgtggttgct                                                                                                                  >  1:109084‑M1/1‑88 (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)                                                                                                                                                   |                                                                                                                                                 TGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGC  >  NC_000913_3_pae_tpiA/1792804‑1793094

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.