New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4450215 =150 (0.820)4 (0.020) 4/248 NT 2.7% coding (254/957 nt) ytfQ galactofuranose ABC transporter periplasmic binding protein
?NC_000913 4450233 = 144 (0.810)coding (272/957 nt) ytfQ galactofuranose ABC transporter periplasmic binding protein
Rejected: Frequency below/above cutoff threshold.

TTGCCTTCGAGGATGTTGTCGGCAGTGACGGTGGTCATATAGAGAGATTTGTCTTTCACATCAATGGAACGATCGAGCAAGAATACCGGGATTTCGGCATCTTTCGCCTCTTTTAATACCGGTTCCCAACCTGTCGCGACCACCGGAGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4450365‑4450215
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgaCAGGTTGGGAACCGGTATTAAAAGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACATCCTCGAAGGCAAGTTGATTGGTGA  >  NC_000913/4450233‑4450377
                                                                                                                                                                                                                                                                                                        
TTGCCTTCGAGGATGTTGTCGGCAGTGACGGTGGTCATATAGAGAGATTTGTCTTTCACATCAATGGAACGATCGAGCAAGAATACCGGGATTTCGGCATCTTTCGCCTCTTTTAATACCGGTTCCCAACCTGTCGCGACCACCGGAGC                                                                                                                                                     <  2:944442/149‑1
TTGCCTTCGAGGATGTTGTCGGCAGTGACGGTGGTCATATAGAGAGATTTGTCTTTCACATCAATGGAACGATCGAGCAAGAATACCGGGATTTCGGCATCTTTCGCCTCTTTTAATACCGGTTCCCAACCTGTCGCGACCACCGGAGC                                                                                                                                                     >  3:768796/1‑149
                                                                                                                                ACCTGTCGCGACCACCGGAGCGACAGGTTGGGAACCGGTATTAAAAGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATC                                                                                    >  1:1715391/1‑86
                                                                                                                                ACCTGTCGCGACCACCGGAGCGACAGGTTGGGAACCGGTATTAAAAGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATC                                                                                    <  2:1715391/86‑1
                                                                                                                                     TCGCGACCACCGGAGCGACAGGTTAGGAACCGGTATTAAAAGAGGCGAAAGA                                                                                                                 <  1:324318/52‑1
                                                                                                                                     TCGCGACCACCGGAGCGACAGGTTAGGAACCGGTATTAAAAGAGGCGAAAGA                                                                                                                 >  2:324318/1‑52
                                                                                                                                                   GCGACAGGTTGGGAACCGGTATTAAAAGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACATCCTCGAAGGCAAGTTGATTGGTGA  >  2:329243/1‑149
                                                                                                                                                     GACAGGTTGGGAACCGGTATTAAAAGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACATCCTCGAAGGCAA              <  3:740115/135‑1
                                                                                                                                                     GACAGGTTGGGAACCGGTATTAAAAGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACATCCTCGAAGGCAA              >  4:740115/1‑135
                                                                                                                                                                                                                                                                                                        
TTGCCTTCGAGGATGTTGTCGGCAGTGACGGTGGTCATATAGAGAGATTTGTCTTTCACATCAATGGAACGATCGAGCAAGAATACCGGGATTTCGGCATCTTTCGCCTCTTTTAATACCGGTTCCCAACCTGTCGCGACCACCGGAGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4450365‑4450215
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgaCAGGTTGGGAACCGGTATTAAAAGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACATCCTCGAAGGCAAGTTGATTGGTGA  >  NC_000913/4450233‑4450377

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