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

TGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCG  >  NC_000913_3_vch_pgi/3134721‑3135002                                                                                                                                               |                                                                                                                                            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                                                                                                                                       <  2:230511‑M1/36‑1 (MQ=255)                                 tgggtatttactttaccggcggcaatttttctttctggtgggctgtactggatagcattgcagttgatttaaccttcctgaAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGaa                                                                                                                                       >  1:230511‑M1/82‑117 (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)                                                                                             cagttgatttaaccttcctgaaaatgcccggtcctggcgatcgggcatttTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgtaaagtaa                                                                    <  2:209813‑M2/74‑1 (MQ=255)                                                                                             cagttgatttaaccttcctgaaaatgcccggtcctggcgatcgggcatttTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgtaaagtaa                                                                    >  1:209813‑M2/51‑124 (MQ=255)                                                                                                          aaaaaTGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGCTCCTGGTTCCCGCTTTACGc                              >  2:122656/1‑149 (MQ=255)                                                                                                           ttcctgaaaatgcccggtcctggcgatcgggcatttTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCa                                                             <  2:230511‑M2/81‑1 (MQ=255)                                                                                                           ttcctgaaaatgcccggtcctggcgatcgggcatttTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCa                                                             >  1:230511‑M2/37‑117 (MQ=255)                                                                                                           aaaaTGGAAATGCCCGATGGCAAGTACCTGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCt                             <  2:171783/149‑1 (MQ=255)                                                                                                              aTGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACTCACGCCATCAGGCTGCTGGTTACCGTTTTAAGCTGTa                          >  1:228325/1‑149 (MQ=255)                                                                                                                   AATGCCCGACCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATCATTTCCATAACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCa                     >  1:384022‑M1/1‑149 (MQ=255)                                                                                                                            tCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGGcc            >  1:25336/1‑149 (MQ=255)                                                                                                                               ccAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTccacc         <  2:377437/149‑1 (MQ=255)                                                                                                                               ccAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGCccacc         >  1:109586/1‑149 (MQ=255)                                                                                                                                cAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCTTccacca        >  1:115562/1‑149 (MQ=255)                                                                                                                                      cgggcatttTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCg  >  1:104212‑M2/10‑149 (MQ=255)                                                                                                                                               |                                                                                                                                            TGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCG  >  NC_000913_3_vch_pgi/3134721‑3135002

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