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 =	4 (0.130)	29 (1.010)	26/146	0.2	92.2%	coding (324/1413 nt)	pykF	pyruvate kinase
	?	REL606	= 1733295	1 (0.030)					coding (331/1413 nt)	pykF	pyruvate kinase

CTTCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1733364‑1733288 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGCTT  <  REL606/1733295‑1733221                                                                                                                                                          CTTCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCNAACAGAAAGGTCAGTAGTGAAACCTTCATACGT                                                                              >  1:778529/1‑76  TTCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTTAGTAGTGAAACCTTCATACGTT                                                                             >  1:10268/1‑76  TTCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT                                                                             >  2:194355/1‑76         CCGATCAGACCATCGTCAACCAGTACGATGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGT                                                                      >  1:89474/1‑76                   CATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTTAAACCCTCATACGTTCATACGTTACCGCAACC                                                            <  1:850270/76‑1                   CATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAAgc                                                            >  2:102717/1‑74                       GTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTT                                                        <  2:966488/76‑1                       GTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTT                                                        <  2:204652/76‑1                         CAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCG                                                      <  2:89474/76‑1                              AGTACGGTGTTGCCAACAGACATGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTT                                                 >  1:509557/1‑76                              AGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTT                                                 <  1:797508/76‑1                              AGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTT                                                 <  2:410433/76‑1                               GTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTG                                                <  1:43275/76‑1                                TACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGC                                               <  1:587952/76‑1                                 ACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCn                                              >  1:842273/1‑75                                 ACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCC                                              >  2:122622/1‑76                                   GGTGTTGCCAACAGACAGGTCAGTAGTCAAANTTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGA                                            >  1:78884/1‑76                                       TTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAAC                                        <  2:131528/76‑1                                        TGCCAACAGACATGGCATTAGTGANACCTTCATACGTTCATACGTTACCGCAACCATTNCGCTGTTGCCGATAACA                                       <  2:793795/76‑1                                          CCAACAGACAGGTCTGTAGTGAAACCTTCATACGTTCATACGTTACCTCAACCATTTCGCTGTTGCCGATAACAGA                                     <  2:937548/76‑1                                               AGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTAT                                >  2:467239/1‑76                                                     GTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGG                          <  1:227837/76‑1                                                        AGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGA                       >  2:18841/1‑76                                                         GTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAA                      <  2:704827/76‑1                                                         GTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAA                      >  1:719861/1‑76                                                             TGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTA                  <  2:719861/76‑1                                                               AAACCTTNATACGTTCATACGTCANCGCAACCATTTCGCTGNTGCCGATATCAGATTTATCAGTGGTGAAAGTAAA                <  2:856292/76‑1                                                                    TTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCT           <  1:649075/76‑1                                                                     TCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTG          <  1:483844/76‑1                                                                      nATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGA         <  1:845701/75‑1                                                                       ATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGAC        >  2:363692/1‑76                                                                          CGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAG     <  1:443191/76‑1                                                                            TTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGCT   >  1:398517/1‑76                                                                             TCATACGTTACCGCAACCATTTCGCTGTTCCCGATAACAGAATTATCNGTGGTGNAAGTAAAGGTCTGACCAGCTn  >  1:805847/1‑75                                                                                                                                                          CTTCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1733364‑1733288 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGCTT  <  REL606/1733295‑1733221

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