Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	REL606	1,733,297	(CGTATGAA)1→2	coding (333/1413 nt)	pykF →	pyruvate kinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	1733288 =	3 (0.140)	19 (0.900)	18/146	0.3	88.5%	coding (324/1413 nt)	pykF	pyruvate kinase
	?	REL606	= 1733295	2 (0.100)					coding (331/1413 nt)	pykF	pyruvate kinase

GATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1733354‑1733288 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTG  <  REL606/1733295‑1733229                                                                                                                                        GATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTA                                                            <  2:166283/76‑1      GACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAGCCTTCATACGTTCATACGTTACCGCA                                                       <  2:571187/76‑1           TCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCAT                                                  <  1:182326/76‑1               CAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCG                                              <  1:170140/76‑1               CAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCG                                              >  1:284956/1‑76                AACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGC                                             >  2:526990/1‑76                  CCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTG                                           >  1:438223/1‑76                     GTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTG                                        >  1:286512/1‑76                     GTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTG                                        <  1:85329/76‑1                       ACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCC                                      >  2:149409/1‑76                           TGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATA                                  <  2:606336/76‑1                                  AACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATT                           <  1:172718/76‑1                                                TAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAA             >  2:6130/1‑76                                                 AGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAG            <  2:628610/76‑1                                                     AAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAA        <  1:157213/76‑1                                                     AAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAA        <  1:157211/76‑1                                                         CTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTC    <  1:32947/76‑1                                                           TCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTG  <  2:305771/76‑1                                                           TCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTG  >  1:95211/1‑76                                                                                                                                        GATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1733354‑1733288 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTG  <  REL606/1733295‑1733229

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