Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913_3_pae_tpiA	1,756,121	Δ557 bp	coding (424‑980/1413 nt)	pykF →	pyruvate kinase I

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_pae_tpiA	1756121	1756677	557	60 [0]	[0] 60	pykF	pyruvate kinase I

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 1756120	0 (0.000)	60 (0.730)	46/284	0.8	100%	coding (423/1413 nt)	pykF	pyruvate kinase I
	?	NC_000913_3_pae_tpiA	1756678 =	0 (0.000)					coding (981/1413 nt)	pykF	pyruvate kinase I

CCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGT  >  NC_000913_3_pae_tpiA/1756530‑1756824                                                                                                                                                     |                                                                                                                                                   tttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaac                                                                                                                                                    >  1:1011203‑M2/149‑149 (MQ=255) tttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaac                                                                                                                                                    >  1:146703‑M2/149‑149 (MQ=255)  ttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaca                                                                                                                                                   <  2:1169742‑M2/2‑1 (MQ=255)   tcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCAt                                                                                                                                                  >  1:155585‑M2/147‑149 (MQ=255)      ccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGc                                                                                                                                               <  2:262568‑M2/6‑1 (MQ=255)      ccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGc                                                                                                                                               <  2:629360‑M2/6‑1 (MQ=255)       cactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtttgatcggtatggaagttaccgccattgaaggtaacacaCATGGCg                                                                                                                                              >  2:418915‑M2/143‑149 (MQ=255)        actgataaatctgttatcggcaacagcgaaatggttgcggtcacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGa                                                                                                                                             <  2:908423‑M2/8‑1 (MQ=255)        actgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGa                                                                                                                                             <  2:155585‑M2/8‑1 (MQ=255)            ataaacctgttatcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCAt                                                                                                                                         <  2:1425313‑M2/12‑1 (MQ=255)             taaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATc                                                                                                                                        <  1:1150319‑M2/13‑1 (MQ=255)              aaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCt                                                                                                                                       >  1:619824‑M2/136‑149 (MQ=255)                atctgttatcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGc                                                                                                                                     <  1:68485‑M2/16‑1 (MQ=255)                 tctgttatcgggaacagttgaatggttgtggtatcgtatgaaggtttcactgctgacctgtgtgttggcaacaccgtacgggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCg                                                                                                                                    <  2:553189‑M2/17‑1 (MQ=255)                    ggtaacggaaacagcgaaatggttacggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAAc                                                                                                                                 <  1:362335‑M2/20‑1 (MQ=255)                        tcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaccaaaCATGGCGACCATCTGCGAACGTac                                                                                                                             >  1:88981‑M2/126‑149 (MQ=255)                                agcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACcgcg                                                                                                                     >  2:46898‑M2/118‑149 (MQ=255)                                    aaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatgtaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGtgat                                                                                                                 <  1:680788‑M2/36‑1 (MQ=255)                                    aaatgggtgcggtaacgtatgaaggcttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGtgat                                                                                                                 <  1:243779‑M2/36‑1 (MQ=255)                                              gggaacttatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGt                                                                                                       <  1:408054‑M2/46‑1 (MQ=255)                                                 aacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGtctc                                                                                                    <  2:1033479‑M2/49‑1 (MQ=255)                                                 aacggatgagggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGtctc                                                                                                    <  1:805543‑M2/49‑1 (MQ=255)                                                             tttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAAt                                                                                        <  1:244380‑M2/61‑1 (MQ=255)                                                             tttcacgactgacctcgctgttggcaacaccgtacgggttgaggatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAAt                                                                                        <  2:914573‑M2/61‑1 (MQ=255)                                                              ttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATg                                                                                       <  2:226178‑M2/62‑1 (MQ=255)                                                                cactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaTATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGAc                                                                                     >  1:1285740‑M2/86‑149 (MQ=255)                                                                  ctactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttcccgccattgaaggtaacaaaCATGGCGACCATCTACGACCGTAACGACCGCGTGATGAACAGCCGGCTCGAGTTCAACAATGACaa                                                                                   >  2:184811‑M2/84‑149 (MQ=255)                                                                          ctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAAc                                                                           >  2:212105‑M2/76‑149 (MQ=255)                                                                          ctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAAc                                                                           >  2:336219‑M2/76‑149 (MQ=255)                                                                          ctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAAc                                                                           >  1:1369268‑M2/76‑149 (MQ=255)                                                                            gtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAActa                                                                         >  2:187175‑M2/74‑148 (MQ=255)                                                                             tcttttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTgc                                                                        <  2:92988‑M2/77‑1 (MQ=255)                                                                                                tggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGGTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATg                                                     >  1:254801‑M2/54‑149 (MQ=255)                                                                                                  gttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGcc                                                   >  1:938721‑M2/52‑149 (MQ=255)                                                                                                    tgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGAGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTAAACAATAACAACAGTAAACTGAGCATTAACGAACCAGTATGCCgt                                                 >  2:565390‑M2/50‑149 (MQ=255)                                                                                                     gacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCgtg                                                <  2:223077‑M2/101‑1 (MQ=255)                                                                                                            gtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGt                                         >  1:53277‑M2/42‑149 (MQ=255)                                                                                                             tctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTTTGCCGTGGTGCCGtt                                        >  1:1162015‑M2/41‑149 (MQ=255)                                                                                                              ctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTAAGCCGTAGTGCCGGTg                                       >  2:1154447‑M2/40‑149 (MQ=255)                                                                                                                 atcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGaaa                                    >  1:275823‑M2/37‑149 (MQ=255)                                                                                                                 atcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGGATGCCGTggtgccgcttcca                                    >  2:1375833‑M2/37‑143 (MQ=255)                                                                                                                  gcggtatggaagttaacgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCCCATTACCGAAGCGGTAGGCCGGGGTGCCGTTGAAAc                                   >  2:896134‑M2/36‑149 (MQ=255)                                                                                                                    ggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAAct                                  >  1:562268‑M2/34‑148 (MQ=255)                                                                                                                    ggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAAct                                  <  2:562268‑M2/115‑1 (MQ=255)                                                                                                                     gtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAActgc                                <  1:46898‑M2/117‑1 (MQ=255)                                                                                                                     gtagggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGCCAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAActgc                                <  2:1011203‑M2/117‑1 (MQ=255)                                                                                                                        tggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGa                             >  1:1223864‑M2/30‑149 (MQ=255)                                                                                                                          gaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTcaa                                                                                            <  1:549273‑M2/57‑1 (MQ=255)                                                                                                                          gaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTcaa                                                                                            >  2:549273‑M2/28‑84 (MQ=255)                                                                                                                          gaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGaaa                           >  2:103238‑M2/28‑149 (MQ=255)                                                                                                                          gaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGaaa                           >  2:651647‑M2/28‑149 (MQ=255)                                                                                                                              ttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACt                       <  1:892408‑M2/126‑1 (MQ=255)                                                                                                                               taccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTg                      <  2:1425960‑M2/127‑1 (MQ=255)                                                                                                                               taccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTg                      <  2:1369268‑M2/127‑1 (MQ=255)                                                                                                                                  cgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGAt                   <  2:1159741‑M2/130‑1 (MQ=255)                                                                                                                                   gtcattgaagttaacaaaCATAGCGACCATCTGCGAACGTGCCGTCCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATg                  <  1:184811‑M2/131‑1 (MQ=255)                                                                                                                                      attgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTAaagaaatgaccgacg               >  2:129312‑M2/16‑142 (MQ=255)                                                                                                                                         gaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGCATGCCGTGGTGCCGTTGAAACTGCTGAACAACTGGATCCTCCGc            >  2:198822‑M2/13‑149 (MQ=255)                                                                                                                                                acaaaCATGGCGACCATCTGCGCACGTACCGACCGCATGATGAACAGCCGTCTCGAGTTCAACTATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAATTACTAAAAAACTGGATGCTCCGCTGATCGt     >  2:1382022‑M2/6‑149 (MQ=255)                                                                                                                                                   aaCATGGCGACCATCTGCGAACGTCCCGACCGCGCGATTAGCAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGt  <  1:758703‑M2/147‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGT  >  NC_000913_3_pae_tpiA/1756530‑1756824

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑