New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_pgi = 184198619 (1.040)3 (0.160) 3/286 2.4 13.1% coding (150/273 nt) ynjH DUF1496 family protein
?NC_000913_3_vch_pgi = 1842012 21 (1.150)coding (124/273 nt) ynjH DUF1496 family protein
Rejected: Frequency below/above cutoff threshold.

TGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATTAATGGAAATGCTAATACTACGGCGAA  >  NC_000913_3_vch_pgi/1841840‑1842120
                                                                                                                                                  |                                                                                                                                      
tgACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAAc                                                                                                                                      >  2:109879/1‑149 (MQ=255)
 gACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACg                                                                                                                                     <  1:109879/149‑1 (MQ=255)
    aTTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgcag                                                                                                                                  <  1:48339/149‑1 (MQ=255)
                      tttATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAacgcagccctgtactcagtgctgc                                                                                                                <  1:291434‑M1/149‑25 (MQ=255)
                            gCTTTACGCGCCGCCAGACCAGCGGCTTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTGCCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGAc                                                                                                          <  1:318664/149‑1 (MQ=255)
                                   gcgcCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCCTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGc                                                                                                   >  1:50910/1‑149 (MQ=255)
                                                         gTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCgg                                                                             >  1:101868/1‑149 (MQ=255)
                                                                        tcatcaCGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAgcgcagccctgtactcagtgctgcgtttatcagga                                                                                                     >  1:231474‑M1/1‑75 (MQ=255)
                                                                        tcatcaCGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAgcgcagccctgtactcagtgctgcgtttatcagga                                                                                                     <  2:231474‑M1/110‑36 (MQ=255)
                                                                                     aTTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGAGGATaa                                                 >  1:135519/1‑149 (MQ=255)
                                                                                           gcagTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCggg                                           <  1:263381/149‑1 (MQ=255)
                                                                                                    cTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTCCGTTTCCCACCACGTCCGGCCGATAATGCGGGTTAGCCAGt                                  >  2:178499/1‑149 (MQ=255)
                                                                                                     ttCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTg                                 <  1:314400/149‑1 (MQ=255)
                                                                                                               aTTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATTAATgg                       >  1:136898/1‑149 (MQ=255)
                                                                                                                 tACCGCCCCTTCAGAATAATTTTGATCCAGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATTAATGGaa                     >  1:111070/1‑149 (MQ=255)
                                                                                                                  aCCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATTAATGGaaa                    >  1:166536/1‑149 (MQ=255)
                                                                                                                          ttCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATTAATGGAAATGCTAAta            >  2:132756/1‑149 (MQ=255)
                                                                                                                          ttCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATCAATGGAAATGCTAAta            >  1:294733/1‑149 (MQ=255)
                                                                                                                             aGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTCCCCGCCTGAACTTAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATTAATGGCAATGCTAAtacta         >  2:118012/1‑149 (MQ=255)
                                                                                                                                ataatTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATTAATGGAAATGCTAATACTACgg      >  1:41467/1‑149 (MQ=255)
                                                                                                                                    ttttGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATTAATGGAAATGCTAATACTACGGCGaa  <  2:136898/149‑1 (MQ=255)
                                                                                                                                    ttttGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATTAATGGAAATGCTAATACTACGGCGaa  >  2:188597/1‑149 (MQ=255)
                                                                                                                                          tCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCttatcaggatcaaaattattctgaaggggcggtaattaaagccgaaggaatactgctgcaatgccagcgtgatga                                   >  1:231474‑M2/1‑35 (MQ=255)
                                                                                                                                          tCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCttatcaggatcaaaattattctgaaggggcggtaattaaagccgaaggaatactgctgcaatgccagcgtgatga                                   <  2:231474‑M2/110‑76 (MQ=255)
                                                                                                                                                  |                                                                                                                                      
TGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATTAATGGAAATGCTAATACTACGGCGAA  >  NC_000913_3_vch_pgi/1841840‑1842120

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