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.

TTTCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTA  >  NC_000913/2331255‑2331523                                                                                                                         |                                                                                                                                                     tttCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttg                                                                                                                                                 <  4:124095‑M2/126‑6 (MQ=255) tttCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttg                                                                                                                                                 >  3:124095‑M2/1‑121 (MQ=255) tttCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGACCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGc                                                                                                                          >  4:79581/1‑149 (MQ=255)       tactTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCAttt                                                                                                                    >  5:65435/1‑149 (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                                                                                                                                                <  3:12188‑M2/82‑7 (MQ=255)                                              tCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttgt                                                                                                                                                >  4:12188‑M2/1‑76 (MQ=255)                                              tCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTttgttg                                                                                                               <  6:16522/115‑1 (MQ=255)                                              tCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTttgttg                                                                                                               >  5:16522/1‑115 (MQ=255)                                                                        aTGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGc                                                   <  6:65435/149‑1 (MQ=255)                                                                           aCATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttgtggataacagctggcccgttgaacagcagaaaaatgtcatgctcggtggtggcgataacgcgctgatgttgcccgagcaggcgagcaatatccggct                                                <  1:85441‑M2/149‑103 (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)                                                                                           cAACGGGCCAGCTGTTATCCACAAAAGACAGttttgtggataacagctggcccgttgaacagcagaaaaatgtcatgctcggtggtggcgataacgcgctgatgttgcccgagcaggcgagcaatatccggctacaaagtagtgaaac                                 >  2:85441‑M2/1‑31 (MQ=255)                                                                                            aaCGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGcagcgtcagcg                               >  5:54013/1‑149 (MQ=255)                                                                                            aaCGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGcagcgtcagc                                >  1:35004/1‑148 (MQ=255)                                                                                            aaCGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGcagcgtcagc                                <  2:35004/148‑1 (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                                                      <  5:51713/109‑1 (MQ=255)                                                                                                             ccACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAAtt                                                      >  6:51713/1‑109 (MQ=255)                                                                                                                caAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCt                            >  2:31991/1‑132 (MQ=255)                                                                                                                 aaaaGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTg          >  3:128880/1‑149 (MQ=255)                                                                                                                         gTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTa  <  6:13266/149‑1 (MQ=255)                                                                                                                         |                                                                                                                                                     TTTCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTA  >  NC_000913/2331255‑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.