New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 4502020	17 (1.360)	4 (0.350)	4/238	1.1	20.8%	coding (1328/1329 nt)	insG	IS4 transposase
	?	NC_000913_3_pae_tpiA	= 4502024	15 (1.310)					coding (1324/1329 nt)	insG	IS4 transposase

GGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATA  >  NC_000913_3_pae_tpiA/4501872‑4502163                                                                                                                                                     |                                                                                                                                                gggCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTa                                                                                                                                                >  2:135087/1‑150 (MQ=255)   tcATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAg                                                                                                                                              >  4:25816/2‑150 (MQ=255)                                tgGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCtgtg                                                                                                                 <  1:135087/150‑1 (MQ=255)                                  gCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTggg                                                                                                               <  3:27907/150‑1 (MQ=255)                                               gTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCCACTGCCTGGGTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGc                                                                                                  <  3:90791/150‑1 (MQ=255)                                                                    aGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTgcttaactgactggcattacagaaccaagat                                                                                                                  <  1:18206‑M1/113‑32 (MQ=255)                                                                    aGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTgcttaactgactggcattacagaaccaagat                                                                                                                  >  2:18206‑M1/1‑82 (MQ=255)                                                                        ttCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGaaag                                                                         >  2:110382/1‑148 (MQ=255)                                                                                 gCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTCCCACTCTCGGGAAGGCCCTTTccc                                                                <  2:68566/150‑1 (MQ=255)                                                                                 gCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTccc                                                                <  3:85626/150‑1 (MQ=255)                                                                                      cAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCcttctt                                                           <  4:48964/150‑1 (MQ=255)                                                                                      cAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCcttctt                                                           <  3:25816/150‑1 (MQ=255)                                                                                        gATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCtctc                                                                                             >  2:40668/1‑114 (MQ=255)                                                                                        gATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCtctc                                                                                             <  1:40668/114‑1 (MQ=255)                                                                                           aaTTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTa                                                                                        <  2:121595/116‑1 (MQ=255)                                                                                           aaTTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTa                                                                                        >  1:121595/1‑116 (MQ=255)                                                                                                             gAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTgcttaactgactggcatt                                                                                                                               <  2:123641‑M1/59‑19 (MQ=255)                                                                                                             gAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTgcttaactgactggcatt                                                                                                                               >  1:123641‑M1/1‑41 (MQ=255)                                                                                                              aaaCGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCttt                                                                                         >  2:18575/1‑96 (MQ=255)                                                                                                              aaaCGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCttt                                                                                         <  1:18575/96‑1 (MQ=255)                                                                                                                           aTCTTGGTTCTGTAATGCCAGTCAGTTAAGCaactgactggcattacagaaccaagatgcctatttcgtttcagcatcaactgttaaattatctggttgccaacgagaatact                                                           >  2:18206‑M2/1‑31 (MQ=255)                                                                                                                           aTCTTGGTTCTGTAATGCCAGTCAGTTAAGCaactgactggcattacagaaccaagatgcctatttcgtttcagcatcaactgttaaattatctggttgccaacgagaatact                                                           <  1:18206‑M2/113‑83 (MQ=255)                                                                                                                                        aaTGCCAGTCAGTTAAGCaactgactggcattacagaaccaagatgcctatttcgtttc                                                                                                    <  2:123641‑M2/59‑42 (MQ=255)                                                                                                                                        aaTGCCAGTCAGTTAAGCaactgactggcattacagaaccaagatgcctatttcgtttc                                                                                                    >  1:123641‑M2/1‑18 (MQ=255)                                                                                                                                               gtcagtTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGtata  <  2:73784/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                GGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATA  >  NC_000913_3_pae_tpiA/4501872‑4502163

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