New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 951207 =135 (1.990)5 (0.080) 5/250 NT 3.8% intergenic (‑127/+65) pflA/pflB pyruvate formate‑lyase 1‑activating enzyme; [formate‑C‑acetyltransferase 1]‑activating enzyme; PFL activase/formate C‑acetyltransferase 1, anaerobic; pyruvate formate‑lyase 1
?NC_000913 951233 = 125 (1.930)intergenic (‑153/+39) pflA/pflB pyruvate formate‑lyase 1‑activating enzyme; [formate‑C‑acetyltransferase 1]‑activating enzyme; PFL activase/formate C‑acetyltransferase 1, anaerobic; pyruvate formate‑lyase 1
Rejected: Frequency below/above cutoff threshold.

GTGTATCTGGCTACGCAGTACGTTTCAACTCGCTGACTAAAGAACAGCAGCAGGACGTTATTACTCGTACCTTCACTCAATCTATGTAATTAGATTTGACTGAAATCGTACAGTAAAAAGCGTACAATAAAGGCTCCACGAAAGTGGGGCCTTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/951360‑951207
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcctttATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGA  >  NC_000913/951233‑951377
                                                                                                                                                                                                                                                                                                           
GTGTATCTGGCTACGCAGTACGTTTCAACTCGCTGACTAAAGAACAGCAGCAGGACGTTATTACTCGTACCTTCACTCAATCTATGTAATTAGATTTGACTGAAATCGTACAGTAAAAAGCGTACAATAAAGGCTCCACGAAAGTGGGG                                                                                                                                                        >  6:151173/1‑149
     TCTGGCTACGCAGTACGTTTCAACTCGCTGACTAAAGAACAGCAGCAGGACGTTATTACTCGTACCTTCACTCAATCTATGTAATTAGATTTGACTGAAATCGTACAGTAAAAAGCGTACAATAAAGGCTCCACGAAAGTGGGGCCTTT                                                                                                                                                   >  1:375198/1‑149
         GCTACGCAGTACGTTTCAACTCGCTGACTAAAGAACAGCAGCAGGACGTTATTACTCGTACCTTCACTCAATCTATGTAATTAGATTTGACTGAAATCGTACAGTAAAAAGCGTACAATAAAGGCTCCACGAAAGTGGGGCCTTT                                                                                                                                                   <  5:151173/145‑1
                                      AAAGAACAGCAGCAGGACGTTATTACTCGTACCTTCACTCAATCTATGTAATTAGATTTGACTGAAATCGTACAGTAAAAAGCGTACAATAAAGGCTCCACGAAAGTGGGGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCA                                                                                                                   >  5:306546/1‑148
                                       AAGAACAGCAGCAGGACGTTATTACTCGTACCTTCACTCAATCTATGTAATTAGATTTGACTGAAATCGTACAGTAAAAAGCGTACAATAAAGGCTCCACGAAAGTGGGGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAA                                                                                                                 <  6:306546/149‑1
                                                                                                    TGAAATCGTACAGTAAAAAGCGTACAATAAAGGCTCCACGAAAGTGGGGCCTTTATTGTACGCTGTTTAATGTACGATTTAAGTCAAATAGAATTACATAGATTTAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGT                                                      >  7:276582/1‑147
                                                                                                                               TAAAGGCTCCACGAAAGTGGGGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGC                            >  5:260554/1‑146
                                                                                                                                  AGGCTCCACGAAAGTGGGGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATA                      >  6:338372/1‑149
                                                                                                                                                    GCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGA                                                                                   >  7:37998/1‑70
                                                                                                                                                    GCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGA                                                                                   <  8:37998/70‑1
                                                                                                                                                     CCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGG   <  5:338372/149‑1
                                                                                                                                                     CCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCG    >  3:168862/1‑148
                                                                                                                                                       TTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGA  >  1:225696/1‑148
                                                                                                                                                                                                                                                                                                           
GTGTATCTGGCTACGCAGTACGTTTCAACTCGCTGACTAAAGAACAGCAGCAGGACGTTATTACTCGTACCTTCACTCAATCTATGTAATTAGATTTGACTGAAATCGTACAGTAAAAAGCGTACAATAAAGGCTCCACGAAAGTGGGGCCTTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/951360‑951207
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcctttATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGA  >  NC_000913/951233‑951377

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