New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2333629	15 (0.490)	4 (0.150)	4/112	2.1	18.7%	pseudogene (774/4605 nt)	yfaS	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
	?	NC_000913	= 2334464	22 (0.850)					coding (496/624 nt)	yfaT	DUF1175 family protein, putative host defense protein
Rejected: Frequency below/above cutoff threshold.

TCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGG  >  NC_000913/2334402‑2334523                                                               |                                                            tCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTcgtgg                                                        <  1:1153544‑M2/68‑6 (MQ=255) tCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTcgtgg                                                        <  1:1152030‑M2/68‑6 (MQ=255)     gtgtCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTcgtggcgct                                                    >  1:439592‑M2/1‑59 (MQ=255)      tgtCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTcgtggcgctt                                                   <  1:658663‑M2/68‑11 (MQ=255)           ttCCATGTCATTTGTTTTTGCTGACTGTCTGCGCTAATTCCGTTGTCAGTTTTCGGGGCGCTTCCGgt                                              <  1:707065/68‑1 (MQ=255)             ccATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGgtgt                                            <  2:332054/68‑1 (MQ=255)             ccATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGgcgt                                            <  2:334016/68‑1 (MQ=255)               aTGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTgg                                          <  2:1152025/68‑1 (MQ=255)               aTGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTgg                                          <  2:1153539/68‑1 (MQ=255)                  tCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTTGTAg                                       <  1:946727/68‑1 (MQ=255)                     tAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATg                                    <  2:21061/68‑1 (MQ=255)                         ttgttgCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGt                                >  1:405819/1‑68 (MQ=255)                             cgcAGACTGACTGCGCGCATTCCGTTGTCAGTTTTcgtggcgcttagcaaagtgtcaggcaaagagct                            >  2:392541‑M2/2‑35 (MQ=255)                              gCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGc                           <  2:88786/68‑1 (MQ=255)                              gCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGc                           <  2:730559/68‑1 (MQ=255)                                    tgactgCGCGCATGCCGTTGGCAGTTTTcgtggcgcttagcaaagtgtcaggcaaagagcttctggtg                     >  1:549548‑M2/1‑28 (MQ=255)                                     gactgCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATcc                    <  1:686306/68‑1 (MQ=255)                                        tgCGCGCATTCCGTTGTCAGTTTTcgtggcgcttagcaaagtgtcaggcaaagagcttctgg                       >  1:637356‑M2/1‑24 (MQ=255)                                        tgCGCGCATTCCGTTGTCAGTTTTcgtggcgcttagcaaagtgtcaggcaaagagcttctgg                       <  2:637355‑M2/62‑39 (MQ=255)                                         gcgcgcATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGGAGATGACGTAACGCCCCATCCAGAc                <  2:844713/68‑1 (MQ=255)                                          cgcgcATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGAcc               <  2:228341/68‑1 (MQ=255)                                          cgcgcATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGAcc               <  1:620325/68‑1 (MQ=255)                                          cgcgcATTCCGTTGTCAGTTTGCGCGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGAcc               <  2:816390/68‑1 (MQ=255)                                            cgcATTCCGTTGTCAGTTTTCGCGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCAt             <  2:148322/68‑1 (MQ=255)                                            cgcATTCCGTTGTCAGTTTTCGCGGCGCTTCCGGTGCGGTAGATGACGTAACGCCCCATCCAGACCAt             <  2:148323/68‑1 (MQ=255)                                            cgcATTCCGTTGTCAGATTTCGTGGCGCTGCCGGGGTGGTTGAAGACGTAACGCCCCATCCAGACCAt             <  2:1028396/68‑1 (MQ=255)                                            cgcAGTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCAt             <  2:842464/68‑1 (MQ=255)                                              cATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTa           <  2:391176/68‑1 (MQ=255)                                                 tCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGt        <  2:490790/68‑1 (MQ=255)                                                       gTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTggg  <  1:907305/68‑1 (MQ=255)                                                               |                                                            TCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGG  >  NC_000913/2334402‑2334523

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