New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_pgi = 187222923 (1.040)3 (0.140) 3/270 2.9 12.0% coding (189/1491 nt) yeaJ putative diguanylate cyclase
?NC_000913_3_vch_pgi = 1872232 22 (1.040)coding (192/1491 nt) yeaJ putative diguanylate cyclase
Rejected: Frequency below/above cutoff threshold.

TGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGC  >  NC_000913_3_vch_pgi/1872085‑1872364
                                                                                                                                                    |                                                                                                                                    
tgtCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGa                                                                                                                                     >  2:369797‑M2/1‑148 (MQ=255)
     ctctgctcACATCTTCCTTCCATATCTTTCAACTTGCCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGatgcgc                                                                                                                                >  1:85459‑M2/4‑143 (MQ=255)
           tGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGAATCCTCCTGTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCtctc                                                                                                                          >  1:348648/1‑149 (MQ=255)
                           tATGTTTGGAGTTGTCGCCAGTGACAGGGCAATGAGTGCCTGTCTGCGCTATATCGTGCTGAAAGCAGTCTGCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGTGCTCTCGCTGCTGAGCTGTCgg                                                                                                          <  2:26441/149‑1 (MQ=255)
                             tttttGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCCCTCGCTGCTGAGCAGTCGGaa                                                                                                        >  2:2369/1‑149 (MQ=255)
                              tttgcaaCTTGTCGCGAGCGAAAGAGCAATGAGTGCCGAACTGCGCGAAATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGACTCCGAGGATTGCCGCGCATGGGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAg                                                                                                       <  1:175810/144‑1 (MQ=255)
                                   aaCTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCg                                                                                                  >  1:199629/1‑149 (MQ=255)
                                                  aCAGCGCAATGAGTGCCTTTTTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGgcgc                                                                                   <  1:162500/149‑1 (MQ=255)
                                                    agagCAATGAGTGCCTAGCTGACGCATATCGTTCAGAAAGCAGACTCCTCCTTTCTTTATGAGAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGgcgcgc                                                                                 <  2:63027/149‑1 (MQ=255)
                                                                             tataTGGCGCAGAAAGCAGACTTCTCTTTTCTTTATGATAAGTATCAGAATCAGAGTATCGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGGCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGc                                                         <  2:131898/149‑1 (MQ=255)
                                                                                   gTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGt                                                             >  1:94151/1‑139 (MQ=255)
                                                                                   gTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGt                                                             <  2:94151/139‑1 (MQ=255)
                                                                                    tGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCTGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATagc                                                  >  2:87613/1‑147 (MQ=255)
                                                                                                    ctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTAAGTCTGCCAATAACACCCATGGCTCAAAcc                                  >  2:57565/1‑149 (MQ=255)
                                                                                                            tatATGATAAGTATCAGAATCAGACTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGcc                          <  1:170986/147‑1 (MQ=255)
                                                                                                                tGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGatgcgcggcaatactctgattctgatacttatcataaagaaaggaggagtctgctttctgcccgagatagcgcagataggcactcattgctctgtcgctggcgacaagttca                      <  1:227417‑M2/149‑113 (MQ=255)
                                                                                                                 gATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATaa                     >  1:325215/1‑149 (MQ=255)
                                                                                                                  aTAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATaaa                    >  2:372685/1‑149 (MQ=255)
                                                                                                                    aaGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAacac                                                >  2:201742/1‑119 (MQ=255)
                                                                                                                    aaGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAacac                                                <  1:201742/119‑1 (MQ=255)
                                                                                                                       tATCAGAATCAGAGTATTGCCGCGCATGTGatgcgcggcaatactctgattctgatacttatcataaagaaaggaggagtctgctttctgcacgatatagcgcagataggcactcattgctctgtcgctgg                                 >  1:6325‑M2/1‑30 (MQ=255)
                                                                                                                       tATCAGAATCAGAGTATTGCCGCGCATGTGatgcgcggcaatactctgattctgatacttatcataaagaaaggaggagtctgctttctgcacgatatagcgcagataggcactcattgctctgtcgctgg                                 <  2:6325‑M2/131‑102 (MQ=255)
                                                                                                                       tATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATAccc               <  2:348648/149‑1 (MQ=255)
                                                                                                                         tcagaatcagaGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCgg             <  1:168562/149‑1 (MQ=255)
                                                                                                                         tcagaatcagaGTATGGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCgg             <  1:169193/149‑1 (MQ=255)
                                                                                                                                  gagTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTAcgcg    <  1:386698/149‑1 (MQ=255)
                                                                                                                                    gTAGTGCCGGGCATGGGATGGGCGCTGGCGCTGCTGAGAAGTCGGAGGTGTCGCGAGAGCAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCAGAAATACCCGGGCTTACGCGGc  <  2:181335/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                    
TGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGC  >  NC_000913_3_vch_pgi/1872085‑1872364

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.