New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3548549102 (0.900)5 (0.060)
+TTTTTTTTTAATGATACGGCGACC
4/150 13.5 6.1% coding (1522/2085 nt) malQ 4‑alpha‑glucanotransferase (amylomaltase)
?NC_000913 3548574 = 103 (0.910)coding (1497/2085 nt) malQ 4‑alpha‑glucanotransferase (amylomaltase)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below cutoff threshold.

TTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3548488‑3548549
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCG  >  NC_000913/3548574‑3548632
                                                              ||||||||||||||||||||||||                                                           
TTCTCGTGGTCGTTTTCGAAATAGAGCATTTTTTAAGAGAACACACCGCCGCTGCGCAGCTTTTTTTTTTTAATGATACGGCGACCACCGAGATCTACAC                                               <  1:1623229/100‑1
 ttttgttGTCGTTTTCTAAATAGAGCATTTTGTATGAGTACACACCGCTGCTGCGCAGTTTTTTTTTTTTAATGATACGGCGACCACCGAGATCTACACC                                              <  1:1842226/93‑1
                CGAAATAGAGCACTTTGTAAGAGTACGCACCGCTGCTGCGCAGTTTTTTTTTTTTAATGATACGGCGACCACCGAGATCTACACCAATCACCATACAGCG                               <  1:1953998/100‑1
                       GAGCTCTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTTTTTTTTTTAATGATACGGCGACCACCGAGATCTACACCAATCACCATACAGCGATGACG                         <  1:2631950/99‑1
                                             CCGCTGCTGCGCAGCTTATTGTTTTTAATGATACGTCGACCACCGAGATCTACACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCG  <  1:1134258/100‑1
                                                              ||||||||||||||||||||||||                                                           
TTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3548488‑3548549
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCG  >  NC_000913/3548574‑3548632

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