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 =	2 (0.050)	33 (0.900)	28/146	0.4	91.7%	coding (324/1413 nt)	pykF	pyruvate kinase
	?	REL606	= 1733295	4 (0.110)					coding (331/1413 nt)	pykF	pyruvate kinase

CATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1733360‑1733288 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAG  <  REL606/1733295‑1733224                                                                                                                                                   CATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCAT                                                                       <  1:947641/76‑1   TACCGATCAGACCATCGTCAACAAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATAC                                                                     >  2:530823/1‑76    ACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACG                                                                    <  2:702037/76‑1      CGATCAGACCATCGTCAACCAGNNCGGTGTTGCCNNGAGACAGGACAGTAGTGAAACCTTCATACGTTCATACGTT                                                                  <  1:1087750/76‑1       GATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTA                                                                 >  2:832183/1‑76           AGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGC                                                             <  1:942686/76‑1               CATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACC                                                         <  2:342019/76‑1                  CGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATT                                                      >  1:663096/1‑76                       ACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCT                                                 <  2:1176271/76‑1                          AGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTT                                              <  2:322325/76‑1                                 TGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATA                                       >  1:356033/1‑76                                  GTTGCCAACAGACANGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAn                                      <  2:169820/76‑2                                   TTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGNTGTTGCCGATgaa                                     >  1:163634/1‑73                                    TGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACA                                    >  2:78046/1‑76                                      CCAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGA                                  <  1:1011293/76‑1                                       CAACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGAT                                 <  1:238603/76‑1                                        AACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATT                                >  1:467945/1‑76                                         ACAGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTT                               <  1:1070606/76‑1                                           AGACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTAT                             >  2:563124/1‑76                                            GACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATn                            <  1:1024103/76‑2                                             ACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCA                           >  1:518396/1‑76                                             ACAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGATGTTGCCGATAACAGATTTATCA                           <  1:721358/76‑1                                              CAGGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAG                          <  1:810431/76‑1                                                GGTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATGTCGCTGTTGCCGATAACAGATTTATCAGTG                        <  1:1095871/76‑1                                                 GTCAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGG                       >  1:91241/1‑76                                                   CAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTG                     <  1:696910/76‑1                                                   CAGTAGTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTG                     <  1:656765/76‑1                                                        GTGAAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGT                <  1:960052/76‑1                                                           AAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAA             >  1:856687/1‑76                                                           AAACCTTCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAA             >  2:762910/1‑76                                                                 TCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTG       >  2:562345/1‑76                                                                 TCATACGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTG       <  1:326830/76‑1                                                                      CGTTCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAG  <  1:813556/76‑1                                                                                                                                                   CATACCGATCAGACCATCGTCAACCAGTACGGTGTTGCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1733360‑1733288 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttCATACGTTACCGCAACCATTTCGCTGTTGCCGATAACAGATTTATCAGTGGTGAAAGTAAAGGTCTGACCAG  <  REL606/1733295‑1733224

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