New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2331370	33 (1.120)	4 (0.140)	4/246	2.9	10.9%	pseudogene (3033/4605 nt)	yfaS	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
	?	NC_000913	= 2331375	34 (1.220)					pseudogene (3028/4605 nt)	yfaS	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Rejected: Frequency below/above cutoff threshold.

GTTTCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTA  >  NC_000913/2331254‑2331523                                                                                                                          |                                                                                                                                                     gTTTCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAACTGTCttttgtggataacagctggcccgttg                                                                                                                            <  2:85441‑M2/148‑27 (MQ=255)  tttCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttg                                                                                                                                                 >  3:124095‑M2/1‑121 (MQ=255)  tttCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttg                                                                                                                                                 <  4:124095‑M2/126‑6 (MQ=255)  tttCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGACCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGc                                                                                                                          >  4:79581/1‑149 (MQ=255)        tactTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCAttt                                                                                                                    >  5:65435/1‑149 (MQ=255)                aGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAACTGTCttttgtggataacagctggcccgttgaacagcagaaaaatgt                                                                                                            >  1:85441‑M2/1‑107 (MQ=255)                    ggATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCt                                                                                                        >  5:17511/1‑149 (MQ=255)                     gATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttgtg                                                                                                                                               <  4:28313‑M2/109‑8 (MQ=255)                     gATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttgtg                                                                                                                                               >  3:28313‑M2/1‑102 (MQ=255)                         ttGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCAc                                                                                                    <  6:17511/148‑1 (MQ=255)                         ttGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGACCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCAc                                                                                                    <  3:79581/148‑1 (MQ=255)                           gCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACtt                                                                                                  >  8:9486/1‑148 (MQ=255)                           gCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACtt                                                                                                  >  5:22276/1‑148 (MQ=255)                             tCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTaa                                                                                                >  5:131783/1‑148 (MQ=255)                                 cTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCAc                                                                                                                              >  1:163139/1‑114 (MQ=255)                                 cTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCAc                                                                                                                              <  2:163139/114‑1 (MQ=255)                                          cAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCt                                                                                  <  7:9486/149‑1 (MQ=255)                                            aCATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTg                                                                                 >  2:68486/1‑148 (MQ=255)                                               tCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttgt                                                                                                                                                >  4:12188‑M2/1‑76 (MQ=255)                                               tCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttgt                                                                                                                                                <  3:12188‑M2/82‑7 (MQ=255)                                               tCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTttgttg                                                                                                               <  6:16522/115‑1 (MQ=255)                                               tCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTttgttg                                                                                                               >  5:16522/1‑115 (MQ=255)                                                                         aTGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGc                                                   <  6:65435/149‑1 (MQ=255)                                                                              aTTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCtt                                              >  6:68935/1‑149 (MQ=255)                                                                               tttttCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCttt                                             <  5:68935/149‑1 (MQ=255)                                                                                ttttCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTtg                                            >  5:44935/1‑149 (MQ=255)                                                                                 tttCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTtg                                            <  6:22276/148‑1 (MQ=255)                                                                                  ttCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTtgt                                           <  1:68486/148‑1 (MQ=255)                                                                                       ctgTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGcagc                                      <  6:118617/148‑1 (MQ=255)                                                                                          ttCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGcagcgtca                                  <  6:44935/149‑1 (MQ=255)                                                                                          ttCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGcagcgt                                    <  6:131783/147‑1 (MQ=255)                                                                                             aaCGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGcagcgtcagcg                               >  5:54013/1‑149 (MQ=255)                                                                                             aaCGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGcagcgtcagc                                <  2:35004/148‑1 (MQ=255)                                                                                             aaCGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGcagcgtcagc                                >  1:35004/1‑148 (MQ=255)                                                                                                      gCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGaaa                                                             <  4:124207/110‑1 (MQ=255)                                                                                                        tGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCAtc                    >  5:13266/1‑149 (MQ=255)                                                                                                         gTTATCCACAAAAGACAGttttgtggataacagctggcccgttgaacagcagaaaaatgtcatgctcggtg                                                                                                 <  8:116965‑M2/71‑54 (MQ=255)                                                                                                         gTTATCCACAAAAGACAGttttgtggataacagctggcccgttgaacagcagaaaaatgtcatgctcggtg                                                                                                 >  7:116965‑M2/1‑18 (MQ=255)                                                                                                           tATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCAtct                   >  1:113/1‑147 (MQ=255)                                                                                                              ccACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAAtt                                                      >  6:51713/1‑109 (MQ=255)                                                                                                              ccACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAAtt                                                      <  5:51713/109‑1 (MQ=255)                                                                                                                 caAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCt                            >  2:31991/1‑132 (MQ=255)                                                                                                                  aaaaGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTg          >  3:128880/1‑149 (MQ=255)                                                                                                                          gTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTa  <  6:13266/149‑1 (MQ=255)                                                                                                                          |                                                                                                                                                     GTTTCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTA  >  NC_000913/2331254‑2331523

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