New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4146335 =84 (1.150)3 (0.040) 3/252 NT 3.6% coding (78/879 nt) pflC putative [formate‑C‑acetyltransferase 2]‑activating enzyme; pyruvate formate‑lyase 1‑activating enzyme
?NC_000913 4146358 = 83 (1.200)coding (101/879 nt) pflC putative [formate‑C‑acetyltransferase 2]‑activating enzyme; pyruvate formate‑lyase 1‑activating enzyme
Rejected: Frequency below/above cutoff threshold.

AAAAGACCACCGTACGAATGCCCTCGCCGTCATTCAGTGAATAACGCTGAATG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4146387‑4146335
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgaatgACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGA  >  NC_000913/4146358‑4146502
                                                                                                                                                                                                      
AAAAGACCACCGTACGAATGCCCTCGCCGTCATTCAGTGAATAACGCTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCA                                                     <  5:265990/147‑1
AAAAGACCACCGTACGAATGCCCTCGCCGTCATTCAGTGAATAACGCTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCA                                                     >  6:265990/1‑147
                           CGTCATTCAGTGAATAACGCTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTG                         >  5:329512/1‑148
                                              CTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGACCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCG     >  7:121473/1‑149
                                               TGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGT    >  8:105430/1‑149
                                               TGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGACCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGT    <  8:121473/149‑1
                                                 AATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGA  <  7:138970/149‑1
                                                   TGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAG                                                                         <  5:230034/76‑1
                                                   TGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAG                                                                         >  6:230034/1‑76
                                                                                                                                                                                                      
AAAAGACCACCGTACGAATGCCCTCGCCGTCATTCAGTGAATAACGCTGAATG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4146387‑4146335
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgaatgACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGA  >  NC_000913/4146358‑4146502

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