New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_pgi	= 2958506	28 (1.260)	3 (0.140)	3/274	2.9	9.5%	coding (379/2889 nt)	ptrA	protease III
	?	NC_000913_3_vch_pgi	= 2958511	30 (1.400)					coding (374/2889 nt)	ptrA	protease III
Rejected: Frequency below/above cutoff threshold.

GCGCCATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGATCTTCCAGCGA  >  NC_000913_3_vch_pgi/2958365‑2958659                                                                                                                                                   |                                                                                                                                                     gcgcCATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATac                                                                                                                                                    <  1:169178‑M2/149‑3 (MQ=255)     cATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGgcg                                                                                                                                                >  1:139900/1‑149 (MQ=255)     cATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGgcg                                                                                                                                                >  1:139916/1‑149 (MQ=255)                          cGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTAc                                                                                                                           >  2:73863/1‑149 (MQ=255)                                   cTCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCa                                                                                                                  >  2:140978/1‑149 (MQ=255)                                    tCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCAt                                                                                                                 >  2:147720/1‑149 (MQ=255)                                    tCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCAt                                                                                                                 >  2:276449/1‑149 (MQ=255)                                     cACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAAGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACCGCGCAGTGCTGGCATTGTGACTACCGCGGTCCAtt                                                                                                                >  2:68941/1‑149 (MQ=255)                                         gtCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTga                                                                                                            <  2:107718/148‑1 (MQ=255)                                                  tttCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGCCATTGTGACTACCGCCGTGCATTTTGAGATATTCgg                                                                                                   >  2:47635/1‑149 (MQ=255)                                                       tggggggCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGa                                                                                              <  2:294562/143‑1 (MQ=255)                                                                         cAATAGCATCGGCCAGGCGGTCTACGGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGACGCAGTGCTAGCATTGTGACTACCGCCGTGCATTTTGAGATCTTGGGCCAGACTGTCAGCCTGCGGGTag                                                                            >  2:144585/1‑148 (MQ=255)                                                                              gCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACttttt                                                                       <  1:60891/149‑1 (MQ=255)                                                                               cATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTc                                                                      <  1:73863/149‑1 (MQ=255)                                                                                aTCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCCCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGACCCGACTGTCAGCCTGCGGGTACTCTTTcc                                                                     >  2:343581/1‑148 (MQ=255)                                                                                            gTCTACCGCACCAAGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCATCCTGCGGGTACTTTTTCGACACCATCAGAc                                                         >  2:180125/1‑149 (MQ=255)                                                                                                 ccGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCgg                                                                                >  1:4701/1‑121 (MQ=255)                                                                                                 ccGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCgg                                                                                <  2:4701/121‑1 (MQ=255)                                                                                                        aGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCaa                                             <  1:276449/149‑1 (MQ=255)                                                                                                        aGGCAAGGCGGCGGTCTCAAGTTCCGGAGAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTccg                                             <  2:167978/149‑3 (MQ=255)                                                                                                                     ttCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAg                                >  1:151076/1‑149 (MQ=255)                                                                                                                       ctcAACTTCCAGATAGAAAGCCGTGCGATacggctttctatctggaagttgagaacgacgccttgcctggtgcggtagaccgcctggccgatgctattgctgaacctttgctcgacaagaaatatg                                                     >  1:351340‑M2/1‑29 (MQ=255)                                                                                                                       ctcAACTTCCAGATAGAAAGCCGTGCGATacggctttctatctggaagttgagaacgacgccttgcctggtgcggtagaccgcctggccgatgctattgctgaacctttgctcgacaagaaatatg                                                     <  2:351340‑M2/126‑98 (MQ=255)                                                                                                                       ctcAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGcc                              <  1:147720/149‑1 (MQ=255)                                                                                                                        tcAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTaa                                        <  2:172822/138‑1 (MQ=255)                                                                                                                        tcAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTaa                                        >  1:172822/1‑138 (MQ=255)                                                                                                                        tcAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGccc                             <  2:151076/149‑1 (MQ=255)                                                                                                                                 aGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCGGGCTGTCAGCCTGCGGGTACTTTTCCAGCGTCGGCATCCTCACATGTTAAAGGTAATTTGCCAGCCCCTGGTACGc                    >  1:89780/1‑149 (MQ=255)                                                                                                                                      gAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTTTCAGCCTGCGGGTACTTTTTCGACCCCATAAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCgg               >  1:33471/1‑149 (MQ=255)                                                                                                                                        aaGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGa             <  1:140978/149‑1 (MQ=255)                                                                                                                                               gggATACGGGCGAGGGCTGGCCGTGTGACTACTGTCGTGTATTTTGAGTTATTCGGCCAGCCTGTCAGCCTGCGGGTACTTTATCGTCCCCCTCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTTCTACGCCTCGGGATCTTCCa      <  2:16023/147‑1 (MQ=255)                                                                                                                                                  aTACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGATCTTCCAGCg   >  1:298324/1‑149 (MQ=255)                                                                                                                                                   tacggctttgtatctggaagttgagaacgacgccttgcgtggtgcggtagaccgcctggccgatgctattgctgaacctttgctctacaagaaatatgccgaacgtgagcgtaatgcggtgaacgctgaattaaccatggcgcgtacgc  >  2:71886‑M2/1‑1 (MQ=255)                                                                                                                                                   tACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGATCTTCCAGCGa  >  1:127465/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCGCCATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGATCTTCCAGCGA  >  NC_000913_3_vch_pgi/2958365‑2958659

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.