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.

GCGTACGCGCCATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGATCTTC  >  NC_000913_3_vch_pgi/2958359‑2958653                                                                                                                                                    |                                                                                                                                                    gcgTACGCGCCATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTAGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCACGCAAGGCGTCGTTCTCAACTTCCAGATACAAAGCCGTa                                                                                                                                                    >  2:71886‑M1/1‑148 (MQ=255)    tACGCGCCATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGa                                                                                                                                                 <  1:159545/149‑1 (MQ=255)       gcgcCATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATac                                                                                                                                              <  1:169178‑M2/149‑3 (MQ=255)           cATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGgcg                                                                                                                                          >  1:139916/1‑149 (MQ=255)           cATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGgcg                                                                                                                                          >  1:139900/1‑149 (MQ=255)                                cGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTAc                                                                                                                     >  2:73863/1‑149 (MQ=255)                                         cTCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCa                                                                                                            >  2:140978/1‑149 (MQ=255)                                          tCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCAt                                                                                                           >  2:276449/1‑149 (MQ=255)                                          tCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCAt                                                                                                           >  2:147720/1‑149 (MQ=255)                                           cACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAAGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACCGCGCAGTGCTGGCATTGTGACTACCGCGGTCCAtt                                                                                                          >  2:68941/1‑149 (MQ=255)                                               gtCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTga                                                                                                      <  2:107718/148‑1 (MQ=255)                                                    cATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTatcgcacggctttctatctggaagttgag                                                                                                                        >  1:351340‑M1/1‑97 (MQ=255)                                                    cATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTatcgcacggctttctatctggaagttgag                                                                                                                        <  2:351340‑M1/126‑30 (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)                                                                                                                                                   gTatcgcacggctttctatctggaagttgagaacgacgccttgcctggtgcggtagaccgcctggccgatgctattgctgaacctttgctcgacaagaaatatgccgaacgtgagcgtaatgcggtgaacgctgaattaaccatggcgc  <  1:169178‑M1/149‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCGTACGCGCCATGGTTAATTCAGCGTTCACCGCATTACGCTCACGTTCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGATCTTC  >  NC_000913_3_vch_pgi/2958359‑2958653

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.