New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3239599 =44 (1.180)0 (0.000)
+TCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCT
6/176 2.8 0.0% intergenic (+54/‑345) alx/sstT putative membrane‑bound redox modulator/sodium:serine/threonine symporter
?NC_000913 3504854 = 34 (0.910)intergenic (+71/+81) frlR/yhfS putative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Rejected: Frequency below/above cutoff threshold.

AGGACTCGGTCAGTCCCCTCGCCCCGTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCAGCATCGGTGCGAAAAATTACCCTCACCCCGACCCTCTCCCTGGAAGGGAGAGGG‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3239713‑3239599
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3504854‑3504827
                                                                                                                   |||||||||| 
AGGACTCGGTCAGTCCCCTCGCCCCGTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCAGCATCGGTGCGAAAAATTACCCTCACCCCGACCCTCTCCCTTCCAGGGAGAGGGTCGGG        >  7:6448/1‑120
AGGACTCGGTCAGTCCCCTCGCCCCGTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCAGCATCGGTGCGAAAAATTACCCTCACCCCGACCCTCTCCCTTCCAGGGAGAGGGTCGGG        <  8:6448/120‑1
             TCCCCTCGCCCCGTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCAGCATCGGTGCGAAAAATTACCCTCACCCCGACCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGG  >  1:66503/1‑113
             TCCCCTCGCCCCGTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCAGCATCGGTGCGAAAAATTACCCTCACCCCGACCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGG  <  2:66503/113‑1
                               GGAGAGGGTTAGGGTGAGGGGAACAGGCCAGCATCGGTGCGAAAAATTACCCTCACCCCGACCCTCTCCCTTCCAGGGAGAGGGTCGGG        >  7:26976/1‑89
                               GGAGAGGGTTAGGGTGAGGGGAACAGGCCAGCATCGGTGCGAAAAATTACCCTCACCCCGACCCTCTCCCTTCCAGGGAGAGGGTCGGG        <  8:26976/89‑1
                                                 GGGAACAGGCCAGCATCGGTGCGAAAAATTACCCTCACCCCGACCCTCTCCCTTCCAGGGAGAGGGT            <  5:134665/67‑1
                                                 GGGAACAGGCCAGCATCGGTGCGAAAAATTACCCTCACCCCGACCCTCTCCCTTCCAGGGAGAGGGT            >  6:134665/1‑67
                                                                                                                   |||||||||| 
AGGACTCGGTCAGTCCCCTCGCCCCGTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCAGCATCGGTGCGAAAAATTACCCTCACCCCGACCCTCTCCCTGGAAGGGAGAGGG‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3239713‑3239599
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3504854‑3504827

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.