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.

GCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATTAATGGAAATGCTAATACTACGGCGAACAATGCTCGACTCACAGGGAACTCCTTAACGTTAT  >  NC_000913_3_vch_pgi/1841868‑1842155
                                                                                                                                                |                                                                                                                                               
gCTTTACGCGCCGCCAGACCAGCGGCTTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTGCCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGAc                                                                                                                                             <  1:318664/149‑1 (MQ=255)
       gcgcCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCCTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGc                                                                                                                                      >  1:50910/1‑149 (MQ=255)
                             gTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCgg                                                                                                                >  1:101868/1‑149 (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)
                                                                                                                         gcagcaCTGAGTACAGGGCTGCGTttatcaggatcaaaattattctgaaggggcggtaattaaagccgaaggaatactgctgcaatgccagcgtgatgacaaaacccttagcactaatccgctggtctggcggcgcgtaaagccataaa                    <  1:291434‑M2/149‑126 (MQ=255)
                                                                                                                           agcaCTGAGGACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTGTACCTCCACGTCCGGACGATAATGCGGTTTAGCCAGTGCAATTAATGGAAATGCTAATACTACGGCGAACAATGCTCGACTCACAggg                  <  1:89340/149‑1 (MQ=255)
                                                                                                                                 gAGTACAGGGCTGCGCTGACTGCGCGCCTGAAGTGAACACTTATGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATTAATGGAAATGCTAAGACTACGGCGAACAATGCTCGACTCACGGGGAACTcg            >  1:256214/1‑148 (MQ=255)
                                                                                                                                          gCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATTAATGGAAATGCTAATACTACGGCGAACAATGCTCGACTCACAGGGAACTCCTTCACGTTa   <  1:146396/149‑1 (MQ=255)
                                                                                                                                          gCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCACTGCGTGCCGGTGATATGCTAATACTACGGCGAACAATGCTCGACTCACAGGGAACTCCTTAACGTTa   >  1:74879/1‑149 (MQ=255)
                                                                                                                                           cTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATTAATGGAAATGCTAATACTACGGCGAACAATGCTCGACTCACCGGGAACTCCTTAACGTTAt  >  1:1351/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGCCAGTGCAATTAATGGAAATGCTAATACTACGGCGAACAATGCTCGACTCACAGGGAACTCCTTAACGTTAT  >  NC_000913_3_vch_pgi/1841868‑1842155

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