New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_pgi	3134834 =	15 (0.680)	4 (0.190)	4/266	2.5	20.5%	intergenic (+11/+38)	yghT/pitB	putative ATP‑binding protein/phosphate transporter
	?	NC_000913_3_vch_pgi	3134863 =	17 (0.820)					intergenic (+40/+9)	yghT/pitB	putative ATP‑binding protein/phosphate transporter

CATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTA  >  NC_000913_3_vch_pgi/3134694‑3134978                                                                                                                                             |                                                                                                                                                 cgatggtggtggacggcggtgggttacagcgtaaaacggtaaccagcatcctgatggcgtgggtatttactttaccggcggcaatttttctttctggtgggctgtactggatagcattgcagttgatttaaccttcctgaAAATGCCCg                                                                                                                                          >  1:104212‑M1/141‑149 (MQ=255)                 aCCTTGGGTGGTGCACAACTGGGTGAAATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACCAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGc                                                                                                                          <  2:24147/149‑1 (MQ=255)                  ccTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCa                                                                                                                         >  1:242629/1‑149 (MQ=255)                            tgCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTGAAAAATGGACATGCCCGATCGCCAGGACCGGGGATTTTCAGGaa                                                                                                               >  1:171783/1‑149 (MQ=255)                            tgCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGAGGGAGAAAAAGCGATGACGGATGTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGaa                                                                                                               <  2:107751/149‑1 (MQ=255)                                            atCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCa                                                                                               >  1:377437/1‑149 (MQ=255)                                                     aaTCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATcc                                                                                      >  1:58461/1‑149 (MQ=255)                                                     aaTCAGCCACTGGAACAGGTGTTGATTGCGTGATGGATGGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATcc                                                                                      >  1:387919/1‑149 (MQ=255)                                                            tgggtatttactttaccggcggcaatttttctttctggtgggctgtactggatagcattgcagttgatttaaccttcctgaAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGaa                                                                                                               >  1:230511‑M1/82‑117 (MQ=255)                                                            tgggtatttactttaccggcggcaatttttctttctggtgggctgtactggatagcattgcagttgatttaaccttcctgaAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGaa                                                                                                               <  2:230511‑M1/36‑1 (MQ=255)                                                               tGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTCCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGccc                                                                            >  2:346152/1‑149 (MQ=255)                                                                   ttactttaccggcggcaatttttctttctggtgggctgtactggatagcattgcagttgatttaaccttcctgaAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTg                                                                                                 <  2:209813‑M1/50‑1 (MQ=255)                                                                   ttactttaccggcggcaatttttctttctggtgggctgtactggatagcattgcagttgatttaaccttcctgaAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTg                                                                                                 >  1:209813‑M1/75‑124 (MQ=255)                                                                                           aGCGAGTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgtaaa                                                <  2:198597/149‑1 (MQ=255)                                                                                              gATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgtaaagta                                             >  1:47876/1‑149 (MQ=255)                                                                                                       ttttATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCAc                                    <  2:242629/149‑1 (MQ=255)                                                                                                              acctcccgTGATTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCTGGCACTTTCAGAAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCa                             <  1:375354/143‑1 (MQ=255)                                                                                                                                     aaaaaTGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGCTCCTGGTTCCCGCTTTACGc      >  2:122656/1‑149 (MQ=255)                                                                                                                                      aaaaTGGAAATGCCCGATGGCAAGTACCTGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCt     <  2:171783/149‑1 (MQ=255)                                                                                                                                         aTGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACTCACGCCATCAGGCTGCTGGTTACCGTTTTAAGCTGTa  >  1:228325/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                 CATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTA  >  NC_000913_3_vch_pgi/3134694‑3134978

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.