New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 = 4468461116 (1.420)4 (0.050) 3/96 NT 3.6% coding (929/1656 nt) treC trehalose‑6‑P hydrolase
?AP009048 = 4468484 100 (1.250)coding (906/1656 nt) treC trehalose‑6‑P hydrolase
Rejected: Frequency below/above cutoff threshold.

AAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/4468391‑4468461
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cATTGTTCCGCCACTGGCAACAAGGAATGCACAACGTAGCATGGAATGCCTTGTTCTGGTGTAACCACGAT  <  AP009048/4468484‑4468415
                                                                                                                                             
AAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTC                                                                        >  1:2503804/1‑71
                        ATCGTGGTTACACCAGAACAAGGCATTCAATGCTACGTTGTGCATTCATTGTTCCGCCACTGGCAACAAGG                                                <  1:594522/71‑1
                          CGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTC                                                                        >  1:3682494/1‑45
                          CGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTC                                                                        >  1:3793626/1‑45
                                            AGGCATTCCATGCTACGTTGTGCATTCATTGTTCCGCCACTGGCAACAAGGAA                                              >  1:1485057/1‑53
                                            AGGCATTCCATGCTACGTTGTGCATTCATTGTTCCGCCACTGGCAACAAGGAA                                              >  1:5830711/1‑53
                                                  TCCATGCTACGTTGTGCATTCATTGTTCCGCCACTGG                                                        <  1:712812/37‑1
                                                                      CATTGTTCCGCCACTGGCAACAAGGAATGCACAACGTAGCATGGAATGCCTTGTTCTGGTGTAACCACGAT  >  1:2113059/1‑71
                                                                      CATTGTTCCGCCACTGGCAACAAGGAATGCACAACGTAGCATGGAATGC                        >  1:2188541/1‑49
                                                                      CATTGTTCCGCCACTGGCAACAAGGAATGCACAACGTAGCATGGAATGC                        >  1:3241408/1‑49
                                                                      CATTGTTCCGCCACTGGCAACAAGGAATGCACAACGTAGCATGGAATGC                        >  1:5975530/1‑49
                                                                      CATTGTTCCGCCACTGGCAACAAGGAATGCACAACGTAGCATGGAATGC                        >  1:771865/1‑49
                                                                      CATTGTTCCGCCACTGGCAACAAGGAATGCACAACGTAGCATGGAATG                         >  1:2952790/1‑48
                                                                      CATTGTTCCGCCACTGGCAACAAGGAATGCACAACGTAGCATGGAATG                         >  1:4013593/1‑48
                                                                      CATTGTTCCGCCACTGGCAACAAGGAATGCACAACGTAGCATGGAATG                         >  1:4262606/1‑48
                                                                      CATTGTTCCGCCACTGGCAACAAGGAATGCACAACGTAGCATGGAATG                         >  1:6577832/1‑48
                                                                                                                                             
AAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/4468391‑4468461
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cATTGTTCCGCCACTGGCAACAAGGAATGCACAACGTAGCATGGAATGCCTTGTTCTGGTGTAACCACGAT  <  AP009048/4468484‑4468415

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.