New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi = 287764624 (1.020)3 (0.130) 3/272 3.1 10.8% intergenic (+28/+923) iap/ygbF aminopeptidase in alkaline phosphatase isozyme conversion/CRISPR adaptation ssRNA endonuclease
?NC_000913_3_pae_pgi = 2877662 27 (1.210)intergenic (+44/+907) iap/ygbF aminopeptidase in alkaline phosphatase isozyme conversion/CRISPR adaptation ssRNA endonuclease
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_pae_pgi/2877492‑2877646
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggagttcCCCCGCCTCTGCGGTAGAACTCCCAG  <  NC_000913_3_pae_pgi/2877662‑2877637
                                                                                                                                                                                     
GCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGCAATTTACAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGG                                  <  1:302396/149‑1
    GACGTAAGACTGGATAATCACCAACGTATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCGATGACGTTATGCGGATACTGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGACGGGTTTGAAAATGGGAGCCGGGAGT                              >  1:372431/1‑149
                                TTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCCCCCGCCTCTGCGGTAGAACTCCCAG  <  1:323445/149‑1
                                                                                                            GTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCCCCCGCCTCTGCGGTAGAACTCC     >  1:134342/1‑70
                                                                                                            GTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCCCCCGCCTCTGCGGTAGAACTCC     <  2:134342/70‑1
                                                                                                                                                                                     
GCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_pae_pgi/2877492‑2877646
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggagttcCCCCGCCTCTGCGGTAGAACTCCCAG  <  NC_000913_3_pae_pgi/2877662‑2877637

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