New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1662393 =64 (0.570)6 (0.060) 5/264 NT 8.9% coding (1838/2424 nt) ynfF putative selenate reductase YnfF
?NC_000913 1662430 = 61 (0.560)coding (1875/2424 nt) ynfF putative selenate reductase YnfF

TCACGGAATGCGCGGAAAGCGACGTAGTGTTCTTCCGGGCATT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1662435‑1662393
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcattCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCG  >  NC_000913/1662430‑1662559
                                                                                                                                                                             
TCACGGAATGCGCGGAAAGCGACGTAGTGTTCTTCCGGGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATG                         >  2:29057/1‑150
 CACGGAATGCGCGGAAAGCGACGTAGTGTTCTTCCGGGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGA                        >  1:1833734/1‑150
   CGGAATGCGCGGAAAGCGACGTAGTGTTCTTCCGGGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAA                      >  1:1056381/1‑150
   CGGAATGCGCGGAAAGCGACGTAGTGTTCTTCCGGGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAA                      >  1:1905252/1‑150
      AATGCGCGGAAAGCGACGTAGTGTTCTTCCGGGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTA                   >  1:1152937/1‑150
                       GTAGTGTTCTTCCGGGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCG  <  2:1833734/150‑1
                                                                                                                                                                             
TCACGGAATGCGCGGAAAGCGACGTAGTGTTCTTCCGGGCATT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1662435‑1662393
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcattCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCG  >  NC_000913/1662430‑1662559

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.