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.210)	32 (0.970)	26/146	0.3	86.7%	coding (324/1413 nt)	pykF	pyruvate kinase
	?	REL606	= 1733295	3 (0.090)					coding (331/1413 nt)	pykF	pyruvate kinase

TCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1733362‑1733288 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGCTT  <  REL606/1733295‑1733221                                                                                                                                                        TCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTC                                                                            >  2:276454/1‑76   CATACCGATCAGTCTATGGTCAAGNAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGNAACCTTCATACGTTCAT                                                                          <  2:889913/76‑1     nACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATAC                                                                        <  1:925776/75‑1      ACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACG                                                                       <  1:75139/76‑1       CCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGT                                                                      <  2:470901/76‑1                     GTCAACCAGTACGGTGTTGCCTACAGACAGGTCAGTAGTGAAACCTTCATANGTTCATACGTTACCGCAACCATTT                                                        <  2:186996/76‑1                           CAGTACGGTGTTGCCAACAGACAGGTCAGCAGTGAAACCTTCATACGTTCATGCGTTACCGCAACCATTGCGCTGT                                                  >  1:936811/1‑76                            AGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTT                                                 <  1:822208/76‑1                              TACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGC                                               <  2:996815/76‑1                                   TGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATA                                          <  1:516961/76‑1                                   TGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATA                                          <  1:464849/76‑1                                   TGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATA                                          >  2:844782/1‑76                                   TGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATA                                          <  1:103507/76‑1                                     TTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGGTGCCGATAAC                                        <  1:522019/76‑1                                         CAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGAT                                    >  1:561939/1‑76                                               ACAGGTCAGTAGTGAAACCTTCGTACGTTCATACGTTACCGCAACCATTTCNCTGTTGCCGATAACAGATTTATCA                              <  2:179603/76‑1                                                 AGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGT                            <  1:477808/76‑1                                                   GTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGG                          >  2:453923/1‑76                                                     CAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTG                        <  1:1008757/76‑1                                                       GTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAA                      <  2:105362/76‑1                                                        TAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAA                     >  1:375564/1‑76                                                            GAAACNTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAA                 <  2:872437/76‑1                                                            GAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAA                 <  1:682816/76‑1                                                            GAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAA                 >  2:79456/1‑76                                                             AAAGCTTCATACGTTCATACGTTACCGCTACCATTTCGCTGTTGCCGATNACAGATTTATCAGTGGTGAAAGTTAA                <  2:935193/76‑1                                                              AACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAG               <  2:305530/76‑1                                                                 CTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTC            >  2:36902/1‑76                                                                   tctTACGTTCATACGTTACCGCAACCNTTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTG          >  2:935818/4‑76                                                                      TACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACC       <  1:249088/76‑1                                                                        CGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAAGAGATTTATCAGTGGTGAAAGTAAAGGTNTGAACAG     <  1:950757/76‑1                                                                        CGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAG     >  2:290708/1‑76                                                                        CGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAAGTCTGACCAG     <  1:139644/76‑1                                                                          TTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGCT   >  1:1007588/1‑76                                                                           TCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGCTT  >  1:1046474/1‑76                                                                                                                                                        TCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1733362‑1733288 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAGCTT  <  REL606/1733295‑1733221

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