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

GCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGAC  >  NC_000913_3_bme_pgi/3566090‑3566357                                                                                                                                            |                                                                                                                                 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)                                                                           gggttaaatctcgaagagctgattgatgaagaaaatgacccaggcctcggtaacggtggcctgggTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTc                                                               >  2:167371‑M2/66‑133 (MQ=255)                                                                           gggttaaatctcgaagagctgattgatgaagaaaatgacccaggcctcggtaacggtggcctgggTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTc                                                               <  1:167371‑M2/68‑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)                                                                                                     tgaagaaaatgacccaggcctcggtaacggtggcctgggTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCa                    >  3:78472‑M2/40‑150 (MQ=255)                                                                                                     tgaagaaaatgacccaggcctcggtaacggtggcctgggTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCa                    >  3:78475‑M2/40‑150 (MQ=255)                                                                                                           aGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCaaa              <  2:149789/150‑1 (MQ=255)                                                                                                                       CCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTTCATCGAc  >  2:168036‑M1/1‑150 (MQ=255)                                                                                                                                            |                                                                                                                                 GCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGAC  >  NC_000913_3_bme_pgi/3566090‑3566357

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.