New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2330390 =	23 (0.890)	3 (0.130)	3/114	2.0	16.8%	pseudogene (4013/4605 nt)	yfaS	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
	?	NC_000913	2330449 =	10 (0.450)					pseudogene (3954/4605 nt)	yfaS	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Rejected: Frequency below/above cutoff threshold.

CTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCG  >  NC_000913/2330382‑2330515                                                                    |                                                                   cTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGc                                                                    <  2:383706/68‑1 (MQ=255)         cctccagttgttttgcctgcgcctgcgggcgcatgggctaaacataagttaactggaggCATTGTCGc                                                            <  1:130651‑M2/9‑1 (MQ=255)         CCTCCAGTTAACTTATGTCTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGc                                                            >  2:23868‑M1/1‑68 (MQ=255)              AGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCatat                                                       <  2:271772‑M1/68‑1 (MQ=255)                 TAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCatat                                                       <  2:885550‑M1/65‑1 (MQ=255)                 TAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCatat                                                       >  1:885546‑M1/1‑65 (MQ=255)                        cctgcgcctgcgggcgcatgggctaaacataagttaactggaggCATTGTCGCCATATATTTCgccag                                             <  2:705523‑M2/24‑1 (MQ=255)                              gCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGtt                                       >  2:596964/1‑68 (MQ=255)                              cctgcgggcgcatgggctaaacataagttaactggaggCATTGTCGCCATATATTTCGCCAGCCAGtt                                       <  1:311613‑M2/30‑1 (MQ=255)                               cccATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTc                                      >  1:846382/1‑68 (MQ=255)                               cccATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTc                                      >  2:654047/1‑68 (MQ=255)                                ccATGCGCCCGCAGGCGCAGGCAAAACAACAGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCa                                     >  1:395183/1‑68 (MQ=255)                                              gcgcaggcAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAgcgcgcgc                       <  2:559713/68‑1 (MQ=255)                                                 caggcAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCATTTCATGGCCAGCGCGCGCTct                    <  2:130651/68‑2 (MQ=255)                                                   ggcAATACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAAt                  <  1:782198/68‑1 (MQ=255)                                                   ggcAATACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAAt                  <  1:938480/68‑1 (MQ=255)                                                   ggcAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAAt                  <  1:321904/68‑1 (MQ=255)                                                            actggaggCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTg         >  2:321807‑M2/9‑68 (MQ=255)                                                                 aggCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAg    >  2:641485‑M2/4‑68 (MQ=255)                                                                   gCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCg  >  1:459495‑M2/2‑68 (MQ=255)                                                                    |                                                                   CTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCG  >  NC_000913/2330382‑2330515

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.