New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 458738686 (0.870)3 (0.030) 3/268 NT 3.2% coding (438/915 nt) mrr methylated adenine and cytosine restriction protein
?NC_000913 = 4587414 98 (1.030)coding (466/915 nt) mrr methylated adenine and cytosine restriction protein
Rejected: Frequency below/above cutoff threshold.

GCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4587233‑4587386
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcagCCACCGCATCACGAAGCTCTTTTAGCGCCTGATCTAACCGATCGTCCGGGCTGCTCTTTGCAAGTTCTTCATGATCGGGAGAGTCGGCTTTCGGATCTAAATCGACGGCATCCGGCCGTGACTTAAGTTTGACATTCACAAAAGCGAAG  <  NC_000913/4587414‑4587266
                                                                                                                                                                                                                                                                                                               
GCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGTATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGCTCGGTTAGA                                                                                                                                                           >  2:1307961/1‑150
GCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGA                                                                                                                                                           >  1:1526756/1‑150
 CCATTGACGGAGCAGGAGACGAACCATGTGGCCTTCGCTTTTGTTAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTATATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGAT                                                                                                                                                          <  1:1630912/150‑1
 CCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTTTGATTGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTATATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCGGGACGATCGGTTAGAT                                                                                                                                                          <  2:271352/150‑1
  CAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATC                                                                                                                                                         >  1:733294/1‑150
   AATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCAGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAATAGCATCCCGGACGATCGGTTAGAACA                                                                                                                                                        >  2:665307/1‑150
                                   aCGCTTTTCTGAATGCCACACTTACGTCACGCCCGGAAGCCGTCGATTTACATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGCCACCGCATCACGAAGCTCTTTTAGCGCCTG                                                                                                                        <  2:215885/149‑1
                                                     AACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGCCACCGCATCACGAAGCTCTTTTAGCGCCTGATCTAACC                                                                                                                <  1:260150/140‑1
                                                     AACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGCCACCGCATCACGAAGCTCTTTTAGCGCCTGATCTAACC                                                                                                                >  2:260150/1‑140
                                                                                                                                                    aATCAGACATCGCATCAATAAGCTCTTTTAGCGTCTGATGTAAACGAGCGTCAGTGCTTCTCTTTGCAAGTTCTTCATGATCGGGAGAGTCGGCTTTCGGATCTAAATAGTCGGCATCCGGCCGTGACTTAAGTTTGACATTCACAAAAG       <  1:1203108/149‑1
                                                                                                                                                     ATCAGCCACCGCATCACGAAGCTCTTTTAGCGCCTGATCTAACCGATCGTCCGGGCTGCTCTTTGCAAGTTCTTCATGATCGGGAGAGTCGGCTTTCGGATCTAAATCGACGGCATCCGGCCGTGACTTAAGTTTGACATTCACAAAAGC      <  1:920877/150‑1
                                                                                                                                                       CAGCCACCGCATCACGAAGCTCTTTTAGCGCCTGATCTAACCGATCGTCCGGGCTGCTCTTTGCAAGTTCTTCATGATCGGGAGAGTCGGCTTTCGGATCTAAATCGACGGCATCCGGCCGTGACTTAAGTTTGACATTCACAAAAGCGA    <  1:1091478/150‑1
                                                                                                                                                        AGCCACCGCATCACTAAGCTCTTTTAGCGCCTGATCTAACCGATCGTCCGGGCTGCTCTTTGCAAGTTCTTCATGATCGGGAGAGTCGGCTTTCGGATCTAAATCGACGGCATCCGGCCGTGACTTAAGTTTGACATTCACTAAAGCGAA   <  2:927032/150‑1
                                                                                                                                                         GCCACCGCATCACGAAGCTCTTTTAGCGCCTGATCTAACCGATCGTCCGGGCTGCTCTTTGCAAGTTCTTCATGATCGGGAGAGTCGGCTTTCGGATCTAAATCGACGGCATCCGGCCGTGACTTAAGTTTGACATTCACAAAAGCGAAG  <  1:963205/150‑1
                                                                                                                                                                                                                                                                                                               
GCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4587233‑4587386
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcagCCACCGCATCACGAAGCTCTTTTAGCGCCTGATCTAACCGATCGTCCGGGCTGCTCTTTGCAAGTTCTTCATGATCGGGAGAGTCGGCTTTCGGATCTAAATCGACGGCATCCGGCCGTGACTTAAGTTTGACATTCACAAAAGCGAAG  <  NC_000913/4587414‑4587266

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