New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_pgi	3375666 =	19 (1.040)	4 (0.230)	4/272	1.9	18.6%	coding (798/1368 nt)	dcuD	putative transporter
	?	NC_000913_3_vch_pgi	3375686 =	17 (0.980)					coding (818/1368 nt)	dcuD	putative transporter
Rejected: Frequency below/above cutoff threshold.

AATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCAGTTTGCCAACGTGGTAACGCTGGTGGTC  >  NC_000913_3_vch_pgi/3375547‑3375831                                                                                                                                            |                                                                                                                                                  aatcaCGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTggtggtggt                                                                                                                                          <  2:111458/149‑1 (MQ=255)            ggCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGa                                                                                                                               <  1:197447/149‑1 (MQ=255)              cAGAGCAAAAAGCTCTCGATAATGTCCCGGCGCTCTATTACGCCATTTTACCTGCGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGt                                                                                                                             <  1:278390/149‑1 (MQ=255)                   aaaaaacgcctgcacatcgtccattgtttcgcgcaagttatgcttgcggaagaactcaacaaacatcgtcacagtcaaactcagtaacatcaccaccaccaGATGCAGTCAGAACTGCATCtgg                                                                                                                                                 <  2:243026‑M1/23‑1 (MQ=255)                   aaaaaacgcctgcacatcgtccattgtttcgcgcaagttatgcttgcggaagaactcaacaaacatcgtcacagtcaaactcagtaacatcaccaccaccaGATGCAGTCAGAACTGCATCtgg                                                                                                                                                 >  1:243026‑M1/102‑124 (MQ=255)                                        ccgccgCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGtt                                                                                                   >  2:186159/1‑149 (MQ=255)                                                      cGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCCGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCGTa                                                                                     >  1:303368/1‑149 (MQ=255)                                                         cATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGAGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTGTGTTGAGTTCTTCCGCAAGCATAACt                                                                                  >  1:210261/1‑149 (MQ=255)                                                             ttATCTGTGATGCCGTTAAGCCGGATGCTCGGCTCGCTGTTCCTCGTCCACGTCGGGCTTATGCAGATAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGCCGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTgcg                                                                              <  2:159430/149‑1 (MQ=255)                                                                  tgtgATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTGTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGaaa                                                                         >  1:200949/1‑149 (MQ=255)                                                                                 ccTGATGCTAGTTTCGATGGACCTCGCCAACGTCGGGCCGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGc                                                          <  1:140703/149‑1 (MQ=255)                                                                                  cTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGcg                                                         >  1:99042/1‑148 (MQ=255)                                                                                  cTGATGCCCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCa                                                         <  1:255360/149‑1 (MQ=255)                                                                                                 actcagtaacatcaccaccaccaGATGCAGTCAGAACTGCATCtggtg                                                                                                                                               <  1:25605‑M1/25‑1 (MQ=255)                                                                                                 actcagtaacatcaccaccaccaGATGCAGTCAGAACTGCATCtggtg                                                                                                                                               >  2:25605‑M1/24‑48 (MQ=255)                                                                                                    ttCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCAt                                       <  1:186159/149‑1 (MQ=255)                                                                                                       cTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATggg                                    >  1:117323/1‑149 (MQ=255)                                                                                                              aCGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCAg                             <  1:304682/149‑1 (MQ=255)                                                                                                                   caccagatgcagttctgactgcatcTGGTGGTGGTGATGTTACTGAGt                                                                                                                             <  1:25605‑M2/23‑1 (MQ=255)                                                                                                                   caccagatgcagttctgactgcatcTGGTGGTGGTGATGTTACTGAGt                                                                                                                             >  2:25605‑M2/26‑48 (MQ=255)                                                                                                                     gCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCATCCTACTATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCAGTTTGCCa                      <  1:144385/149‑1 (MQ=255)                                                                                                                     ccagatgcagttctgactgcatcTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGtttttt                                               >  1:243026‑M2/24‑124 (MQ=255)                                                                                                                     ccagatgcagttctgactgcatcTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGtttttt                                               <  2:243026‑M2/101‑1 (MQ=255)                                                                                                                      cTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGCGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCCGTTTGccca                     >  2:185917/1‑147 (MQ=255)                                                                                                                              GTCAGAACTGCATTTGGTGGTGGTGATGTTACTGAGTTTGGCTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCAGTTTGCCAACGTGGTaa             >  2:72644‑M1/1‑149 (MQ=255)                                                                                                                               TCAGAACTGAATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCAGTTTGCCAACGTGGTAAc            >  1:13159‑M1/1‑149 (MQ=255)                                                                                                                                 aGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGGCGTTTGCCAGCGTGGTAACGc          >  1:267293/1‑149 (MQ=255)                                                                                                                                    actgcatcTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTCAGTTCTTCCGCAAGCATTACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCAGTTTGCCAACGTGGTAACGCtgg       <  1:299042‑M2/141‑1 (MQ=255)                                                                                                                                         atcTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTCAGCAGTTTGCCCACGTGTtaacgctgtgtgtc  >  2:48005‑M2/4‑143 (MQ=255)                                                                                                                                            |                                                                                                                                                  AATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCAGTTTGCCAACGTGGTAACGCTGGTGGTC  >  NC_000913_3_vch_pgi/3375547‑3375831

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.