New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	3566202 =	14 (0.970)	4 (0.300)	3/244	1.6	24.1%	coding (380/2448 nt)	glgP	glycogen phosphorylase
	?	NC_000913_3_bme_pgi	3566229 =	12 (0.890)					coding (353/2448 nt)	glgP	glycogen phosphorylase

TTCCAGCCAGTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAA  >  NC_000913_3_bme_pgi/3566060‑3566345                                                                                                                                               |                                                                                                                                                ttCCAGCCAGTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGcc                                                                                                                                          >  1:93785/1‑150 (MQ=255)         agtagtCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTAc                                                                                                                                  <  1:26773/150‑1 (MQ=255)                               gCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttcttc                                                                                                            >  1:149522/1‑150 (MQ=255)                                tggccgtacgctctccaacgccatgttgtcgctaggaatttacgaagatgtacagggcgcactggaagcgatggggttaaatctcgaagagctgattgatgaagaaaatgaCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtca                                                                                                           >  3:78475‑M1/112‑150 (MQ=255)                                tggccgtacgctctccaacgccatgttgtcgctaggaatttacgaagatgtacagggcgcactggaagcgatggggttaaatctcgaagagctgattgatgaagaaaatgaCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtca                                                                                                           >  3:78472‑M1/112‑150 (MQ=255)                                              gTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTc                                                                                             <  2:149522/150‑1 (MQ=255)                                                                       tAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGt                                                                    >  1:6811/1‑150 (MQ=255)                                                                         gCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCACTCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTgc                                                                  >  1:149789/1‑150 (MQ=255)                                                                           gaagatgtacagggcgcactggaagcgatggggttaaatctcgaagagctgattgatgaagaaaatgaCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAAccc                                                                                 >  2:167371‑M1/69‑133 (MQ=255)                                                                           gaagatgtacagggcgcactggaagcgatggggttaaatctcgaagagctgattgatgaagaaaatgaCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAAccc                                                                                 <  1:167371‑M1/65‑1 (MQ=255)                                                                                             ccGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCa                                                                      >  3:86336/1‑126 (MQ=255)                                                                                             ccGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCa                                                                      <  4:86336/126‑1 (MQ=255)                                                                                                               ccAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATg                             >  2:144297/1‑149 (MQ=255)                                                                                                               ccAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATg                             <  1:144297/149‑1 (MQ=255)                                                                                                                      aTCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGga                     <  2:6811/150‑1 (MQ=255)                                                                                                                         aaGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTc                                                                                             <  1:3310/75‑1 (MQ=255)                                                                                                                         aaGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTc                                                                                             >  2:3310/1‑75 (MQ=255)                                                                                                                                         aGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCaaa  <  2:149789/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                TTCCAGCCAGTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAA  >  NC_000913_3_bme_pgi/3566060‑3566345

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