New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 4609514	33 (0.910)	4 (0.120)	4/270	4.6	11.7%	coding (185/1590 nt)	prfC	peptide chain release factor RF‑3
	?	NC_000913_3_pae_tpiA	= 4609530	29 (0.840)					coding (201/1590 nt)	prfC	peptide chain release factor RF‑3
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCATTATTTCTCACCCGGACGCCGGTAAGACTACCATCACCGAGAAGGTGCTGCTGTTCGGACAGGCCATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATG  >  NC_000913_3_pae_tpiA/4609379‑4609662                                                                                                                                        |                                                                                                                                                     ccATCATTTCTCACCCGGACGCCGGTAAGACAACCATCACCGAGAAGGTGAAGCTTGTCGCACAGGCCATTCAGACCGCCGGGACAGGAAAAGGCCGAAGCTCCACCCAGCACGCGACGTCGGAGTGGCCGTAGCTGGAAAAGCAgtgt                                                                                                                                         >  2:165097/1‑149 (MQ=255)   attattTCTCACCCGGACGCCGGTAAGACTACCATCACCGAGAAGGTGCTGCTGTTCGGACAGGCCATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTgg                                                                                                                                       >  1:148149/1‑149 (MQ=255)    ttattTCTCACCCGGACGCCAGTAAGACTACCATCACCGAGAAGGGGCTGCGGCTGGGACATAAAGTGCAGAGCACCGGGAAGGTAAAGGGACGTGGGGACACCGAGCACACTAATTCGGACTCGATGGCGATGGACAAGCAGCGTggg                                                                                                                                      >  2:510074/1‑149 (MQ=255)    ttagtTTTCACCCGGACGGCGGGAAGACTACCATCACCGAGAAGGTGCTGCTGTTCTGAAAGGTCATTCAGAACGCCGGTACAGTAAAGGGCCGTGGTTGCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTggg                                                                                                                                      <  1:205233/149‑1 (MQ=255)    gtggttCTCACCCGGACGTCGGTAAGACTACCAGCACCGAGAAGGTGCTGCTGTTGGGACAGGCCATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGGACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTggg                                                                                                                                      <  2:259286/145‑1 (MQ=255)             aCCCGGACGCCGGTAAGACTACCATCACCGAGAAGGTGCTGCTGTTCGGACAGGCCATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCAtt                                                                                                                             >  1:400235/1‑149 (MQ=255)               ccGGACGCCGGTAAGACTACCATCACCGAGAAGGTGCTGCTGTTCGGACAGGCCATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATtac                                                                                                                           >  1:77863/1‑149 (MQ=255)                 ggACGCCGGTAAGACTACCATCACCGAGAAGGTGCTGCTGTTCGGACAGGCCATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATtacta                                                                                                                         <  2:94410/149‑1 (MQ=255)                        ggTAAGACTCCCCGCACCGAGAAGGTGCTGCTGTTCGGACAGGCCATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTTGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCtgt                                                                                                                  <  2:284053/149‑1 (MQ=255)                         gTAAGACTACCATCACCGAGAAGGTGCTGCTGTTCGGACAGGCCATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCtgtg                                                                                                                 <  1:59803/149‑1 (MQ=255)                                     tCACCGAGAAGGTGCTGCTGTTCGGACAGGCCATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCTGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCg                                                                                                     >  1:573233/1‑149 (MQ=255)                                       aCCGAGAAGGTGCTGCTGTTCGGACAGGCCATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTa                                                                                                   <  2:573233/149‑1 (MQ=255)                                                    tgctgTTCGGACAGGCCATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGCTTCCCTg                                                                                      >  1:453072/1‑149 (MQ=255)                                                    tgctgTTCGGACAGGCCATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTACGTATCACGATTGCCTg                                                                                      >  1:434889/1‑149 (MQ=255)                                                     gctgTTCGGACAGGCCATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTgg                                                                                     >  1:270270/1‑149 (MQ=255)                                                         ttCGGAGAGGCCATTGAGGCCGCGGGTACAGTCAGAGGCGGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTGCCGTAGCACGATTGCCTGGTTaa                                                                                 <  2:205912/149‑1 (MQ=255)                                                            ggACAGGCCATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCt                                                                              <  2:77863/149‑1 (MQ=255)                                                                  gCCATTCAGACCGCCGGTACAGTAAAAGGCAGTGGTTCCAACGAGCACGCTAAGGCGGACTGGAGGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATACCGATTGCCTGGTTAACCTGGTCGa                                                                        >  2:635607/1‑149 (MQ=255)                                                                   ccATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGac                                                                       >  1:430580/1‑149 (MQ=255)                                                                           aCCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCggg                                                               <  2:400235/149‑1 (MQ=255)                                                                           aCCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCggg                                                               <  1:544892/149‑1 (MQ=255)                                                                                     cAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGAc                                                     <  2:148149/149‑1 (MQ=255)                                                                                       gTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGAcgt                                                   <  2:563713/149‑3 (MQ=255)                                                                                       gTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCCCCATTACTACGTCTGTGATGAAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCCGGGCACGAAGACtt                                                   >  2:172661/1‑149 (MQ=255)                                                                                              gCCGGGGCTACAACCAGCACGATGAGTCGAACGCGAAGGAAATGCAAAAGAAGCGGACGATCTCCAGCACCACAGCGGTGCTGGAGTTCACGGACCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACCTCTCGGaa                                            <  2:18132/149‑1 (MQ=255)                                                                                                  ttgTTCCAACCAGCACGCAAAGTAGGGCCGGAGGGAGATGGAAACGCTCCGGGGTAGCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATa                                        <  2:103200/147‑1 (MQ=255)                                                                                                   ggTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTACGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATAc                                       <  2:434889/149‑1 (MQ=255)                                                                                                                     aaGTCGGACTGGATGGAGATGGAAAAGCAGCATGGGATCTGCATTACTACGTCTGGGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCAGGGGCACGAAGACTTCTCGGCAGATACCTATCGTACCCTGAcggc                     >  2:441643/1‑149 (MQ=255)                                                                                                                        tCGTACTGCATGGAGTTTGAAGAGCAGTGTGGGATCTCTATTACTACGTCTGTGATGCGGCTTACGCATCTCGATGGCCCCGTTAACCCGCTCGACACCCCGGGGCACGCAGACCTCTCGGTAGATACCTATCGTACCCTGACGGCGGt                  <  2:70125/149‑1 (MQ=255)                                                                                                                           ggcgggaggggGATGGAAGAGCTGGGTGGGAAGTGTATTACTAAGTGTGTGATGCAGGCTCCGCATCACGATTGCCTGGTTAACCTGCTCGACGCCCCGGGGCACGAACCCTTCCTGGCAGATACCTATCTGATCCGGACGGCGGTGGa               <  2:635582/145‑1 (MQ=255)                                                                                                                                       atggaAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGa                                                                        >  1:621345/1‑80 (MQ=255)                                                                                                                                       atggaAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGa                                                                        <  2:621345/80‑1 (MQ=255)                                                                                                                                        tggaAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGACGCAGTTTCCGTATCACGATTGCCTGCTTAACCTGCTCGCCACCCCGGCGCACGAAGCCTTCTGGCAAGATACCCATCGTACCCTCCCGCCGGTGGCCTCCTCCatgatg  >  2:236824/1‑149 (MQ=255)                                                                                                                                        |                                                                                                                                                     CCATTATTTCTCACCCGGACGCCGGTAAGACTACCATCACCGAGAAGGTGCTGCTGTTCGGACAGGCCATTCAGACCGCCGGTACAGTAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATG  >  NC_000913_3_pae_tpiA/4609379‑4609662

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