New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_pgi	= 1872229	23 (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.

AGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTG‑CGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCG  >  NC_000913_3_vch_pgi/1872083‑1872377                                                                                                                                                    |                                                                                                                                                     aGTGTCATTGTCGTGACATCTTCCTTCCTCATTTTTGACCTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTG‑CGCTATATCGTTCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCAtgt                                                                                                                                                    <  1:197962/149‑1 (MQ=255)   tgtCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTG‑CGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGa                                                                                                                                                  >  2:369797‑M2/1‑148 (MQ=255)        ctctgctcACATCTTCCTTCCATATCTTTCAACTTGCCGCCAGCGACAGAGCAATGAGTGCCTATCTG‑CGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGatgcgc                                                                                                                                             >  1:85459‑M2/4‑143 (MQ=255)              tGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTG‑CGCTATATCGTGCAGAAAGCAGAATCCTCCTGTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCtctc                                                                                                                                       >  1:348648/1‑149 (MQ=255)                              tATGTTTGGAGTTGTCGCCAGTGACAGGGCAATGAGTGCCTGTCTG‑CGCTATATCGTGCTGAAAGCAGTCTGCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGTGCTCTCGCTGCTGAGCTGTCgg                                                                                                                       <  2:26441/149‑1 (MQ=255)                                tttttGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTG‑CGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCCCTCGCTGCTGAGCAGTCGGaa                                                                                                                     >  2:2369/1‑149 (MQ=255)                                 tttgcaaCTTGTCGCGAGCGAAAGAGCAATGAGTGCCGAACTG‑CGCGAAATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGACTCCGAGGATTGCCGCGCATGGGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAg                                                                                                                    <  1:175810/144‑1 (MQ=255)                                    tGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTG‑CGCTATCTCGGGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATcacatgcgcggcaatactctgattctgatacttatca                                                                                                                 <  1:227417‑M1/149‑38 (MQ=255)                                      aaCTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTG‑CGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCg                                                                                                               >  1:199629/1‑149 (MQ=255)                                               ccAGCGACAGAGCAATGAGTGCCTATCTG‑CGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATcacatgcgcggcaatactctgattctgata                                                                                                                        <  2:6325‑M1/131‑31 (MQ=255)                                               ccAGCGACAGAGCAATGAGTGCCTATCTG‑CGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATcacatgcgcggcaatactctgattctgata                                                                                                                        >  1:6325‑M1/1‑101 (MQ=255)                                                     aCAGCGCAATGAGTGCCTTTTTG‑CGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGgcgc                                                                                                <  1:162500/149‑1 (MQ=255)                                                       agagCAATGAGTGCCTAGCTGACGC‑ATATCGTTCAGAAAGCAGACTCCTCCTTTCTTTATGAGAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGgcgcgc                                                                                              <  2:63027/149‑1 (MQ=255)                                                                                tataTGGCGCAGAAAGCAGACTTCTCTTTTCTTTATGATAAGTATCAGAATCAGAGTATCGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGGCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGc                                                                      <  2:131898/149‑1 (MQ=255)                                                                                      gTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGt                                                                          <  2:94151/139‑1 (MQ=255)                                                                                      gTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGt                                                                          >  1:94151/1‑139 (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)                                                                                                                                               gcgcATcacatgcgcggcaatactctgattctgatacttatcataaagaaaggaggagtctgctttctgcacgatatagcgcagataggcactcattgctctgtcgctggcggcaagttgaaagatatggaaggaagatgtgagcagag       >  1:85459‑M1/1‑6 (MQ=255)                                                                                                                                                    tcacatgcgcggcaatactctgattctgatacttatcataaagaaaggaggagtctgctttctgcacgatatagcgcagataggcactcattgctctgtcgctggcgacaagttcaaaaataaggaaggaagatgtcagcacaatgaca  >  2:369797‑M1/1‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     AGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTG‑CGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCG  >  NC_000913_3_vch_pgi/1872083‑1872377

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.