New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 463832982 (0.890)8 (0.090) 5/168 NT 8.8% coding (152/1353 nt) creD inner membrane protein
?NC_000913 = 4638348 89 (1.030)coding (171/1353 nt) creD inner membrane protein
Rejected: Frequency below cutoff threshold.

GTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4638226‑4638329
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcggCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCAT  <  NC_000913/4638348‑4638250
                                                                                                                                                                                                           
GTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCA                                                                                                         >  2:397281/1‑100
 TATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCA                                                                                                         >  2:260870/1‑99
  ATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGC                                                                                                       >  2:659503/1‑100
   tgGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGCTTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGC                                                                                                       <  1:745927/97‑1
   TTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCG                                                                                                      >  2:1116125/1‑100
                           CTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGCCCAACGAGTTTTTGCGGCCCGC                                                                              <  1:1234769/100‑1
                            TGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGCCCAACGAGTTTTTGCGGCCCGCT                                                                             <  1:1290112/100‑1
                            TGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGCCCAACGAGTTTTTGCGGCCCGCT                                                                             <  2:651539/100‑1
                             GATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGCCCAACGAGTTTTTGCGGCCCGCTn                                                                            <  2:976917/100‑2
                             GATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGCCCAACGAGTTTTTGCGGCCCGCTG                                                                            <  2:2314367/100‑1
                              ATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGCCCAACGAGTTTTTGCGGCCCGCTGG                                                                           <  1:1667830/100‑1
                               TTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGT                                                                          <  2:2021997/100‑1
                                TCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGCCCAACGAGTTTTTGCGGCCCGCTGGT                                                                          <  1:1246513/99‑1
                                                                                                   AGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCA      >  2:783376/1‑100
                                                                                                   AGCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAAT        >  2:1041778/1‑98
                                                                                                    GCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAG     <  2:179718/100‑1
                                                                                                    GCGGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAG     <  1:2506233/100‑1
                                                                                                      GGCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCA   >  2:697574/1‑100
                                                                                                       GCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCAT  >  2:1301540/1‑100
                                                                                                                                                                                                           
GTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4638226‑4638329
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcggCCCAACGAGTTTTTGCGGCCCGCTGGTACTTTGGCGAATCGCATCTTCCACATCGCTACGGTAATCAGCACGTTCGACAATCACCTGCCGAATCAGCAT  <  NC_000913/4638348‑4638250

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