New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	3982818 =	28 (0.950)	4 (0.140)	4/266	3.5	12.4%	intergenic (+7/‑140)	proM/aslB	tRNA‑Pro/putative AslA‑specific sulfatase‑maturating enzyme
	?	NC_000913_3_pae_pgi	3982838 =	30 (1.080)					intergenic (+27/‑120)	proM/aslB	tRNA‑Pro/putative AslA‑specific sulfatase‑maturating enzyme
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTG  >  NC_000913_3_pae_pgi/3982690‑3982982                                                                                                                                                     |                                                                                                                                                 ccacGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGctt                                                                                                                                                  >  1:334487/1‑149 (MQ=255)    cgccattAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGCAGCGCAGCTTGGTAGCGCAACTGGGTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGtttttt                                                                                                                                               <  2:237105/144‑1 (MQ=255)       cTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGtt                                                                                                                                            <  1:58059/149‑1 (MQ=255)        tttAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGCCGGGGTTTTTTGttt                                                                                                                                           >  2:163778/1‑149 (MQ=255)                aTTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCa                                                                                                                                   >  2:50000/1‑149 (MQ=255)                    aaCTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGACCCCCGCTTCGACGGGGTGTTTTGTTTTCTGTGCATTTc                                                                                                                               >  1:355726/1‑149 (MQ=255)                    aaCTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACACCGCTTCTGCGGGGTTTTTTGTTTTCTGTGCATTTc                                                                                                                               >  1:355683/1‑149 (MQ=255)                               cAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCAccctccc                                                                                                                    <  2:98452/149‑1 (MQ=255)                                aaaaaGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACGCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCAccctccct                                                                                                                   >  1:42071/1‑149 (MQ=255)                                  aaaGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGCCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTc                                                                                                                 >  2:510665/1‑149 (MQ=255)                                      cAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGCACTTTCGGGAACCCAGATGGTGCGTGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCaa                                                                                                             <  1:123225/149‑1 (MQ=255)                                               ggCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGGCGGAGGTTCGAATCGTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGccc                                                                                                    >  2:58809/1‑149 (MQ=255)                                               ggCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTAGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGccc                                                                                                    <  2:310484/149‑1 (MQ=255)                                                gCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATaaa                                                                                                         <  1:284698/143‑1 (MQ=255)                                                gCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATaaa                                                                                                         >  2:284698/1‑143 (MQ=255)                                                              ttGGTAGCGCAACTGGTTGGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCTTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGc                                                                                     <  2:421651/149‑1 (MQ=255)                                                               tGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGcc                                                                                    <  1:50000/149‑1 (MQ=255)                                                                 gTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGcccc                                                                                  >  2:199957/1‑149 (MQ=255)                                                                 gTAGCGCAACTGGTTTGGGACCAGTGCGCCGGCACTGCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAAGTCATTGcccc                                                                                  >  1:126006/1‑149 (MQ=255)                                                                                      cAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTATCTCTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTaa                                                             >  2:222876/1‑149 (MQ=255)                                                                                         tGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGAACAAGCTGTGATGTGCTTAAATCCTCGCTCCCTTCGCAATAAACGGCCTTAATAACACATTGCCCCACGGTATGATTTCGCCCTTAACGt                                                          >  1:253322/1‑149 (MQ=255)                                                                                           gctccgtggTTCGAATCCTTCCTCTCCGTCCAATTTTTAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCTTCACCCTCCCTTCGCAATAAACGCCCTTAATAACTCATTGCCCCACGATATGATTTCGCCCTTAACGTAt                                                        <  2:334487/142‑1 (MQ=255)                                                                                             tCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTg                                                      >  2:153607/1‑149 (MQ=255)                                                                                               ggAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGaa                                                    <  2:284694/149‑1 (MQ=255)                                                                                                           ccTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAg                                        <  1:264729/149‑1 (MQ=255)                                                                                                               tctcGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCaa                                    <  2:253322/149‑1 (MQ=255)                                                                                                                 tcGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAgg                                  >  1:344569/1‑149 (MQ=255)                                                                                                                         cAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTCAACGTATTGCAGGATGACTTCAGGCACGAAGCGCCCa                          >  2:120083/1‑149 (MQ=255)                                                                                                                                    cgcGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAg               <  2:175705/147‑1 (MQ=255)                                                                                                                                      cGCTTCGGCGGGGTTTTTTGGTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGt             <  1:163778/149‑1 (MQ=255)                                                                                                                                                 ggTTGTTTGTTTTCTGTGCATTTCGTCACCCGCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTg  <  2:513870/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTG  >  NC_000913_3_pae_pgi/3982690‑3982982

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