New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 4176601 =71 (1.210)3 (0.060) 3/90 NT 4.5% coding (892/984 nt) gspK general secretory pathway component, cryptic
?AP009048 4176625 = 64 (1.240)coding (868/984 nt) gspK general secretory pathway component, cryptic
Rejected: Frequency below/above cutoff threshold.

CCACTGGTTGCTCAGGTGAAAAGGCATCTCTTTCCCTATAGCCGCTATTTTACGCTGAGTACGGAAAGTATCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/4176676‑4176601
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcagcgTAAAATAGCGGCTATAGGGAAAGAGATGCCTTTTCACCTGAGCAACCAGTGGCTTCACACCGGAGAAATC  >  AP009048/4176625‑4176694
                                                                                                                                                  
CCACTGGTTGCTCAGGTGAAAAGGCATCTCTTTCCCTATAGCCGCTATTTTACGCTGAGTACGGAAAGTAT                                                                             <  1:898924/71‑1
CCACTGGTTGCTCAGGTGAAAAGGCATCTCTTTCCCTATAGCCGCTATTTTACGCTGAGTACGGAAAGTAT                                                                             <  1:236756/71‑1
  ACTGGTTGCTCAGGTGAAAAGGCATCTCTTTCCCTATAGCCGCTATTTTACGCTGAGTACGGAAAGTATC                                                                            >  1:3356633/1‑70
    TGGTTGCTCAGGTGAAAAGGCATCTCTTTCCCTATAGCCGCTATTTTACGCTGAGTACGGAAAGTATCAGC                                                                         >  1:3088997/1‑71
    TGGTTGCTCAGGTGAAAAGGCATCTCTTTCCCTATAGCCGCTATTTTACGCTGAGTACGGAAAGTATCAGC                                                                         >  1:4499384/1‑71
      GTTGCTCAGGTGAAAAGGCATCTCTTTCCCTATAGCCGCTATTTTACGCTGAGTACGGAAAGTATCAGC                                                                         <  1:2684831/69‑1
      GTTGCTCAGGTGAAAAGGCATCTCTTTCCCTATAGCCGCTATTTTACGCTGAGTACGGAAAGTAT                                                                             >  1:3638252/1‑65
            CAGGTGAAAAGGCATCTCTTTCCCTATAGCCGCTATTTTACGCTGAGTACGGAAAGTATCA                                                                           >  1:4927863/1‑61
                TGAAAAGGCATCTCTTTCCCTATAGCCGCTATTTTACGCTGAGTACGGAAAGTAT                                                                             <  1:428285/55‑1
                               TTCCCTATAGCCGCTATTTTACGCTGAGTACGGAAAGTATCAG                                                                          >  1:4050676/1‑43
                                     ATAGCCGCTATTTTACGCTGAGTACGGAAAGTATC                                                                            >  1:1791861/1‑35
                                                     GCTGAGTACGGAAAGTATCAGCGTAAAATAGCGGCTATAGGGAAAGAGATGCCTTTTCACCT                                 <  1:5061744/62‑1
                                                     GCTGAGTACGGAAAGTATCAGCGTAAAATAGCGGCTATAGGGAAAGAGAT                                             <  1:4458178/50‑1
                                                       TGAGTACGGAAAGTATCAGCGTAAAATAGCGGCTATAGGGA                                                    <  1:3068501/41‑1
                                                                      TCAGCGTAAAATAGCGGCTATAGGGAAAGAGATGCCTTTTCACCTGAGCA                            <  1:3744828/50‑1
                                                                      TCAGCGTAAAATAGCGGCTATAGGGAAAGAGATGCCTTTTCAC                                   <  1:5022219/43‑1
                                                                      TCAGCGTAAAATAGCGGCTATAGGGAAAGAGATGCCTTTTC                                     <  1:2823352/41‑1
                                                                       CAGCGTAAAATAGCGGCTATTGGGAAAGAGATGCCTTTTCACCTGAGCAACCAG                       <  1:4995638/54‑1
                                                                        AGCGTAAAATAGCGGCTATAGGGAAAGAGATGCCTTTTCACCTGAGCAACCAGTGGCTTCACACCGGAGAA     <  1:4043947/71‑1
                                                                        AGCGTAAAATAGCGGCTATAGGGAAAGAGATGCCTTTTCACCTGAGCAACCAGTG                     >  1:2425729/1‑55
                                                                        AGCGTAAAATAGCGGCTATAGGGAAAGAGATGCCTTTTCACCTGAGCAA                           >  1:146776/1‑49
                                                                           GTAAAATAGCGGCTATAGGGAAAGAGATGCCTTTTCACCTGAGCAACCAGTGGCTTCACACCGGAGAAATC  <  1:1597629/71‑1
                                                                                                                                                  
CCACTGGTTGCTCAGGTGAAAAGGCATCTCTTTCCCTATAGCCGCTATTTTACGCTGAGTACGGAAAGTATCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/4176676‑4176601
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcagcgTAAAATAGCGGCTATAGGGAAAGAGATGCCTTTTCACCTGAGCAACCAGTGGCTTCACACCGGAGAAATC  >  AP009048/4176625‑4176694

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