New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3858205	26 (0.680)	4 (0.120)	4/240	4.3	13.0%	coding (200/1494 nt)	yidJ	sulfatase/phosphatase superfamily protein
	?	NC_000913	= 3858228	30 (0.870)					coding (177/1494 nt)	yidJ	sulfatase/phosphatase superfamily protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGAC  >  NC_000913/3858074‑3858350                                                                                                                                    |                                                                                                                                                  tGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCg                                                                                                                                  <  1:84934/149‑1 (MQ=255)      aaaTAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTaaa                                                                                                                             >  6:109857/1‑149 (MQ=255)                    gTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCgg                                                                                                                                 <  8:109788/131‑1 (MQ=255)                    gTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCgg                                                                                                                                 >  7:109788/1‑131 (MQ=255)                            gCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTa                                                                                                       >  8:4750/1‑149 (MQ=255)                                  tGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAAc                                                                                                  <  7:4750/148‑1 (MQ=255)                                    ccGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCg                                                                                                <  5:109857/148‑1 (MQ=255)                                         gTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCg                                                                                                >  6:58907/1‑143 (MQ=255)                                         gTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCg                                                                                                <  5:58907/143‑1 (MQ=255)                                                                 cTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGaa                                                                  >  1:48355/1‑149 (MQ=255)                                                                        tAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGgtgt                                                                                    <  6:37257/124‑1 (MQ=255)                                                                        tAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGgtgt                                                                                    >  5:37257/1‑124 (MQ=255)                                                                                       gagaTGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTc                                                           >  6:127314/1‑134 (MQ=255)                                                                                       gagaTGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTc                                                           <  5:127314/134‑1 (MQ=255)                                                                                       gagaTGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTAt                                             >  5:139244/1‑148 (MQ=255)                                                                                         gaTGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCaa                                          >  6:178205/1‑149 (MQ=255)                                                                                             tttttGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGaa                                                                  <  6:98075/121‑1 (MQ=255)                                                                                                    cTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTAt                                             <  7:9727/135‑1 (MQ=255)                                                                                                    cTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTAt                                             >  8:9727/1‑135 (MQ=255)                                                                                                      ggCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTAtt                             <  2:48355/149‑1 (MQ=255)                                                                                                         gcgACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAg                          >  3:68392/1‑149 (MQ=255)                                                                                                           gACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCg                        <  6:139244/149‑1 (MQ=255)                                                                                                                  ttgGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCg                        <  4:114202/142‑1 (MQ=255)                                                                                                                  ttgGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCg                        >  3:114202/1‑142 (MQ=255)                                                                                                                                 gACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGAc  >  2:10563/1‑149 (MQ=255)                                                                                                                                 gACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGAc  >  8:65994/1‑149 (MQ=255)                                                                                                                                    |                                                                                                                                                  TGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGAC  >  NC_000913/3858074‑3858350

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