New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_pgi	= 2729634	15 (0.680)	3 (0.140)	3/274	2.9	19.7%	noncoding (1524/1542 nt)	rrsG	16S ribosomal RNA of rrnG operon
	?	NC_000913_3_vch_pgi	4209658 =	10 (0.470)					noncoding (1512/1542 nt)	rrsE	16S ribosomal RNA of rrnE operon
Rejected: Frequency below/above cutoff threshold.

TGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACG  >  NC_000913_3_vch_pgi/2729489‑2729780                                                                                                                                                  |                                                                                                                                                   tGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAg                                                                                                                                                 <  2:324792/149‑1 (MQ=255)                 gAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTAc                                                                                                                                 <  1:104048/149‑1 (MQ=39)                  aaaGCGAACGGCCTTATTCTCTTCACCCTCACTCCCAACGTGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGtaggggaacctgcggttgga                                                                                                                                <  2:119027‑M1/149‑21 (MQ=255)                        aatggCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTa                                                                                                                          <  2:359899/146‑1 (MQ=255)                        aaCGGGCTTACTCTCTTCAGCCTCACTCCCAACTCGTAAACGCCTTGCGATTTCCTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTTATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTa                                                                                                                          <  2:12024/149‑1 (MQ=14)                        aaCGGCTGTATGCTCTTCAGCCTCGCTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGGGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTa                                                                                                                          <  2:325219/149‑1 (MQ=14)                        aaCGGCCTCATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAATGACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTa                                                                                                                          <  2:325104/149‑1 (MQ=255)                         aCGGCCTTATTCCCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTAc                                                                                                                         <  1:61767/149‑1 (MQ=255)                                       ttCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGATGATCCAACCGCAGGTTCCGCTCCGGTTACCTTGTTACGACTTCCCCCCAGt                                                                                                           >  2:315131/1‑149 (MQ=255)                                            ccTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGa                                                                                                      <  2:353143/149‑1 (MQ=255)                                                aCTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATca                                                                                                  <  1:148430/149‑1 (MQ=255)                                                          cgTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCCGTCATGAATCACAAAGGGGTa                                                                                        >  2:124057/1‑149 (MQ=255)                                                                        tATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGaa                                                                          <  2:331628/149‑1 (MQ=35)                                                                          ttCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAgg                                                                        >  2:11474/1‑149 (MQ=33)                                                                          gtcccTTTTCATCGGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAgg                                                                        <  1:128597/145‑1 (MQ=11)                                                                             aCTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTa                                                                     >  2:194995/1‑149 (MQ=31)                                                                                                                             aggtaaggAGGTGATCCAACCGtaggggaacctgcggt                                                                                                                                    <  2:273245‑M1/38‑17 (MQ=255)                                                                                                                             aggtaaggAGGTGATCCAACCGtaggggaacctgcggt                                                                                                                                    >  1:273245‑M1/1‑22 (MQ=255)                                                                                                                                                ccGtagggtaccctgcggttggatcacctccttaccttaaagaagcgttctttgcagtgctcacacagattgtctgatagaaagtgaaaagcaaggcgtcttgcgaagcagactgatacgtccccttcgtctagaggcccaggacaccg  <  1:99785‑M1/149‑147 (MQ=255)                                                                                                                                                 cGtaggggaacctgcggttggatcacctccttaccttaaagaagcgttctttgcagtgctcacacagattgtctgatagaaagtgaaaagcaaggcgtcttgcgaagcagactgatacgtccccttcgtctagaggcccag         >  2:25703‑M1/1‑2 (MQ=255)                                                                                                                                                 cGtaggggaacctgcggttggatcacctccttaccttaaagaagcgttctttgcagtgctcacacagattgtctgatagaaagtgaaaagcaaggcgtcttgcgaagcagactgatacgtccccttcgtctagaggcccag         <  1:25703‑M1/141‑140 (MQ=255)                                                                                                                                                  |                                                                                                                                                   TGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACG  >  NC_000913_3_vch_pgi/2729489‑2729780

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.