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 =	13 (0.230)	57 (1.030)	38/146	0.4	87.3%	coding (324/1413 nt)	pykF	pyruvate kinase
	?	REL606	= 1733295	4 (0.070)					coding (331/1413 nt)	pykF	pyruvate kinase

TCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1733362‑1733288 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGCTT  <  REL606/1733295‑1733221                                                                                                                                                        aCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTC                                                                            >  1:1024269/2‑76     TACCGATCAGACCATCGTCAACCAGTACGGGGTTGCCAACAGACAGGGCAGTAGTGAAACCTTCATACGTTCATAC                                                                        >  2:988796/1‑76      ACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAAACAGGTCAGTAGTGAAACCTTCATACGTTCATACG                                                                       >  1:836583/1‑76        CGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTT                                                                     >  1:1088459/1‑76         GATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTA                                                                    <  2:1139772/76‑1         GATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCAGACGTTCATACGTTA                                                                    <  1:983582/76‑1          ATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTAC                                                                   >  2:614759/1‑76          ATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTAC                                                                   >  2:344172/1‑76            CAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACNTTCATACGTTACCG                                                                 <  1:1000316/76‑1             AGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGC                                                                >  1:206265/1‑76                   TCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACGAT                                                          <  1:678977/76‑1                    CGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATT                                                         >  1:503165/1‑76                    CGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATT                                                         >  1:456868/1‑76                       CAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCG                                                      >  1:1002901/1‑76                          CCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTG                                                   <  2:894071/76‑1                           CAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGT                                                  >  1:553030/1‑76                            AGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTT                                                 <  2:1444139/76‑1                             GTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTG                                                <  2:1505932/76‑1                                   TGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATATGTTACCGCAACCATTTCGCTGTTGCCGATA                                          >  1:1731856/1‑76                                   TGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATA                                          >  2:514430/1‑76                                   TGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATA                                          <  1:1337532/76‑1                                    GTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAA                                         <  2:810577/76‑1                                    GTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAA                                         <  2:198399/76‑1                                     TTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAAC                                        <  2:1396340/76‑1                                     TTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAAC                                        >  1:827961/1‑76                                     TTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAAC                                        >  1:1452506/1‑76                                      cGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCAGATCACA                                       >  2:200091/2‑76                                      TGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACA                                       <  1:1458173/76‑1                                       GCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAG                                      <  2:1185041/76‑1                                        CCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGA                                     <  2:994105/76‑1                                         CAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGAT                                    >  2:1482569/1‑76                                           ACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTT                                  <  2:804197/76‑1                                           ACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTT                                  <  1:1043383/76‑1                                           ACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTT                                  <  2:460683/76‑1                                           ACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTT                                  <  2:460664/76‑1                                            CAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTA                                 >  2:256953/1‑76                                             AGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTAT                                <  2:864016/76‑1                                                   GTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGG                          <  2:826454/76‑1                                                     CAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTG                        <  2:1049567/76‑1                                                     CAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTG                        <  2:1076178/76‑1                                                          GTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAATGGTGAAAGT                   <  1:1697203/76‑1                                                           TGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTA                  <  2:1673082/76‑1                                                           TGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTA                  >  1:218228/1‑76                                                            GAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAA                 <  1:2998/76‑1                                                            GAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAA                 <  1:1310014/76‑1                                                                CCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGT             >  2:507448/1‑76                                                                  TTCATACGTTCATACGTTACCGCAACCATTTAGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCT           <  1:94280/76‑1                                                                   TCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTG          <  1:1307939/76‑1                                                                    CATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGA         <  2:1024269/76‑1                                                                    CATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGA         >  1:150476/1‑76                                                                    CATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGA         <  1:723981/76‑1                                                                         GTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAGGTAAAGGTCTGACCAGC    >  2:1334097/1‑76                                                                         GTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGC    <  2:1387649/76‑1                                                                         GTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGC    <  1:123123/76‑1                                                                         GTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGC    >  1:351332/1‑76                                                                         GTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGC    <  1:185135/76‑1                                                                         GTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGC    >  1:1821782/1‑76                                                                           TCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGCTT  <  1:1524263/76‑1                                                                                                                                                        TCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1733362‑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 < 34 ≤ 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.