New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_tpiA 2920036 =21 (0.840)4 (0.170) 4/272 2.9 16.3% noncoding (2920036/4641682 nt) repeat_region repeat region
?NC_000913_3_vch_tpiA 2920054 = 21 (0.890)noncoding (2920054/4641682 nt) repeat_region repeat region
Rejected: Frequency below/above cutoff threshold.

GCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCTTTAAAGCCGTAGCGATCCTGTGAGGCTTCC  >  NC_000913_3_vch_tpiA/2919890‑2920174
                                                                                                                                                  |                                                                                                                                          
acagcgaagtcgttgtgcgtctggcggaagcctcacaggatagctacggctttaaagatttcaaacttaagggcggcgtgttacctggcgagcaagaaatcgacactgttcgtgcattgaagaaacgcttcccggatgcgcggattAcc                                                                                                                                          >  1:252213‑M1/147‑149 (MQ=255)
     gaagccgttgtgcgtctggcggaagcctcacaggatcgctacggctttaaagatttcaaacttaagggcggcgtgttacctggcgagcaagaaatcgacactgttcgtgcattgaagaaacgcttcccggatgcgcgcattACCGTtgc                                                                                                                                     <  2:208644‑M1/8‑2 (MQ=255)
           gCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCaa                                                                                                                               <  1:402428/149‑1 (MQ=255)
                          aaaTTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCAt                                                                                                                >  1:326614/1‑149 (MQ=255)
                                cGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCTGCATAGGTAAGAACATTATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGa                                                                                                          <  1:12163/149‑1 (MQ=255)
                                        cTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTCAGAACATCCTTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTttc                                                                                                  <  1:202629/149‑1 (MQ=255)
                                             cGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCaa                                                                                             <  2:326614/149‑1 (MQ=255)
                                               tCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATg                                                                                           >  2:411995/1‑149 (MQ=255)
                                                  ggAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCAc                                                                                        >  1:389107/1‑149 (MQ=255)
                                                  ggAGAAGCCCTGTTCTGCGCAGCATGGATCTTCGGCATGGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATAAAGCAGCCATGCACCGATGTGCTTAACGTTAATCTGCGCATCCGGACAGCGTTTCTTCAATacac                                                                                        >  2:393054/1‑149 (MQ=255)
                                                         cccTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAgt                                                                                 >  2:283525/1‑149 (MQ=255)
                                                           cTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAgtgt                                                                               <  2:393955/149‑1 (MQ=255)
                                                            tGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTc                                                                              >  1:302947/1‑149 (MQ=255)
                                                               aacttaagggcggcgtgttacctggcgagcaagaaatcgacactgttcgtgcattgaagaaacgcttcccggatgcgcggattACCGTTGGGATCAACGGTAATc                                                                                                                       <  1:184039‑M1/22‑1 (MQ=255)
                                                               aacttaagggcggcgtgttacctggcgagcaagaaatcgacactgttcgtgcattgaagaaacgcttcccggatgcgcggattACCGTTGGGATCAACGGTAATc                                                                                                                       >  2:184039‑M1/84‑105 (MQ=255)
                                                                            gATCTTCGGCATATGTTAGTACATCATTCAGCTTTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCa                                                              <  1:196601/149‑1 (MQ=255)
                                                                                           tAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGccc                                               <  1:283525/149‑1 (MQ=255)
                                                                                             aGAACATCATTCAGCCTATTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCtt                                             <  2:344944/149‑1 (MQ=255)
                                                                                                   tcgacactgttcgtgcattgaagaaacgcttcccggatgcgcggattACCGTTGGGATCAACGGTAATCCGCGCATCCGGGa                                                                                                          <  1:362654‑M1/35‑1 (MQ=255)
                                                                                                   tcgacactgttcgtgcattgaagaaacgcttcccggatgcgcggattACCGTTGGGATCAACGGTAATCCGCGCATCCGGGa                                                                                                          >  2:362654‑M1/48‑82 (MQ=255)
                                                                                                   tcatTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGc                                                                                                       <  1:395213/85‑1 (MQ=255)
                                                                                                   tcatTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGc                                                                                                       >  2:395213/1‑85 (MQ=255)
                                                                                                     atTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGCTAACACGCCGCCCTTAAGTTTGa                                     >  2:32134/1‑149 (MQ=255)
                                                                                                        cAGCCCTTTGCATAAAGAACTGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGAGTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGaaag                                  >  2:380848/1‑148 (MQ=255)
                                                                                                                    tAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCTTTAAAGCCGt                      >  2:81641/1‑149 (MQ=255)
                                                                                                                               cTTCATCAAGCAGCCATGCACCGTTGGGCTCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCTTTAAAGCCGTAGCGATCCtgt           >  2:205491/1‑149 (MQ=255)
                                                                                                                                        gcagcCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCTTTAAAGCCGTAGCGATCCTGTGAGGCTTcc  <  2:302947/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
GCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGGTAATCCGCGCATCCGGGAAGCGTTTCTTCAATGCACGAACAGTGTCGATTTCTTGCTCGCCAGGTAACACGCCGCCCTTAAGTTTGAAATCTTTAAAGCCGTAGCGATCCTGTGAGGCTTCC  >  NC_000913_3_vch_tpiA/2919890‑2920174

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.