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

ATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGT  >  NC_000913_3_pae_tpiA/1756535‑1756824                                                                                                                                                |                                                                                                                                                   ccactgataaatctgttctctggaacagcggacttgtgtcgggaacgtaggaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGc                                                                                                                                               <  1:220910‑M2/6‑1 (MQ=255)  cactgataaatctgttaccgccaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCg                                                                                                                                              >  2:182621‑M2/143‑149 (MQ=255)   actgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGa                                                                                                                                             >  2:398358‑M2/142‑149 (MQ=255)            tctgttatcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgcccttgaaggtaccaaaCATGGCGCCCATCTCCg                                                                                                                                    >  2:349649‑M2/133‑149 (MQ=255)             ctgttatcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGa                                                                                                                                   <  1:433735‑M2/18‑1 (MQ=255)                  atcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattggaggtaacaaaCATGGCGACCATCTGCGAACGTa                                                                                                                              >  1:142967‑M2/127‑149 (MQ=255)                   tcggcaacagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTac                                                                                                                             <  1:389913‑M2/24‑1 (MQ=255)                         acagcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgcaattgaaggtaacaaaCATGGCGACCATCTTCGAACGTaccaaccg                                                                                                                       >  1:96011‑M2/120‑149 (MQ=255)                           agcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtgtggaagttaccggcattgaaggtaacaaaCATGGCGACCATCTGCgactgtaccgggcggg                                                                                                                     >  2:156077‑M2/118‑143 (MQ=255)                           agcgaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACcgcg                                                                                                                     >  1:442153‑M2/118‑149 (MQ=255)                              gaaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCCACCGTACCGACCGCGtga                                                                                                                  >  2:399898‑M2/115‑149 (MQ=255)                               aaatggttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGtgat                                                                                                                 >  2:463707‑M2/114‑149 (MQ=255)                                    gttgcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCAGGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACa                                                                                                            <  1:190114‑M2/41‑1 (MQ=255)                                    gttgcggtaacgtatgaaggtttcactaccgagctggaaatgggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACa                                                                                                            >  1:137549‑M2/109‑149 (MQ=255)                                       gcggtaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGcc                                                                                                         <  1:180143‑M2/44‑1 (MQ=255)                                              cgtatcgaggtttcactcctgacctgtctgtgggcaacaccgtacgggtcgacgatggtctgttcggtatggaaggtaccgccactgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAgccgtcccaa                                                                                                  <  2:387952‑M2/51‑5 (MQ=255)                                                 atgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGCAATGAACAGCCGTCGAGAGtt                                                                                               >  2:478416‑M2/96‑149 (MQ=255)                                                    aaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTcaa                                                                                            >  2:195846‑M2/93‑149 (MQ=255)                                                    aaggtttcactactgaccagtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCAGGGCGACCATCTGCGAACGTCCCGACCGCGTGCTGAACAGCCGTCGCGAGttcta                                                                                            >  2:451647‑M2/93‑147 (MQ=255)                                                      ggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTcaaca                                                                                          <  1:340350‑M2/59‑1 (MQ=255)                                                        tttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAAt                                                                                        <  1:554630‑M2/61‑1 (MQ=255)                                                          tcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGa                                                                                      <  1:195846‑M2/63‑1 (MQ=255)                                                              tactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAAc                                                                                  >  1:66354‑M2/83‑149 (MQ=255)                                                                 tgacctgtctgttggcaacaccgtactggttaaagaagttctgatcggaatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGCCTCGAGTTCAACCATGACAACcgc                                                                               >  2:123729‑M2/80‑148 (MQ=255)                                                                            ttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGaa                                                                                                                                  <  2:288160‑M2/19‑1 (MQ=255)                                                                            ttggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGaa                                                                                                                                  >  1:288160‑M2/69‑87 (MQ=255)                                                                            tgggcaacaccgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGTCGACCATCTGCGAACGTACCGACCGCGTGATGACCAGCCGTCTCGAGCTCAACAATGACAACCGTAAACTACCCAt                                                                    >  2:364790‑M2/69‑149 (MQ=255)                                                                                      cgtactggttgacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACGGCGCATTACCGAGGCg                                                          >  1:149542‑M2/59‑149 (MQ=255)                                                                                                gacgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCgtg                                                <  1:85265‑M2/101‑1 (MQ=255)                                                                                                  cgatggtctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCgtggt                                              <  1:398358‑M2/103‑1 (MQ=255)                                                                                                        tctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGtt                                        <  1:64473‑M2/109‑1 (MQ=255)                                                                                                        tctgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGtt                                        >  2:535539‑M2/41‑149 (MQ=255)                                                                                                          tgatcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGa                                      <  1:535539‑M2/111‑1 (MQ=255)                                                                                                            atcggtatggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGGTGaaa                                    >  1:16217‑M2/37‑149 (MQ=255)                                                                                                                   tggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGa                             >  1:300142‑M2/30‑149 (MQ=255)                                                                                                                   tggaagttaccgccattgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCATGGAGCAGATGAAactgctgc                             >  2:376634‑M2/30‑148 (MQ=255)                                                                                                                                 attgaaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTc               <  2:300142‑M2/134‑1 (MQ=255)                                                                                                                                     aaggtaacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACACTGACAACCGTAAACTGCGCATTACCGAAGCGGGATGCCGTGGTCACGTTGAACCTGCTGAAAAACTGGCTGCTCCGCt           >  2:333533‑M2/12‑149 (MQ=255)                                                                                                                                          aacaaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTCTGCAGTGGTGCCGTTGAAACTGCTGAAAAACGGGATGCTCCCCTGATCg      >  2:324483‑M2/7‑149 (MQ=255)                                                                                                                                             aaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCAtt                                                                   >  2:72569‑M2/4‑85 (MQ=255)                                                                                                                                             aaaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCAtt                                                                   <  1:72569‑M2/82‑1 (MQ=255)                                                                                                                                              aaCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCCCCGCTGATCGTGGt  >  1:9427‑M2/3‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                   ATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGT  >  NC_000913_3_pae_tpiA/1756535‑1756824

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