New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_tpiA = 260329526 (0.800)5 (0.160) 5/268 3.6 17.3% noncoding (2603295/4641682 nt) repeat_region repeat region
?NC_000913_3_vch_tpiA = 2603299 23 (0.750)noncoding (2603299/4641682 nt) repeat_region repeat region
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCT  >  NC_000913_3_vch_tpiA/2603148‑2603424
                                                                                                                                                   |                                                                                                                                 
tAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGt                                                                                                                                  <  2:26994/149‑1 (MQ=255)
   tgtgCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCTGGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTtgg                                                                                                                               <  1:178827/149‑1 (MQ=255)
        ttCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTtggctggc                                                                                                                          <  1:516920/149‑1 (MQ=255)
            tAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCg                                                                                                                      <  2:237168/149‑1 (MQ=255)
             aaacaTGTACGTTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGTTCGTCGCTATGTTACTGCTTGTGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGa                                                                                                                     <  1:116036/149‑1 (MQ=255)
                  tGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATa                                                                                                                <  2:378884/149‑1 (MQ=255)
                   gTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGCCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCGCTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATaa                                                                                                               <  2:259635/149‑1 (MQ=255)
                    tACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACTCTCAAAAGGCACACGAAGTCCCAAATACGATGATCGTCCCCATAGTACAACTCGCTGCACTCTGCGTATTAATTCCGCTTACTCCCAGTCAGCTAACACCAAAGATAAt                                                                                                              >  2:419506/1‑149 (MQ=255)
                          agtACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCATCGAAGATAATGCatat                                                                                                        <  2:62529/147‑1 (MQ=255)
                                        tgCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAAAACCAAATACGATGATCGTCGCCATACTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATACATATTGCCAATGCGATTa                                                                                          >  1:92682/1‑149 (MQ=255)
                                                                       caGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCata                                                                                                         >  1:523670/1‑103 (MQ=255)
                                                                       caGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCata                                                                                                         <  2:523670/103‑1 (MQ=255)
                                                                                       cGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCTTTTTTTTACTTATACTCTTACCGCCACTGTCGCCAGCGGAATAGCACTTGTAccc                                           <  1:12891/149‑1 (MQ=255)
                                                                                               gtcgCCAGGCGACTTCTCGCTGCACTCTTCGTATGAATGCCGCTTAGTGCCAGTGGGCTGGCACCGAAGATAATGCATATGGCCCATGCGTTTACCAATACCCCTCCCGCCACGTTCTCCAGCGGAATAGCACTTGTACCCGTCACGtt                                   <  2:114996/149‑1 (MQ=255)
                                                                                                tcgCCATGCTACTGCTCGCGGCACTCTGCGTATTTATTGCGCTTAGTGCCAGTTGGCTGGCATCGAATATAATGCATATTGCCCTTTCGTTTACCAATATTCTTCCCGCTTTTGTCGTCAGCGGAATAGCACTTGTACCCGGCACGttt                                  <  1:110554/149‑1 (MQ=255)
                                                                                                   ccATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCAc                                                   <  1:310487/129‑1 (MQ=255)
                                                                                                   ccATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCAc                                                   >  2:310487/1‑129 (MQ=255)
                                                                                                        cTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAATGCCAGTTGGCTAGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGACGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATAcaca                          >  1:319087/1‑149 (MQ=255)
                                                                                                              cTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGTACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCTCCACTGTCGCCAGCGGAATAGCACTGGTACCCGGCACGTTTCATACACAGGTCAc                    <  1:476560/149‑1 (MQ=255)
                                                                                                                 gcgGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCAcccc                 >  1:580768/1‑149 (MQ=255)
                                                                                                                      aCTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGAAACAAAGATAATGCATATTGCCAATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCGGGCACGTTTCATACACAGGTCACCCCCTCAt            >  1:357802/1‑149 (MQ=255)
                                                                                                                          tGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGACACCGAAAATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACt        >  1:47222/1‑149 (MQ=255)
                                                                                                                           gccgaTTAATTCCGCTTAGGCCCAGGTGCCTGGTACCCAAGACAATGCTTTTCGCCCTTGCTTTCTCCAATACCCCCCCCCCCGCTGCCGCCAGCGGTAGACCACTTTTACCCGGCACGTTTCATACACAGGTCCCCCCCTAACTACTg       <  2:312129/145‑1 (MQ=255)
                                                                                                                                ttAATTGCGCTTAGGGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCt  <  2:508250/149‑1 (MQ=255)
                                                                                                                                        gcTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTa                                              >  2:36484/1‑97 (MQ=255)
                                                                                                                                        gcTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTa                                              <  1:36484/97‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                 
TAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCT  >  NC_000913_3_vch_tpiA/2603148‑2603424

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