New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 3916335	16 (1.030)	4 (0.270)	4/232	1.4	21.1%	coding (1041/1383 nt)	atpD	F1 sector of membrane‑bound ATP synthase, beta subunit
	?	NC_000913_3_pae_tpiA	= 3916349	15 (1.030)					coding (1027/1383 nt)	atpD	F1 sector of membrane‑bound ATP synthase, beta subunit

TCAGACAGTTCATCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTCCTGACCAACCACCAGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCGGGTCAACGGCCGGGTAGATACCCAGAGACGCGATCTGACGGCTCAGTACCACGGTTGCGTCAAGGTGCGCAAAGGTGGTTGCCGGAGACGGGTCAGTCAAGTCATCCGCAGGTACGTATACTGCC  >  NC_000913_3_pae_tpiA/3916222‑3916491                                                                                                                                |                                                                                                                                               tcAGACAGTTCATCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTCCTGACCAACCACCAGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCGGGTCa                                                                                                                          >  1:182043/1‑150 (MQ=255)              ccATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTCCTGACCAACCACCAGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCGGGTCAACGGCCGGGTAGa                                                                                                             <  2:182043/150‑1 (MQ=255)                                       tttCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTCCTGACCAACCACCAGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCGGGTCAACGGCCGGGTAGATACCCAGAGACGCGATCTGACGGCt                                                                                    >  2:44822/1‑150 (MQ=255)                                             ttCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTCCTGACCAACCACCAGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCGGGTCAACGGCCGGGTAGATACCCAGAGACGCGATCTGACGGCTCAGTAc                                                                              >  3:95543/1‑150 (MQ=255)                                                 tGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTCCTGACCAACCACCAGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCGGGTCAACGGCCGGGTAGATACCCAGAGACGCGATCTGACGGCTCAGTACCACg                                                                          >  2:3229/1‑150 (MQ=255)                                                        gTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTCCTGACCAACCACCAGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCGGGTCAACGGCCGGGTAGATACCCAGAGACGCGATCTGACGGCTCAGTACCACGGTTGCGt                                                                   >  4:57888/1‑150 (MQ=255)                                                                     cTGAACGCCACGCGCGGTGTCGTAGTGTTCCTGACCAACCACCAGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCggg                                                                                                                             <  4:7660/79‑1 (MQ=255)                                                                     cTGAACGCCACGCGCGGTGTCGTAGTGTTCCTGACCAACCACCAGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCggg                                                                                                                             >  3:7660/1‑79 (MQ=255)                                                                                  gcgGTGTCGTAGTGTTCCTGACCAACCACCAGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCGGGTCAACGGCCGGGTAGATACCCAGAGACGCGATCTGACGGCTCAGTACCACGGTTGCGTCAAGGTGCGCAAAGGTGGTTGCCGga                                         <  3:57888/150‑1 (MQ=255)                                                                                                ttCCTGACCAACCACCAGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCGGGTCAACGGCCGGGTAGATACCCAGAGACGCGATCTGACGGCTCAGTACCACGGTTGCGTCAAGGTGCGCAAAGGTGGTTGCCGGAGACGGGTCagtcaa                           <  1:44822/150‑1 (MQ=255)                                                                                                 tCCTGACCAACCACCAGCGGGTCCAGCTGACGgctggtggttggtcaggaacactacgacaccgcgcgtggcgttcagtccatcctgcaacgttatc                                                                               <  3:101867‑M2/97‑66 (MQ=255)                                                                                                 tCCTGACCAACCACCAGCGGGTCCAGCTGACGgctggtggttggtcaggaacactacgacaccgcgcgtggcgttcagtccatcctgcaacgttatc                                                                               >  4:101867‑M2/1‑32 (MQ=255)                                                                                                 tCCTGACCAACCACCAGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCGGGTCAACGGCCGGGTAGATACCCAGAGACGCGATCTGACGGCTCAGTACCACGGTTGCGTCAAGGTGCGCAAAGGTGGTTGCCGGAGACGGGTCagtcaag                          <  1:3229/150‑1 (MQ=255)                                                                                                  ccTGACCAACCACCAGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCGGGTCAACGGCCGGGTAGATACCCAGAGACGCGATCTGACGGCTCAGTACCACGGTTGCGTCAAGGTGCGCAAAGGTGGTTGCCGGAGACGGGTCagtcaagt                         >  2:26195/1‑150 (MQ=255)                                                                                                    tGACCAACCACCAGCGGGTCCAGCTGACGgctggtggttggtcaggaacactacgacaccgcgcgtggcgttcagtccatcctgcaacgttatcaggaactgaaagacatcatcgccatcctgg                                                 >  2:20855‑M2/1‑29 (MQ=255)                                                                                                    tGACCAACCACCAGCGGGTCCAGCTGACGgctggtggttggtcaggaacactacgacaccgcgcgtggcgttcagtccatcctgcaacgttatcaggaactgaaagacatcatcgccatcctgg                                                 <  1:20855‑M2/124‑96 (MQ=255)                                                                                                              ccaGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCGGGTCAACGGCCGGGTAGATACCCAGAGACGCGATCTGACGGCTCAGTACCACGGTTGCGTCAAGGTGCGCAAAGGTGGTTGCCGGAGACgg                                    >  2:84765/1‑127 (MQ=255)                                                                                                              ccaGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCGGGTCAACGGCCGGGTAGATACCCAGAGACGCGATCTGACGGCTCAGTACCACGGTTGCGTCAAGGTGCGCAAAGGTGGTTGCCGGAGACgg                                    <  1:84765/127‑1 (MQ=255)                                                                                                                         aGCTGACGGCTGGTGGAGTCCAGCGGGTCAACGGCCGGGTAGATACCCAGAGACGCGATCTGACGGCTCAGTACCACGGTTGCGTCAAGGTGCGCAAAGGTGGTTGCCGGAGACGGGTCAGTCAAGTCATCCGCAGGTACGTATACTGcc  >  4:3035/1‑150 (MQ=255)                                                                                                                                |                                                                                                                                               TCAGACAGTTCATCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTCCTGACCAACCACCAGCGGGTCCAGCTGACGGCTGGTGGAGTCCAGCGGGTCAACGGCCGGGTAGATACCCAGAGACGCGATCTGACGGCTCAGTACCACGGTTGCGTCAAGGTGCGCAAAGGTGGTTGCCGGAGACGGGTCAGTCAAGTCATCCGCAGGTACGTATACTGCC  >  NC_000913_3_pae_tpiA/3916222‑3916491

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.