New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3973849 =125 (1.110)3 (0.030) 3/298 12.8 2.4% coding (242/882 nt) rffH glucose‑1‑phosphate thymidylyltransferase
?NC_000913 3973867 = 126 (1.150)coding (260/882 nt) rffH glucose‑1‑phosphate thymidylyltransferase
Rejected: Position hash score below cutoff.

CTTCGGTGCGCGCCGCAACATGACGCAGCTTCGGACTGAAGCCCTGACCGAAGAAGATGTTATCGCCCAGCACCAGACAAGAAGGTTCACCATTAAGGAAGGTTTCACCGATGATAAAGGCCTGCGCCAGACCGTCCGGGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3973992‑3973849
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctggCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCG  >  NC_000913/3973867‑3974065
                                                                                                                                                                                                                                                                                                                                                       
CTTCGGTGCGCGCCGCAACATGACGCAGCTTCGGACTGAAGCCCTGACCGAAGAAGATGTTATCGCCCAGCACCAGACAAGAAGGTTCACCATTAAGGAAGGTTTCACCGATGATAAAGGCCTGCGCCAGACCGTCCGGGC                                                                                                                                                                                                            <  1:523290/141‑1
CTTCGGTGCGCGCCGCAACATGACGCAGCTTCGGACTGAAGCCCTGACCGAAGAAGATGTTATCGCCCAGCACCAGACAAGAAGGTTCACCATTAAGGAAGGTTTCACCGATGATAAAGGCCTGCGCCAGACCGTCCGGGC                                                                                                                                                                                                            <  1:580700/141‑1
         GCGCAGCAACATGACGCAGCTTCGGACGGAAGCCCTGACCGAAGAAGATGTTATCGCCCAGCACCAGACAAGAAGGTTCACCATTAAGGAAGGTTTCACCGATGATAAAGGCCTGCGCCAGACCGTCCGGGC                                                                                                                                                                                                            <  3:589387/132‑1
          CGCCGCAACATGACGCAGCTTCGGACTGAAGCCCTGACCGAAGAAGAGGTTATCGCCCAGCACCAGACAAGAAGGTTCACCATTAAGGAAGGTTTCACCGATGATAAAGGCCTGCGCCAGACCGTCCGGGCT                                                                                                                                                                                                           >  3:567969/1‑132
                                                 GAAGAAGATGTTATCGCCCAGCACCAGACAAGAAGGTTCACCATTAAGGAAGGTTTCACCGATGATAAAGGCCTGCGCCAGACCGTCCGGGCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGG                                                                                                                                                                    <  1:552072/132‑1
                                                 GAAGAAGATGTTATCGCCCAGCACCAGACAAGAAGGTTCACCATTAAGGAAGGTTTCACCGATGATAAAGGCCTGCGCCAGACCGTCCGGGCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGG                                                                                                                                                                    >  3:552064/1‑132
                                                                 CCCAGGACCAGACAAGAAGGTTCACCATTAAGGAAGGTTTCACCGATGATAAAGGCATGCGCCAGACCGTCCGGGCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGC                                                                                                                                           <  1:161412/141‑1
                                                                                                                                             TGGCGCAGGCCTTTATCATCGGTGAAACCTGCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGG                                                                        <  3:619548/132‑1
                                                                                                                                               GCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCG  >  4:167311/1‑200
                                                                                                                                               GCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTT                                                             >  1:729552/1‑141
                                                                                                                                                                                                                                                                                                                                                       
CTTCGGTGCGCGCCGCAACATGACGCAGCTTCGGACTGAAGCCCTGACCGAAGAAGATGTTATCGCCCAGCACCAGACAAGAAGGTTCACCATTAAGGAAGGTTTCACCGATGATAAAGGCCTGCGCCAGACCGTCCGGGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3973992‑3973849
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctggCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCG  >  NC_000913/3973867‑3974065

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