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	42 [0]	[0] 42	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)	42 (1.280)	37/286	0.1	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

CCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCT  >  NC_000913_3_pae_tpiA/1755978‑1756267                                                                                                                                               |                                                                                                                                                    ccacTGATAAATCTGTTCTCTGGAACAGCGGACTTGTGTCGGGAACGTAGGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggc                                                                                                                                               <  1:220910‑M1/149‑7 (MQ=255)  cacTGATAAATCTGTTACCGCCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcg                                                                                                                                              >  2:182621‑M1/1‑142 (MQ=255)   acTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcga                                                                                                                                             >  2:398358‑M1/1‑141 (MQ=255)            tCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCCTTGAAGGTACCAAAcatggcgcccatctccg                                                                                                                                    >  2:349649‑M1/1‑132 (MQ=255)             cTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcga                                                                                                                                   <  1:433735‑M1/149‑19 (MQ=255)                  aTCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGGAGGTAACAAAcatggcgaccatctgcgaacgta                                                                                                                              >  1:142967‑M1/1‑126 (MQ=255)                   tCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtac                                                                                                                             <  1:389913‑M1/149‑25 (MQ=255)                         aCAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCAATTGAAGGTAACAAAcatggcgaccatcttcgaacgtaccaaccg                                                                                                                       >  1:96011‑M1/1‑119 (MQ=255)                           aGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTGTGGAAGTTACCGGCATTGAAGGTAACAAAcatggcgaccatctgcgactgtaccgggcggg                                                                                                                     >  2:156077‑M1/1‑117 (MQ=255)                           aGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcg                                                                                                                     >  1:442153‑M1/1‑117 (MQ=255)                              gAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgccaccgtaccgaccgcgtga                                                                                                                  >  2:399898‑M1/1‑114 (MQ=255)                               aaaTGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgat                                                                                                                 >  2:463707‑M1/1‑113 (MQ=255)                                    gTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcagggcgaccatctgcgaacgtaccgaccgcgtgatgaaca                                                                                                            <  1:190114‑M1/149‑42 (MQ=255)                                    gTTGCGGTAACGTATGAAGGTTTCACTACCGAGCTGGAAATGGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaaca                                                                                                            >  1:137549‑M1/1‑108 (MQ=255)                                       gCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagcc                                                                                                         <  1:180143‑M1/149‑45 (MQ=255)                                              cGTATCGAGGTTTCACTCCTGACCTGTCTGTGGGCAACACCGTACGGGTCGACGATGGTCTGTTCGGTATGGAAGGTACCGCCACTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtcccaa                                                                                                  <  2:387952‑M1/149‑52 (MQ=255)                                                 aTGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgcaatgaacagccgtcgagagtt                                                                                               >  2:478416‑M1/1‑95 (MQ=255)                                                    aaGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaa                                                                                            >  2:195846‑M1/1‑92 (MQ=255)                                                    aaGGTTTCACTACTGACCAGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcagggcgaccatctgcgaacgtcccgaccgcgtgctgaacagccgtcgcgagttcta                                                                                            >  2:451647‑M1/1‑92 (MQ=255)                                                      ggTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaaca                                                                                          <  1:340350‑M1/149‑60 (MQ=255)                                                        tttCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaat                                                                                        <  1:554630‑M1/149‑62 (MQ=255)                                                          tCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatga                                                                                      <  1:195846‑M1/149‑64 (MQ=255)                                                              tactGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaac                                                                                  >  1:66354‑M1/1‑82 (MQ=255)                                                                 tGACCTGTCTGTTGGCAACACCGTACTGGTTAAAGAAGTTCTGATCGGAATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgcctcgagttcaaccatgacaaccgc                                                                               >  2:123729‑M1/1‑79 (MQ=255)                                                                            ttGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaa                                                                                                                                  <  2:288160‑M1/87‑20 (MQ=255)                                                                            ttGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaa                                                                                                                                  >  1:288160‑M1/1‑68 (MQ=255)                                                                            tgggcaACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatgtcgaccatctgcgaacgtaccgaccgcgtgatgaccagccgtctcgagctcaacaatgacaaccgtaaactacccat                                                                    >  2:364790‑M1/3‑68 (MQ=255)                                                                                      cGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaacggcgcattaccgaggcg                                                          >  1:149542‑M1/1‑58 (MQ=255)                                                                                                gACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtg                                                <  1:85265‑M1/149‑102 (MQ=255)                                                                                                  cGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggt                                              <  1:398358‑M1/149‑104 (MQ=255)                                                                                                        tCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgtt                                        <  1:64473‑M1/149‑110 (MQ=255)                                                                                                        tCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgtt                                        >  2:535539‑M1/1‑40 (MQ=255)                                                                                                          tGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttga                                      <  1:535539‑M1/149‑112 (MQ=255)                                                                                                            aTCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccggtgaaa                                    >  1:16217‑M1/1‑36 (MQ=255)                                                                                                                   tGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaactgctga                             >  1:300142‑M1/1‑29 (MQ=255)                                                                                                                   tGGAAGTTACCGCCATTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccatggagcagatgaaactgctgc                             >  2:376634‑M1/1‑29 (MQ=255)                                                                                                                                 aTTGAAGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaactgctgaaaaactggatgctc               <  2:300142‑M1/149‑135 (MQ=255)                                                                                                                                     aaGGTAACAAAcatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacactgacaaccgtaaactgcgcattaccgaagcgggatgccgtggtcacgttgaacctgctgaaaaactggctgctccgct           >  2:333533‑M1/1‑11 (MQ=255)                                                                                                                                          aacaaacatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtctgcagtggtgccgttgaaactgctgaaaaacgggatgctcccctgatcg      >  2:324483‑M1/1‑6 (MQ=255)                                                                                                                                             aaacatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatt                                                                   >  2:72569‑M1/1‑3 (MQ=255)                                                                                                                                             aaacatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcatt                                                                   <  1:72569‑M1/85‑83 (MQ=255)                                                                                                                                              aacatggcgaccatctgcgaacgtaccgaccgcgtgatgaacagccgtctcgagttcaacaatgacaaccgtaaactgcgcattaccgaagcggtatgccgtggtgccgttgaaactgctgaaaaactggatgccccgctgatcgtggt  >  1:9427‑M1/1‑2 (MQ=255)                                                                                                                                               |                                                                                                                                                    CCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACCGCCATTGAAGGTAACAAAGTTATCTGTAAAGTGCTGAACAACGGTGACCTGGGCGAAAACAAAGGTGTGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCT  >  NC_000913_3_pae_tpiA/1755978‑1756267

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