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 =	8 (0.220)	48 (1.340)	39/146	0.0	91.6%	coding (324/1413 nt)	pykF	pyruvate kinase
	?	REL606	= 1733295	1 (0.030)					coding (331/1413 nt)	pykF	pyruvate kinase

CTTCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1733364‑1733288 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAG  <  REL606/1733295‑1733224                                                                                                                                                       CTTCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGT                                                                           <  1:894114/76‑1   TCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTC                                                                         <  1:308345/76‑1     CATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCAT                                                                       <  1:108500/76‑1     CATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCAT                                                                       <  1:1127482/76‑1      ATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATA                                                                      >  2:630088/1‑76            ATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTAC                                                                <  1:927435/76‑1              CAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCG                                                              >  2:885290/1‑76              CAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCG                                                              >  2:848789/1‑76                aACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCA                                                            <  2:108252/75‑1                 ACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCCTACGTTACCGCAA                                                           <  1:77386/76‑1                  CCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAAC                                                          <  1:496098/76‑1                    ATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATGCGTTCATACGTTACCGCAACCA                                                        >  2:32876/1‑76                    ATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCA                                                        >  1:177960/1‑76                     TCGTCAATCAGTACGGTGTTACNNACAGACAGGTNNATAGTGAACCCTTCGTACGTTCATACGTTACCGCAACCAT                                                       <  1:1049220/76‑1                     TCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCAT                                                       <  2:576013/76‑1                      CGTCAACCAATACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCACAACCATT                                                      <  1:32876/76‑1                              AGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTT                                              <  1:349080/76‑1                                   GGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGA                                         >  1:117283/1‑76                                    GTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATGCGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGAT                                        <  1:819000/76‑1                                    GTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGAT                                        <  2:706769/76‑1                                     TGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATA                                       <  1:758031/76‑1                                      GTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAA                                      <  1:46413/76‑1                                      GTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAA                                      <  2:582838/76‑1                                       TTGCCAACAGACAGGTCGGTAGTGAAACCTTCATACGTTCATACGTTACCGCAATCATTTCGCTGTTGCCGATAAC                                     <  1:131643/76‑1                                         GCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATGCGTTACCGCAACCATTTCGCTGTTGCCGATAACAG                                   <  1:318808/76‑1                                          CCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGA                                  <  2:219597/76‑1                                               AGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACACATTTAT                             <  1:71313/76‑1                                                GACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATC                            >  1:966380/1‑76                                                 ACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCA                           >  2:243290/1‑76                                                   AGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATATCGCTGTTGCCGATAACAGATTTATCAGT                         <  2:52222/76‑1                                                    GGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTG                        >  2:699172/1‑76                                                         GTAGTGAAACCTTCNTACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAn                   <  2:167140/76‑2                                                         GTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGCGAA                   <  2:513125/76‑1                                                          TAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAA                  >  1:485933/1‑76                                                          TAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAA                  <  1:253630/76‑1                                                           AGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAG                 >  1:830722/1‑76                                                           AGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAG                 >  1:635866/1‑76                                                           AGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAG                 <  1:572672/76‑1                                                            GTGAAACCTTCATACGTTCATANNTTACCGCAACNNGTTCGCTGCTGCCGTTAACAGATTTATCAGTGGTGAAAGT                <  1:1053067/76‑1                                                            GTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGT                >  1:940594/1‑76                                                             TGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTA               <  2:1120569/76‑1                                                              GAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAA              >  1:630522/1‑76                                                              GAAACCTTCATACGGTCATACGTTACCGTAANCATTTCGCGGTTGCNGAT‑ACNGATTTATCAGTGGTGAAAGTAAA             <  2:148197/76‑1                                                               AAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAA             >  2:287658/1‑76                                                                   CTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTC         >  1:842147/1‑76                                                                      CATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGA      >  2:205003/1‑76                                                                      CATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGA      >  1:549299/1‑76                                                                         ACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCA   <  1:752981/76‑1                                                                          CGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAG  <  2:848992/76‑1                                                                                                                                                       CTTCCATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1733364‑1733288 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAG  <  REL606/1733295‑1733224

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