New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3698459 =93 (0.990)4 (0.040) 4/242 NT 4.2% coding (246/1680 nt) bcsG DUF3260 family cellulose production inner membrane protein
?NC_000913 3698508 = 92 (1.030)coding (295/1680 nt) bcsG DUF3260 family cellulose production inner membrane protein
Rejected: Frequency below/above cutoff threshold.

ACAGGAATAACCAGGCCACTAATAAAACAAAAATGGCCCCAATCATCTGCCAGTTAATAAAGCGTGTGACAAGGTCGATTAAATAATCGGTACTGAACCCCGCCACCTGCGAACCCTGGCTCATTATGCTTTCCGGGCCAGGCAACCAGGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3698610‑3698459
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caggtgGCGGGGTTCAGTACCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTATTCCTGTCACAATGGATTCGCATTA  >  NC_000913/3698508‑3698628
                                                                                                                                                                                                                                                                                 
ACAGGAATAACCAGGCCACTAATAAAACAAAAATGGCCCCAATCATCTGCCAGTTAATAAAGCGTGTGACAAGGTCGATTAAATAATCGGTACTGAACCCCGCCCCCTGCCACCCCTGGCTCATTATGCTTGCCGGGCCCGGCACCCAG                                                                                                                              >  2:309514/1‑149
ACAGGAATAACCAGGCCACTAATAAAACAAAAATGGCCCCAATCATCTGCCAGTTAATAAAGCGCGTGACAAGGTCGATTAAATAATCGGTACTGAACCCCGCCACCTGCGAACCCTGGCTCATTATGCTTTCCGGGCCAGGCAACCAG                                                                                                                              <  2:146558/149‑1
 CAGGAATAACCAGGCCACTAATAAAACAAAAATGGCCCCAATCATCTGCCAGTTAATAAAGCGTGTGACAAGGTCGATTAAATAATCGGTACTGAACCCCGCCACCTGCGAACCCTGGCTCATTATGCTTTCCGGGCCAGGCAACCAGG                                                                                                                             >  1:149685/1‑149
  AGGAATAACCAGGCCACTAATAAAACAAAAATGGCCCCAATCATCTGCCAGTTAATAAAGCGTGTGACAAGGTCGATTAAATAATCGGTACTGAACCCCGCCACCTGCGAACCCTGGCTCATTATGCTTTCCGGGCCAGGCAACCAGGT                                                                                                                            <  8:223067/149‑1
  AGGAATAACCAGGCCACTAATAAAACAAAAATGGCCCCAATCATCTGCCAGTTAATAAAGCGTGTGACAAGGTCGATTAAATAATCGGTACTGAACCCAGCCACCTGCGAACCCTGGCTCATTATGCTTTCCGGGCCAGGCAACCAGGT                                                                                                                            <  3:40241/149‑1
          CCAGGCCACTAATAAAACAAAAATGGCCCCAATCATCTGCCAGTTAATAAAGCGTGTGACAAGGTCGATTAAATAATCGGTACTGAACCCCGCCACCTGCGATCCCTGGCTCATTATGCTTTCCGGGCCAGGCAACC                                                                                                                                <  2:364653/137‑1
          CCAGGCCACTAATAAAACAAAAATGGCCCCAATCATCTGCCAGTTAATAAAGCGTGTGACAAGGTCGATTAAATAATCGGTACTGAACCCCGCCACCTGCGACCCCTGGCTCATTATGCTTTCCGGGCCAGGCAACC                                                                                                                                >  1:364653/1‑137
          CCAGGCCACTAATAAAACAAAAATGGCCCCAATCATCTGCCAGTTAATAAAGCGTGTGACAAGGTCGATTAAATAATCGGTACTGAACCCCGCCACCTGCGAACCCTGGCTCATTATGCTTTCCGGGCCAGGCAACCAG                                                                                                                              >  6:426731/1‑139
          CCAGGCCACTAATAAAACAAAAATGGCCCCAATCATCTGCCAGTTAATAAAGCGTGTGACAAGGTCGATTAAATAATCGGTACTGAACCCCGCCACCTGCGAACCCTGGCTCATTATGCTTTCCGGGCCAGGCAACCAG                                                                                                                              <  5:426731/139‑1
                                                                        GGTCGATTAAATAATCGGTACTGAACCCCGCCACCTGCGAACCCTGGCTCATTATGCTTTCCGGGCCAGGCAACCAGGT                                                                                                                            <  1:464983/79‑1
                                                                        GGTCGATTAAATAATCGGTACTGAACCCCGCCACCTGCGAACCCTGGCTCATTATGCTTTCCGGGCCAGGCAACCAGGT                                                                                                                            >  2:464983/1‑79
                                                                                             TGAACCCCGCCACCTGCGAACCCTGGCTCATTATGCTTTCCGGGCCAGGCAACCAGGTGGCGGGGTTCAGTACCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCC                                 >  2:365310/1‑149
                                                                                                  CCCGCCACCTGCGAACCCTGGCTCATTATGCTTTCCGGGCCAGGCAACCAGGTGGCGGGGTTCAGTACCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGT                             >  8:474441/1‑148
                                                                                                      CCACCTGCGAACCCTGGCTCATTATGCTTTCCGGGCCAGGCAACCAGGTGGCGGGGTTCAGTACCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTATTC                        <  1:365310/149‑1
                                                                                                                             ATGCTTTCCGGGCCAGGCAACCAGGTGGCGGGGTTCAGTACCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTATTCCTGTCACAATGGATTCGCATTA  <  7:474441/148‑1
                                                                                                                                                  CAGGTGGCGGGGTTCAGTACCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTATTCCTGTCACAAT              >  2:312082/1‑115
                                                                                                                                                  CAGGTGGCGGGGTTCAGTACCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTATTCCTGTCACAAT              <  1:312082/115‑1
                                                                                                                                                                                                                                                                                 
ACAGGAATAACCAGGCCACTAATAAAACAAAAATGGCCCCAATCATCTGCCAGTTAATAAAGCGTGTGACAAGGTCGATTAAATAATCGGTACTGAACCCCGCCACCTGCGAACCCTGGCTCATTATGCTTTCCGGGCCAGGCAACCAGGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3698610‑3698459
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caggtgGCGGGGTTCAGTACCGATTATTTAATCGACCTTGTCACACGCTTTATTAACTGGCAGATGATTGGGGCCATTTTTGTTTTATTAGTGGCCTGGTTATTCCTGTCACAATGGATTCGCATTA  >  NC_000913/3698508‑3698628

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