New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 639508 =39 (1.180)4 (0.120) 3/586 4.1 10.2% coding (564/564 nt) ahpC alkyl hydroperoxide reductase, C22 subunit
?NC_000913 639635 = 32 (0.990)intergenic (+127/‑118) ahpC/ahpF alkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below cutoff threshold.

TCCGAAGTATTTCGCGTCAGAGCAAGGCGGCTAGCCTGTGAATCCCCGGGAGCTTACACAAGTAAGTGACCGGGGTGAGCGGGCGACGCCAACGCCGCTATGGCGTGAAAGACGAAGGAAATT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/639630‑639508
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggaaattCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAGGAACTGTTGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGT  >  NC_000913/639635‑639931
                                                                                                                                                                                                                                                                                                                                                                                                                                    
TCCGAAGTATTTCGCGTCAGAGCAAGGCGGCTAGCCTGTGAATCCCCGGGAGCTTACACAAGTAAGTGACCGGGGTGAGCGGGCGACGCCAACGCCGCTATGGCGTGAAAGACGAAGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACCCAAATATGAAAACTAAACTAAAGGCTTACCTTGAGAAATTGACCAAGCC                                                                                                                          >  2:3202/1‑300
TCCGAAGTATTTCGCGTCAGAGCAAGGCGGCTAGCCTGTGAATCCCCGGGAGCTTACACAAGTAAGTGACCGGGGTGAGCGGGCGACGCCAACGCCGCTATGGCGTGAAAGACGAAGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGACGCCGCTCGCATGATGATGTTTAATGCCCAGGAGTTAAACCTGCTCGACATCATTTTGAAAACTACACTCACAGCTTACCTTGCAAAATTCACCAAca                                                                                                                           >  2:52360/1‑297
              CGTCAGAGCAAGGCGGCTAGCCTGTGAATCCCCGGGAGCTTACACAAGTAAGTGACCGGGGTGAGCGGGCGACGCCAACGCCGCTATGGCGTGAAAGACGAAGGAAATTCACCTAATGCTTCGGGTGCTGCGGCACCCGATGTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTG                                                                                                            <  1:52363/300‑1
                                  CCTGTGAATCCCCGGGAGCTTACACAAGTAAGTGACCGGGGTGAGCGGGCGACGCCAACGCCGCTATGGCGTGAAAGACGAAGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGC                                                                                         <  1:3202/299‑1
                                                                                                                       AATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAGGAACTGTTGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGT  >  1:256577/1‑301
                                                                                                                                                                                                                                                                                                                                                                                                                                    
TCCGAAGTATTTCGCGTCAGAGCAAGGCGGCTAGCCTGTGAATCCCCGGGAGCTTACACAAGTAAGTGACCGGGGTGAGCGGGCGACGCCAACGCCGCTATGGCGTGAAAGACGAAGGAAATT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/639630‑639508
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggaaattCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAGGAACTGTTGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGT  >  NC_000913/639635‑639931

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