New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_pgi	306122 =	22 (1.000)	3 (0.150)	3/266	2.8	13.6%	coding (696/1644 nt)	ecpD	polymerized tip adhesin of ECP fibers
	?	NC_000913_3_vch_pgi	306162 =	18 (0.890)					coding (656/1644 nt)	ecpD	polymerized tip adhesin of ECP fibers
Rejected: Frequency below/above cutoff threshold.

TCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACA  >  NC_000913_3_vch_pgi/305974‑306260                                                                                                                                                     |                                                                                                                                           tgcccgtcggcagcagctttacaatgaccatcaatgcctgccagacctctgtgaactatgacgccagcagcggcgcacgctgtaaggatcaggcctccggcaactggtatgttcgcaacgtcacccatacgaaagcagcaaatctacgg                                                                                                                                            >  1:289767‑M1/149‑149 (MQ=255)         ggcagcagctttacaatgaccatcaatgcctgccagacctccgttagctgtgacgccagcagcggcgcacgctgtaaggatcaggcctccggcaactggtatgttcgcaacgtcacccatacgaaagcagcaaatctacgGTTGATaaa                                                                                                                                    >  1:168896‑M1/141‑149 (MQ=255)           cagcagctttacaatgaccatcaatgcctgccagacctctgtgaactatgacgccagcagcggcgcacgctgtaaggatcaggcctccggcaactggtatgttcgcaacgtcacccatacgaaagcagcaaatctacgGTTGATAAATa                                                                                                                                  <  2:196436‑M1/11‑1 (MQ=255)           cagcagctttacaatgaccatcaatgcctgccagacctctgtgaactatgacgccagcagcggcgcacgctgtaaggatcaggcctccggcaactggtatgttcgcaacgtcacccatacgaaagcagcaaatctacgGTTGATAAATa                                                                                                                                  >  1:233462‑M1/139‑149 (MQ=255)              atGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACTACTTAATCGTGAATGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACtt                                                                                                                               >  2:44770/1‑149 (MQ=255)                                   tCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGCGTGGGTATTTa                                                                                                          >  2:50348/1‑149 (MQ=255)                                      gtttgtttGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGTGTATTTATCa                                                                                                       >  2:254521/1‑149 (MQ=255)                                      gtttgtttGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCa                                                                                                       >  1:314625/1‑149 (MQ=255)                                       tgccagacctctgtgaactatgacgccagcagcggcgcacgctgtaaggatcaggcctccggcaactggtatgttcgcaacgtcacccatacgaaagcagcaaatctacgGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCaa                                                                                                      >  1:145856‑M1/111‑149 (MQ=255)                                                aGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGAtt                                                                                             >  2:372683/1‑149 (MQ=255)                                                    tataCTTAACCATCTTACAACTACATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTAATCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGAGTATTTATCAACCGTAGATTtgct                                                                                         <  1:328086/149‑1 (MQ=255)                                                              cATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCCGCGAGTGGGTATTTATCAACCGTAGATTTCCTGCTTTCGTAt                                                                               >  2:93041/1‑149 (MQ=255)                                                              cATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTAt                                                                               >  1:224349/1‑149 (MQ=255)                                                                  ttACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGGATTTATCAACCGTAGGTTGGCTGCTTTCGTATggct                                                                           >  2:38635/1‑147 (MQ=255)                                                                      aaCTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATTAACCGTAGATTTGCTGCTTTCGTTTGGGTGTCg                                                                       >  2:205358/1‑149 (MQ=255)                                                                                  aagggCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCGTGCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATAc                                                           <  2:315321/146‑1 (MQ=255)                                                                                   taaggatcaggcctccggcaactggtatgttcgcaacgtcacccatacgaaagcagcaaatctacgGTTGATAAATACTTCCGCCAGCGAGTGGGTAtt                                                                                                            >  2:148506‑M1/67‑99 (MQ=255)                                                                                   taaggatcaggcctccggcaactggtatgttcgcaacgtcacccatacgaaagcagcaaatctacgGTTGATAAATACTTCCGCCAGCGAGTGGGTAtt                                                                                                            <  1:148506‑M1/33‑1 (MQ=255)                                                                                    cGGCTGCCGATGGTGTGAGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTATCGTCAAGGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCa                                                         <  1:170634/149‑1 (MQ=255)                                                                                          ccGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGCCGTTGCGAACATACCAGTTGcc                                                   >  2:115527/1‑149 (MQ=255)                                                                                            gATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATAACAGTTGccag                                                 >  1:96619/1‑147 (MQ=255)                                                                                                      gTCCGGGAGTGGGCGTTCCCTTCGCCCAGGGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATc                                       <  2:337146/149‑1 (MQ=255)                                                                                                      gTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGATGACGGTGCGCACATACCAGTTGCCCGAGGCCTGATc                                       >  1:332543/1‑149 (MQ=255)                                                                                                         cGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGATACTCCGTCGCTGTTGCTAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGCTTTGCTGGTTCCCTATGGTTGTGGTTGCGAACATACCAGTTCCCCGAGGCCTCATCCta                                    >  2:275410/1‑148 (MQ=255)                                                                                                            cAGTCGGCGTTCCCTTCGCCCAAAGTCGATACACCGTCGCTGTTGTTAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTCGCCGGAGGCCTGATCCttcaa                                 >  1:307311/1‑146 (MQ=255)                                                                                                                       cccTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGcgccgc                      <  1:78476/149‑1 (MQ=255)                                                                                                                                       ggTCCTCCGTCGCTGTTGATAGATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGtt      <  2:247773/149‑1 (MQ=255)                                                                                                                                           ctgctgcGCTGTTGATAACTACTTCCGCCAGCGAGTGTGTATTTATCAACCGTAGATTTGCCGCTGTCGTATGGGTGACGTTGCGAACATACCAGTTGGCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTcaca  <  2:55763/143‑1 (MQ=255)                                                                                                                                           cTCCGGCGCTGTCGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTcaca  <  2:268832/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                           TCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACA  >  NC_000913_3_vch_pgi/305974‑306260

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