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.

CCGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCG  >  NC_000913_3_vch_pgi/4209529‑4209803                                                                                                                                  |                                                                                                                                                  aaagcgaacggccttattctcttcaccctcactcccaacgtgtaaacgccttgctattcacttttcatcagacaatctgtgtgagcactgcaaagtacgcttctttaaggtaaggaggtgatccaaccgTAGGGGAACCTGCGGTTGGa                                                                                                                                <  2:119027‑M2/20‑1 (MQ=255)                                                                   cTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGaaa                                                              >  2:355944/1‑149 (MQ=32)                                                                        ggAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTGACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGATAAGCaa                                                         <  2:341314/149‑1 (MQ=17)                                                                                        tttGTGATTCATGACTGGGGTGAAGTCATAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGTCGTCTTGCGAAGCa                                         >  2:385179/1‑149 (MQ=17)                                                                                          tgtgATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGACAGTGAAAAGCAAGGCGTCTTGCGAAGCAGa                                       <  1:214098/149‑1 (MQ=17)                                                                                                tCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATa                                 >  1:76780/1‑149 (MQ=32)                                                                                                tCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATa                                 >  1:98195/1‑149 (MQ=32)                                                                                                 cATGACTGGGGTGAAGTCGTAACAAGGTAACCGTATGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTTCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAACGCAAGGCGTCTTGCTAAGCAGACTGATAc                                >  2:141545/1‑149 (MQ=255)                                                                                                            aggtaaggaggtgatccaaccgTAGGGGAACCTGCGGt                                                                                                                                    >  1:273245‑M2/23‑38 (MQ=255)                                                                                                            aggtaaggaggtgatccaaccgTAGGGGAACCTGCGGt                                                                                                                                    <  2:273245‑M2/16‑1 (MQ=255)                                                                                                                gTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTcc                            >  1:280935/1‑138 (MQ=33)                                                                                                                gTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTcc                            <  2:280935/138‑1 (MQ=33)                                                                                                                        aaGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAg         >  1:180691/1‑149 (MQ=32)                                                                                                                               ccgTAGGGTACCCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCg  <  1:99785‑M2/146‑1 (MQ=255)                                                                                                                                cgTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAg         <  1:25703‑M2/139‑1 (MQ=255)                                                                                                                                cgTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAg         >  2:25703‑M2/3‑141 (MQ=255)                                                                                                                                  |                                                                                                                                                  CCGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCG  >  NC_000913_3_vch_pgi/4209529‑4209803

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