New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 374495363 (0.830)3 (0.040) 3/250 NT 4.4% coding (150/1278 nt) yiaN 2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
?NC_000913 = 3744979 72 (1.010)coding (176/1278 nt) yiaN 2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Rejected: Frequency below/above cutoff threshold.

TCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3744802‑3744953
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gaacggaATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACA  <  NC_000913/3744979‑3744832
                                                                                                                                                                                                                                                                                                            
TCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAAC                                                                                                                                                         >  7:200926/1‑149
TCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAAC                                                                                                                                                         >  8:136968/1‑149
                                        tcTTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGTCACAAACGCTGGTGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCA                                                                                                                 <  7:239851/147‑1
                                                    TTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCAC                                                                                                                   <  7:68064/135‑1
                                                    TTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCAC                                                                                                                   >  8:68064/1‑135
                                                                                                                                                     GGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAA    >  3:325765/1‑149
                                                                                                                                                      GAATCGCCCGCAGGGAGAAGCTATCTGCTCCGTTCACCAGCTTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACTCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAAC   <  8:347691/149‑1
                                                                                                                                                       AATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACA  <  1:310717/149‑1
                                                                                                                                                       AATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACC          <  3:23627/141‑1
                                                                                                                                                       AATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACC          >  4:23627/1‑141
                                                                                                                                                                                                                                                                                                            
TCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3744802‑3744953
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gaacggaATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACA  <  NC_000913/3744979‑3744832

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