New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4586784 =153 (0.820)4 (0.020) 4/226 NT 2.5% intergenic (‑23/‑165) hsdR/mrr endonuclease R Type I restriction enzyme/methylated adenine and cytosine restriction protein
?NC_000913 4586809 = 169 (0.950)intergenic (‑48/‑140) hsdR/mrr endonuclease R Type I restriction enzyme/methylated adenine and cytosine restriction protein
Rejected: Frequency below/above cutoff threshold.

GCAGAATCAGGTAGATGTTTTTCGGCTACTATAGCACTACAAAAATAGACGAACACGTTAGAAATGAGTCAGTTGCTGTGACCGTGGTCATTGCCCGGAAAGGTACAGAAAGCTAAGATGAGATGTTATGGGCCTTAAATATTTGGACAGGCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4586937‑4586784
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aggcccATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGACGGT  >  NC_000913/4586809‑4586956
                                                                                                                                                                                                                                                                                                              
GCAGAATCAGGTAGATGTTTTTCGGCTACTATAGCACTACAAAAATAGACGAACACGTTAGAAATGAGTCAGTTGCTGTGACCGTGGTCATTGCCCGGAAAGGTACAGAAAGCTAAGATGAGATGTTATGGGCCTTAAATATTTGGACA                                                                                                                                                           <  2:1458075/149‑1
   GAATCAGGTAGATGTTTTTCGGCTACTATAGCACTACAAAAATAGACGAACACGTTAGAAATGAGTCAGTTGCTGTGACCGTGGTCATTGCCCGGAAAGGTACAGAAAGCTAAGATGAGATGTTATGGGCCTTAAATATTTGGACAGGC                                                                                                                                                        <  4:628968/149‑1
   GAATCAGGTAGATGTTTTTCGGCTACTATAGCACTACAAAAATAGACGAACACGTTAGAAATGAGTCAGTTGCTGTGACCGTGGTCATTGCCCGGAAAGGTACAGAAAGCTAAGATGAGATGTTATGGGCCTTAAATATTTGGACAGGC                                                                                                                                                        >  3:1503399/1‑149
     ATCAGGTAGATGTTTTTCGGCTACTATAGCACTACAAAAATAGACGAACACGTTAGAAATGAGTCAGTTGCTGTGACCGTGGTCATTGCCCGGAAAGGTACAGAAAGCTAAGATGAGATGTTATGGGCCTTAAATATTTGGACAGGCCC                                                                                                                                                      <  4:1037809/149‑1
     ATCAGGTAGATGTTTTTCGGCTACTATAGCACTACAAAAATAGACGAACACGTTAGAAATGAGTCAGTTGCTGTGACCGTGGTCATTGCCCGGAAAGGTACAGAAAGCTAAGATGAGATGTTATGGGCCTTAAATATTTGGACAGGCCC                                                                                                                                                      <  2:175060/149‑1
           TAGATGTTTTTCGGCTACTATAGCACTACAAAAATAGACGAACACGTTAGAAATGAGTCAGTTGCTGTGACCGTGGTCATTGCCCGGAAAGGTACAGAAAGCTAAGATGAGATGTTATGGGCCTTAAATATTTGGACA                                                                                                                                                           >  4:766802/1‑138
           TAGATGTTTTTCGGCTACTATAGCACTACAAAAATAGACGAACACGTTAGAAATGAGTCAGTTGCTGTGACCGTGGTCATTGCCCGGAAAGGTACAGAAAGCTAAGATGAGATGTTATGGGCCTTAAATATTTGGACA                                                                                                                                                           <  3:766802/138‑1
                               TAGCACTACAAAAATAGACGAACACGTTAGAAATGAGTCAGTTGCTGTGACCGTGGTCATTGCCCGGAAAGGTACAGAAAGCTAAGATGAGATGTTATGGGCCTTAAATATTTGGACAGG                                                                                                                                                         >  1:1732966/1‑120
                               TAGCACTACAAAAATAGACGAACACGTTAGAAATGAGTCAGTTGCTGTGACCGTGGTCATTGCCCGGAAAGGTACAGAAAGCTAAGATGAGATGTTATGGGCCTTAAATATTTGGACAGG                                                                                                                                                         <  2:1732966/120‑1
                                        AAAAATAGACGAACACGTTAGAAATGAGTCAGTTGCTGTGACCGTGGTCATTGCCCGGAAAGGTACAGAAAGCTAAGATGAGATGTTATGGGCCTTAAATATTTGGACAGGC                                                                                                                                                        >  3:548447/1‑112
                                        AAAAATAGACGAACACGTTAGAAATGAGTCAGTTGCTGTGACCGTGGTCATTGCCCGGAAAGGTACAGAAAGCTAAGATGAGATGTTATGGGCCTTAAATATTTGGACAGGC                                                                                                                                                        <  4:548447/112‑1
                                                                                                                             TTATGGGCCTTAAATATTTGGACAGGCCCATAACATCTTATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGA                                                                                         >  3:1177881/1‑90
                                                                                                                             TTATGGGCCTTAAATATTTGGACAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGA                                                                                         <  4:1177881/90‑1
                                                                                                                                      TTAAATATTTGGACAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGAC                                                                                                           <  2:394537/63‑1
                                                                                                                                      TTAAATATTTGGACAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGAC                                                                                                           >  1:394537/1‑63
                                                                                                                                                      GCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGAC     <  3:1238618/149‑1
                                                                                                                                                       CCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGACG    >  2:1865817/1‑149
                                                                                                                                                       CCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGACG    >  2:1282427/1‑149
                                                                                                                                                        CCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGC                                                                                                                 >  3:1042989/1‑39
                                                                                                                                                        CCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGC                                                                                                                 <  4:1042989/39‑1
                                                                                                                                                         CATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGACGGT  >  4:1626444/1‑149
                                                                                                                                                         CATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGACGGT  >  1:14614/1‑149
                                                                                                                                                         CATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTA            <  3:1305645/139‑1
                                                                                                                                                         CATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTA            >  4:1305645/1‑139
                                                                                                                                                                                                                                                                                                              
GCAGAATCAGGTAGATGTTTTTCGGCTACTATAGCACTACAAAAATAGACGAACACGTTAGAAATGAGTCAGTTGCTGTGACCGTGGTCATTGCCCGGAAAGGTACAGAAAGCTAAGATGAGATGTTATGGGCCTTAAATATTTGGACAGGCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4586937‑4586784
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aggcccATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTTTTGTAGTGCTATAGTAGCCGAAAAACATCTACCTGATTCTGCAAGGATGTACTATGACGGT  >  NC_000913/4586809‑4586956

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