New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4450215 =226 (0.870)4 (0.020) 4/282 NT 1.8% coding (254/957 nt) ytfQ galactofuranose ABC transporter periplasmic binding protein
?NC_000913 4450233 = 222 (0.880)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                                                                                                                                                     <  6:944442/149‑1
TTGCCTTCGAGGATGTTGTCGGCAGTGACGGTGGTCATATAGAGAGATTTGTCTTTCACATCAATGGAACGATCGAGCAAGAATACCGGGATTTCGGCATCTTTCGCCTCTTTTAATACCGGTTCCCAACCTGTCGCGACCACCGGAGC                                                                                                                                                     >  7:768796/1‑149
                                                                                                                                ACCTGTCGCGACCACCGGAGCGACAGGTTGGGAACCGGTATTAAAAGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATC                                                                                    >  5:1715391/1‑86
                                                                                                                                ACCTGTCGCGACCACCGGAGCGACAGGTTGGGAACCGGTATTAAAAGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATC                                                                                    <  6:1715391/86‑1
                                                                                                                                     TCGCGACCACCGGAGCGACAGGTTAGGAACCGGTATTAAAAGAGGCGAAAGA                                                                                                                 <  5:324318/52‑1
                                                                                                                                     TCGCGACCACCGGAGCGACAGGTTAGGAACCGGTATTAAAAGAGGCGAAAGA                                                                                                                 >  6:324318/1‑52
                                                                                                                                                   GCGACAGGTTGGGAACCGGTATTAAAAGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACATCCTCGAAGGCAAGTTGATTGGTGA  >  6:329243/1‑149
                                                                                                                                                     GACAGGTTGGGAACCGGTATTAAAAGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACATCCTCGAAGGCAA              <  7:740115/135‑1
                                                                                                                                                     GACAGGTTGGGAACCGGTATTAAAAGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACATCCTCGAAGGCAA              >  8: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.