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

TTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCT  >  NC_000913_3_pae_tpiA/1755973‑1756267                                                                                                                                                    |                                                                                                                                                    tttCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAc                                                                                                                                                    >  1:1011203‑M1/1‑148 (MQ=255) tttCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAc                                                                                                                                                    >  1:146703‑M1/1‑148 (MQ=255)  ttCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAca                                                                                                                                                   <  2:1169742‑M1/149‑3 (MQ=255)   tCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcat                                                                                                                                                  >  1:155585‑M1/1‑146 (MQ=255)      ccacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggc                                                                                                                                               <  2:262568‑M1/149‑7 (MQ=255)      ccacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggc                                                                                                                                               <  2:629360‑M1/149‑7 (MQ=255)       cacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTTTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACACAcatggcg                                                                                                                                              >  2:418915‑M1/1‑142 (MQ=255)        acTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTCACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcga                                                                                                                                             <  2:908423‑M1/149‑9 (MQ=255)        acTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcga                                                                                                                                             <  2:155585‑M1/149‑9 (MQ=255)            aTAAACCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccat                                                                                                                                         <  2:1425313‑M1/149‑13 (MQ=255)             tAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatc                                                                                                                                        <  1:1150319‑M1/149‑14 (MQ=255)              aaaTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatct                                                                                                                                       >  1:619824‑M1/1‑135 (MQ=255)                aTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgc                                                                                                                                     <  1:68485‑M1/149‑17 (MQ=255)                 tCTGTTATCGGGAACAGTTGAATGGTTGTGGTATCGTATGAAGGTTTCACTGCTGACCTGTGTGTTGGCAACACCGTACGGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcg                                                                                                                                    <  2:553189‑M1/149‑18 (MQ=255)                    ggtaacggaaaCAGCGAAATGGTTACGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaac                                                                                                                                 <  1:362335‑M1/144‑21 (MQ=255)                        tCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTACCAAAcatggcgaccatctgcgaacgtac                                                                                                                             >  1:88981‑M1/1‑125 (MQ=255)                                aGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcg                                                                                                                     >  2:46898‑M1/1‑117 (MQ=255)                                    aaaTGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGTAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgat                                                                                                                 <  1:680788‑M1/149‑37 (MQ=255)                                    aaaTGGGTGCGGTAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgat                                                                                                                 <  1:243779‑M1/149‑37 (MQ=255)                                              gggaacttATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgt                                                                                                       <  1:408054‑M1/146‑47 (MQ=255)                                                 aaCGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctc                                                                                                    <  2:1033479‑M1/149‑50 (MQ=255)                                                 aaCGGATGAGGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctc                                                                                                    <  1:805543‑M1/149‑50 (MQ=255)                                                             tttCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaat                                                                                        <  1:244380‑M1/149‑62 (MQ=255)                                                             tttCACGACTGACCTCGCTGTTGGCAACACCGTACGGGTTGAGGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaat                                                                                        <  2:914573‑M1/149‑62 (MQ=255)                                                              ttCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatg                                                                                       <  2:226178‑M1/149‑63 (MQ=255)                                                                cACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAtatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgac                                                                                     >  1:1285740‑M1/1‑85 (MQ=255)                                                                  ctactGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTCCCGCCATTGAAGGTAACAAAcatggcgaccatctacgaccgtaacgaccgcgtgatgaacagccggctcgagttcaacaatgacaa                                                                                   >  2:184811‑M1/1‑83 (MQ=255)                                                                          ctgtctgtTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaac                                                                           >  2:212105‑M1/1‑75 (MQ=255)                                                                          ctgtctgtTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaac                                                                           >  2:336219‑M1/1‑75 (MQ=255)                                                                          ctgtctgtTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaac                                                                           >  1:1369268‑M1/1‑75 (MQ=255)                                                                            gtctgtTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaacta                                                                         >  2:187175‑M1/1‑73 (MQ=255)                                                                             tctttTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgc                                                                        <  2:92988‑M1/149‑78 (MQ=255)                                                                                                tGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgaggtcaacaatgacaaccgtaaactgcgcattaccgaagcggtatg                                                     >  1:254801‑M1/1‑53 (MQ=255)                                                                                                  gTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgcc                                                   >  1:938721‑M1/1‑51 (MQ=255)                                                                                                    tGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgagaacgtaccgaccgcgtgatgaacagccgtctcgagttaaacaataacaacagtaaactgagcattaacgaaccagtatgccgt                                                 >  2:565390‑M1/1‑49 (MQ=255)                                                                                                     gACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtg                                                <  2:223077‑M1/149‑102 (MQ=255)                                                                                                            gTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgt                                         >  1:53277‑M1/1‑41 (MQ=255)                                                                                                             tCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtttgccgtggtgccgtt                                        >  1:1162015‑M1/1‑40 (MQ=255)                                                                                                              cTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtaagccgtagtgccggtg                                       >  2:1154447‑M1/1‑39 (MQ=255)                                                                                                                 aTCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaa                                    >  1:275823‑M1/1‑36 (MQ=255)                                                                                                                 aTCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcgggatgccgtggtgccgcttcca                                    >  2:1375833‑M1/1‑36 (MQ=255)                                                                                                                  gcgGTATGGAAGTTAACGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcccattaccgaagcggtaggccggggtgccgttgaaac                                   >  2:896134‑M1/2‑35 (MQ=255)                                                                                                                    ggTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaact                                  >  1:562268‑M1/1‑33 (MQ=255)                                                                                                                    ggTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaact                                  <  2:562268‑M1/148‑116 (MQ=255)                                                                                                                     gTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaactgc                                <  1:46898‑M1/149‑118 (MQ=255)                                                                                                                     gTAGGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgccaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaactgc                                <  2:1011203‑M1/149‑118 (MQ=255)                                                                                                                        tGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaactgctga                             >  1:1223864‑M1/1‑29 (MQ=255)                                                                                                                          gAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaactgctgaaa                           >  2:103238‑M1/1‑27 (MQ=255)                                                                                                                          gAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaactgctgaaa                           >  2:651647‑M1/1‑27 (MQ=255)                                                                                                                          gAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaa                                                                                            <  1:549273‑M1/84‑58 (MQ=255)                                                                                                                          gAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaa                                                                                            >  2:549273‑M1/1‑27 (MQ=255)                                                                                                                              ttACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaactgctgaaaaact                       <  1:892408‑M1/149‑127 (MQ=255)                                                                                                                               tACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaactgctgaaaaactg                      <  2:1425960‑M1/149‑128 (MQ=255)                                                                                                                               tACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaactgctgaaaaactg                      <  2:1369268‑M1/149‑128 (MQ=255)                                                                                                                                  cGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaactgctgaaaaactggat                   <  2:1159741‑M1/149‑131 (MQ=255)                                                                                                                                   gtcatTGAAGTTAACAAAcatagcgaccatctgcgaacgtgccgtccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaactgctgaaaaactggatg                  <  1:184811‑M1/147‑132 (MQ=255)                                                                                                                                      aTTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaactgctaaagaaatgaccgacg               >  2:129312‑M1/1‑15 (MQ=255)                                                                                                                                         gAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggcatgccgtggtgccgttgaaactgctgaacaactggatcctccgc            >  2:198822‑M1/1‑12 (MQ=255)                                                                                                                                                acaaacatggcgaccatctgcgcacgtaccgaccgcatgatgaacagccgtctcgagttcaactatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaattactaaaaaactggatgctccgctgatcgt     >  2:1382022‑M1/1‑5 (MQ=255)                                                                                                                                                   aacatggcgaccatctgcgaacgtcccgaccgcgcgattagcagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaactgctgaaaaactggatgctccgctgatcgtggt  <  1:758703‑M1/149‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCT  >  NC_000913_3_pae_tpiA/1755973‑1756267

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: ‑