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 =	7 (0.150)	61 (1.310)	43/146	0.1	90.5%	coding (324/1413 nt)	pykF	pyruvate kinase
	?	REL606	= 1733295	6 (0.130)					coding (331/1413 nt)	pykF	pyruvate kinase

TTCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1733363‑1733288 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGCTT  <  REL606/1733295‑1733221                                                                                                                                                         TTCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT                                                                             <  1:560056/76‑1   CCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCA                                                                           >  2:418988/1‑76   CCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCA                                                                           <  2:239535/76‑1    CATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCAT                                                                          >  1:1079135/1‑76    CATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCAT                                                                          <  1:1321025/76‑1       ACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACG                                                                       <  1:445673/76‑1         CGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTT                                                                     <  2:724020/76‑1         CGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTT                                                                     >  1:670459/1‑76         CGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTT                                                                     <  2:1048164/76‑1            TCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACC                                                                  >  2:308791/1‑76             CAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCG                                                                 <  1:1113676/76‑1             CAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCG                                                                 >  1:806049/1‑76             CAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCG                                                                 <  1:461068/76‑1               GACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCA                                                               >  1:816870/1‑76                ACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAA                                                              <  1:831173/76‑1                 CCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAAC                                                             <  1:699326/76‑1                 CCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAAC                                                             <  2:1363348/76‑1                  CATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACC                                                            <  2:664656/76‑1                    TCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCAT                                                          >  1:195446/1‑76                    TCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCAT                                                          >  2:471473/1‑76                      GTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTT                                                        >  2:1440324/1‑76                       TCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTC                                                       >  2:1203130/1‑76                              GTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTG                                                >  1:1069642/1‑76                              GTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTG                                                >  1:228128/1‑76                                ACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCC                                              >  2:1415890/1‑76                                 CGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCG                                             <  2:1019145/76‑1                                    TGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATA                                          >  2:457806/1‑76                                     GTTGCCAGCAGACAGGTCAGTAGTGAAACCTGCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAA                                         <  2:1140605/76‑1                                     GTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAA                                         <  2:432412/76‑1                                     GTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAA                                         <  2:319148/76‑1                                     GTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAA                                         >  2:1457373/1‑76                                     GTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAA                                         >  2:1368878/1‑76                                      TTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAAC                                        <  1:622848/76‑1                                      TTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAAC                                        <  2:1210146/76‑1                                       TGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACA                                       <  1:697954/76‑1                                       TGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACA                                       <  1:646920/76‑1                                       TGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACA                                       >  2:1294676/1‑76                                        GCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAG                                      <  1:1199518/76‑1                                          CAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGAT                                    <  2:931578/76‑1                                            ACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTT                                  <  1:67846/76‑1                                                ACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCA                              <  1:34623/76‑1                                                ACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCA                              >  1:1340346/1‑76                                                   GGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTG                           <  1:125700/76‑1                                                    GTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGG                          <  2:284862/76‑1                                                          AGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAG                    <  2:505573/76‑1                                                          AGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAG                    <  2:1503177/76‑1                                                            TGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTA                  <  2:1087626/76‑1                                                             GAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAA                 <  2:1503329/76‑1                                                             GAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAA                 <  2:238516/76‑1                                                             GAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAA                 <  1:237525/76‑1                                                              AAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAA                <  1:395056/76‑1                                                              AAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAA                <  2:420225/76‑1                                                                    TCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTG          >  1:1470545/1‑76                                                                     CATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGGTAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGA         <  1:652717/76‑1                                                                     CATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGA         >  2:1065758/1‑76                                                                     CATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGA         >  2:1065769/1‑76                                                                     CATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGA         >  2:1080938/1‑76                                                                     CATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGA         >  1:1208479/1‑76                                                                        ACGTTCATACGTTACCGTAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCA      <  2:1365364/76‑1                                                                        ACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCA      >  1:941413/1‑76                                                                        ACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCA      >  1:931339/1‑76                                                                         CGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAG     <  1:308791/76‑1                                                                            TCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGCTT  >  1:664584/1‑76                                                                                                                                                         TTCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1733363‑1733288 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGCTT  <  REL606/1733295‑1733221

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