New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_pgi = 72066718 (0.980)3 (0.190) 3/252 2.1 15.5% coding (85/258 nt) ybfK uncharacterized protein
?NC_000913_3_vch_pgi = 720649 17 (1.060)coding (67/258 nt) ybfK uncharacterized protein
Rejected: Frequency below/above cutoff threshold.

TCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATGG  >  NC_000913_3_vch_pgi/720503‑720796
                                                                                                                                                  |                                                                                                                                                   
tcGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAAta                                                                                                                                                   >  1:323491‑M2/1‑147 (MQ=255)
                  gTGCCAGAAAAAAGCAAAGTTTCCCCTCTGATAAATGAATGAGGTAGATTTCTTTTTTCAGAATGCAGAAACCCTTACCATTCAGTTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGTCATATGCCTGTAtt                                                                                                                                 <  2:78783‑M1/149‑3 (MQ=255)
                    gCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATAc                                                                                                                               <  1:63593/149‑1 (MQ=255)
                           aaaaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGt                                                                                                                        >  1:19251/1‑149 (MQ=255)
                              aGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAg                                                                                                                     <  1:209340/149‑1 (MQ=255)
                                                ataaataaATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCGGTACGCTGaa                                                                                                   >  1:32237/1‑149 (MQ=255)
                                                     taaatCAGGTATATCGCTTTTTTCAGAATTCAGAAACCCTTACCCGTCAGGTGGCTGATAATACAAGGGGTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGcc                                                                                              <  1:19421/149‑1 (MQ=255)
                                                                              gAATGCAGAAACCCTTACCAGGCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACATGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTGGTCGGTAGACTAACGGGCAtt                                                                     <  2:307801/149‑1 (MQ=255)
                                                                                         cccTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTTAAGAAGTGaa                                                          <  1:153613/149‑1 (MQ=255)
                                                                                               ccAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGa                                                    >  2:23341/1‑149 (MQ=255)
                                                                                                             gATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGGAGAAGTCAAAGTCGAAACTATTGCGAttg                                      >  2:27047/1‑147 (MQ=255)
                                                                                                             gATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGa                                      <  1:23341/149‑1 (MQ=255)
                                                                                                                  taCAAGGGTTAACAGAGAATGAGCACTGCGGTATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCaaa                                 <  1:311730/149‑1 (MQ=255)
                                                                                                                  taCAAGGGTTAACAGAGAATGAGCACTGCGGAAtacaggcatatgcctgtattccgcagtgctcattctctgttaacc                                                                                                        >  2:118756‑M2/1‑33 (MQ=255)
                                                                                                                  taCAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAttccgcagtgctcattctctgttaacc                                                                                                        >  1:118756‑M1/1‑51 (MQ=255)
                                                                                                                     aaGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAATTGCCCTTAGTCGGTAGACTAACGGGCATCGAAGAAGTGAAAGTCGAAGCTATTGCGAt                                        <  1:238763/139‑1 (MQ=255)
                                                                                                                     aaGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATCGAAGAAGTGAAAGTCGAAGCTATTGCGAt                                        >  2:238763/1‑139 (MQ=255)
                                                                                                                        ggTTAACAGAGAATGAGCACTGCGGAAtacaggcatatgcctgtattccgcagtgctcattctctgttaacccttgta                                                                                                  >  1:118756‑M2/1‑27 (MQ=255)
                                                                                                                        ggTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAttccgcagtgctcattctctgttaacccttgta                                                                                                  >  2:118756‑M1/1‑45 (MQ=255)
                                                                                                                        ggTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGTCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCat                           <  2:52901/149‑1 (MQ=255)
                                                                                                                               agagAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCataa                         <  2:289054/144‑1 (MQ=255)
                                                                                                                               agagAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCataa                         >  1:289054/1‑144 (MQ=255)
                                                                                                                                      gAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTcc             >  1:231652/1‑149 (MQ=255)
                                                                                                                                                 aatacaggcatatgactgtattccgcagtgctcattctctgttaacccttgtattatcagccaactgaatggtaagggtttctgcattctgaaaaaagaaatctacctcattcatttatcagaggggaaactttgcttttttctggcac  <  2:78783‑M2/149‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATGG  >  NC_000913_3_vch_pgi/720503‑720796

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