New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 3746905 =60 (1.350)4 (0.090) 3/272 6.5 7.0% coding (812/987 nt) yiaO 2,3‑diketo‑L‑gulonate‑binding periplasmic protein
?NC_000913_3_bme_pgi 3746922 = 50 (1.180)coding (829/987 nt) yiaO 2,3‑diketo‑L‑gulonate‑binding periplasmic protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTCAGCAGATCAGCTCCCTGCGGCACATCTTTAACAAACATGTCGCGAACCTGATTCCCCAGTGCGTCGCTAAAGGCTTTGCGGTCGAGATCGGTGATGACTTCCACGCCCGCTTCTTTCATGCCGTCGATGATTTTTTGCTGATCTTCAGCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_bme_pgi/3747058‑3746905
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcagcaaAAAATCATCGACGGCATGAAAGAAGCGGGCGTGGAAGTCATCACCGATCTCGACCGCAAAGCCTTTAGCGACGCACTGGGGAATCAGGTTCGCGACATGTTTGTTAAAGATGTGCCGCAGGGAGCTGATCTGCTGAAAGCCG  >  NC_000913_3_bme_pgi/3746922‑3747063
                                                                                                                                                                                                                                                                                                        
TTCAGCAGATCAGCGCCCTGCGGCACATCTTTAACAAACATGTCGCGAACTTGATTACCCATTGCGTCTCTAAAGGCTTTGCGGTCGAGATCGGTGATGACTTCCACGCCCGCTTCTTTCATGCCGTCGATGATTTTTTGCTGATCTTC                                                                                                                                                     <  1:757751/149‑1
     CAGATCAGCTCCCTGCGGCACATCTTTAACAAACATGTCGCGAACCTGATTCCCCAGTGCGACGCTAAAGGCTTTGCGGTCGAGATCGGTGATGACTTCCACGCCCGCTTCTTTCATGCCGTCGATGATTTTTTGCTGATCTTCAGCAA                                                                                                                                                >  2:673848/1‑149
                                            GCGAACCTGATTCCCCAGTGCGTCGCTAAAGGCTTTGCGGTCGAGATCGGTGATGACTTCCACGCCCGCTTCTTTCATGCCGTCGATGATTTTTTGCTGATCTTCAGCA                                                                                                                                                 >  1:283651/1‑109
                                            GCGAACCTGATTCCCCAGTGCGTCGCTAAAGGCTTTGCGGTCGAGATCGGTGATGACTTCCACGCCCGCTTCTTTCATGCCGTCGATGATTTTTTGCTGATCTTCAGCA                                                                                                                                                 <  2:283651/109‑1
                                                                                                                              TCGATGATTTTTTGCTGATCTTCAGCAAAAAATCATCGACGGCATGAAAGAAGCGGGCGTGGAAGTCATCACCGATCTCGAC                                                                                          <  1:398142/82‑1
                                                                                                                              TCGATGATTTTTTGCTGATCTTCAGCAAAAAATCATCGACGGCATGAAAGAAGCGGGCGTGGAAGTCATCACCGATCTCGAC                                                                                          >  2:398142/1‑82
                                                                                                                                     tgTTTTGCTGATCTTCAGCAAAAAATCATCGACGGCATGAAAGAAGCGGGCGTGGAAGTCATCACCGATCTCGACCGCAAAGCCTTTAGCGACGCACTGGGGAATCAGGTTCGCGACATGTTTGTTAAAGATGTGCCGCAGGGAGCTGA                >  1:543381/3‑149
                                                                                                                                        TTTGCTGATCTTCAGCAAAAAATCATCGACGGCATGAAAGAAGCGGGCGTGGAAGTCATCACCGATCTCGACCGCAAAGCCTTTAGCGACGCACTGGGGAATCAGGTTCGCGACATGTTTGTTAAAGATGTGCCGCAGGGAGCTGATCT             <  2:543381/149‑1
                                                                                                                                                   TCAGAAAAAAATCACCGACGGCATGAAATAAGCGGGCGTGGAAGTCATCACCGATCTCGACCGCAAAGCCTTTATCGACGCACTGGGGAATCAGGTTCGCGACATGTTTGTTAAAGATGTGCCGCAGGGAGCTGATCTGCTGAAAGCCG  <  1:437685/149‑1
                                                                                                                                                                                                                                                                                                        
TTCAGCAGATCAGCTCCCTGCGGCACATCTTTAACAAACATGTCGCGAACCTGATTCCCCAGTGCGTCGCTAAAGGCTTTGCGGTCGAGATCGGTGATGACTTCCACGCCCGCTTCTTTCATGCCGTCGATGATTTTTTGCTGATCTTCAGCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_bme_pgi/3747058‑3746905
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcagcaaAAAATCATCGACGGCATGAAAGAAGCGGGCGTGGAAGTCATCACCGATCTCGACCGCAAAGCCTTTAGCGACGCACTGGGGAATCAGGTTCGCGACATGTTTGTTAAAGATGTGCCGCAGGGAGCTGATCTGCTGAAAGCCG  >  NC_000913_3_bme_pgi/3746922‑3747063

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