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.

GTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATGGGGCATAGTGCGGCTGTTA  >  NC_000913_3_vch_pgi/720521‑720814
                                                                                                                                                  |                                                                                                                                                   
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                                                                                                                          >  1:118756‑M1/1‑51 (MQ=255)
                                                                                                   aaGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAATTGCCCTTAGTCGGTAGACTAACGGGCATCGAAGAAGTGAAAGTCGAAGCTATTGCGAt                                                          <  1:238763/139‑1 (MQ=255)
                                                                                                   aaGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATCGAAGAAGTGAAAGTCGAAGCTATTGCGAt                                                          >  2:238763/1‑139 (MQ=255)
                                                                                                      ggTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAttccgcagtgctcattctctgttaacccttgta                                                                                                                    >  2:118756‑M1/1‑45 (MQ=255)
                                                                                                      ggTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGTCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCat                                             <  2:52901/149‑1 (MQ=255)
                                                                                                             agagAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCataa                                           >  1:289054/1‑144 (MQ=255)
                                                                                                             agagAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCataa                                           <  2:289054/144‑1 (MQ=255)
                                                                                                                    gAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTcc                               >  1:231652/1‑149 (MQ=255)
                                                                                                                                    aGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATGGGGCAt               >  1:167535/1‑149 (MQ=255)
                                                                                                                                                 tAttccgcagtgctcattctctgttaacccttgtattatcagccacctgactggtaagggtttctgcattctgaaaaaagcaatctacctgatttatttatcagcggggtaactttgcttttttccggcacgatcgatttctcattcga  >  1:323491‑M1/1‑2 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGAGTTGTCCGCCTTATATGGGGCATAGTGCGGCTGTTA  >  NC_000913_3_vch_pgi/720521‑720814

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.