New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_pgi = 310653920 (1.090)3 (0.170) 3/274 2.3 13.3% coding (826/1257 nt) nupG nucleoside transporter
?NC_000913_3_vch_pgi = 3106551 20 (1.150)coding (838/1257 nt) nupG nucleoside transporter
Rejected: Frequency below/above cutoff threshold.

TGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGCGCATTCGACTTCTT  >  NC_000913_3_vch_pgi/3106417‑3106687
                                                                                                                          |                                                                                                                                                    
tgaaaaGCTTAGAGTAAGAACCGATGTTGGCCAGCAGCTTTATCTTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTTTTCATTATGACCATGCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAatg                                                                                                                            <  2:165084/144‑1 (MQ=255)
  cacaTCTTCGACAAAGATCCGATGTTTGCCAGCAGCTCTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTCCATTCTGACCATCCCGTTCTTCTTAACCCGCTACGGTATTAAAAACGTaataat                                                                                                                          >  2:246972/1‑149 (MQ=255)
           gACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTCCGGTATTAAGAACGTAGTGATGATCAGtct                                                                                                                 >  2:180987/1‑147 (MQ=255)
            aCAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCCGTAtt                                                                                                                >  1:96893/1‑149 (MQ=255)
                       aTGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGAAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGAt                                                                                                     >  1:306327/1‑149 (MQ=255)
                         gTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATcc                                                                                                   <  2:245231/149‑1 (MQ=255)
                              ccAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGt                                                                                              >  2:215964/1‑149 (MQ=255)
                                agcagcTTTATTGTGGAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGttt                                                                                            <  2:71897/149‑1 (MQ=255)
                                  cagcTTGATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGCATTGTGGCGTGGATCCTGCGTTTTg                                                                                          <  1:293777/149‑1 (MQ=255)
                                     cTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCTCTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTgcgc                                                                                       <  1:190737/149‑1 (MQ=255)
                                     cTTTATTGGGCGCGATGCGTAAATCATTGTGTCGATTTCGCAGATTTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTAttaagaagaacgggatggtcagaatgaacagggtttcagagatctgcgaaa                                                                                       >  2:128685‑M2/1‑98 (MQ=255)
                                      tttATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATtgtg                                                                                                             <  1:58115/126‑1 (MQ=255)
                                      tttATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATtgtg                                                                                                             >  2:58115/1‑126 (MQ=255)
                                        tATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAtaccgtagcggcttaagaagaacg                                                                                                                              <  1:92618‑M1/107‑25 (MQ=255)
                                        tATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAtaccgtagcggcttaagaagaacg                                                                                                                              >  2:92618‑M1/1‑83 (MQ=255)
                                         aTTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGtt                                                                                   <  1:296021/149‑1 (MQ=255)
                                                                catGTCGATTTTGCAGAGCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGt                                                            <  2:53959/149‑1 (MQ=255)
                                                                        tttCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAtaccgtagcggcttaagaagaacgggatggtcagaatgaacagggtttcagaaatctgcgaaatcgacacaatgatttacgcatcgcgcccaataaag                                                    >  2:128685‑M1/1‑51 (MQ=255)
                                                                                aTCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGCTGATCAGTATTGTGGCGTGGATCCTGCGTTTTTCGCTGTTTGCTTACGGCGACCCGCCCCCGTTCGGTACTGTActgct                                            >  1:4207/1‑149 (MQ=255)
                                                                                            cTGTTCATTCTGATCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCgat                                <  1:37555/149‑1 (MQ=255)
                                                                                                tCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATc                            <  1:167471/149‑1 (MQ=255)
                                                                                                   ttCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGtt                         >  1:189349/1‑149 (MQ=255)
                                                                                                             cccGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTgcg               >  1:229186/1‑149 (MQ=255)
                                                                                                               cGTTCTTCTTAAGCCGCTACGGTAttaagaagaacgggatggtcagaatgaacagggtttcagagatctgcgaaatcgacatgatgattgacgcatgctgcacaata                                                       <  1:92618‑M2/107‑84 (MQ=255)
                                                                                                               cGTTCTTCTTAAGCCGCTACGGTAttaagaagaacgggatggtcagaatgaacagggtttcagagatctgcgaaatcgacatgatgattgacgcatgctgcacaata                                                       >  2:92618‑M2/1‑24 (MQ=255)
                                                                                                                          aGCCGCTACGGTATTAAGAAGGTAATGCCGATCAGTATTGTGGCGTGGCTCCTGCGTTTTGCGCTGTTTGCTTACGACGACCCGACTCCGTTCGGTACTGTACGTCTGGTACGGTCGATGCTCCTTTACCGATGCGCACTCCActtctt  >  2:32615/1‑149 (MQ=255)
                                                                                                                          |                                                                                                                                                    
TGCACAGCTTCGACAAAGATCCGATGTTTGCCAGCAGCTTTATTGTGCAGCATGCGTCAATCATCATGTCGATTTCGCAGATCTCTGAAACCCTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGCGCATTCGACTTCTT  >  NC_000913_3_vch_pgi/3106417‑3106687

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.