Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913_3_pae_tpiA	1,483,617	Δ5 bp	coding (557‑561/3903 nt)	hrpA →	putative ATP‑dependent helicase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_pae_tpiA	1483617	1483621	5	24 [0]	[0] 24	hrpA	putative ATP‑dependent helicase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 1483616	0 (0.000)	24 (0.660)	22/282	0.8	100%	coding (556/3903 nt)	hrpA	putative ATP‑dependent helicase
	?	NC_000913_3_pae_tpiA	1483622 =	0 (0.000)					coding (562/3903 nt)	hrpA	putative ATP‑dependent helicase

TGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACT  >  NC_000913_3_pae_tpiA/1483473‑1483753                                                                                                                                                |                                                                                                                                          tGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgat                                                                                                                                      >  1:597666‑M1/1‑144 (MQ=255)     gAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgatgcag                                                                                                                                  <  2:549881‑M1/149‑10 (MQ=255)        ccGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgatgcagtac                                                                                                                               >  1:405493‑M1/1‑137 (MQ=255)                       gCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGCTCACGAAAGTGCACAAACGATGTACAACCTGATAACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgatgcagtacgacactatcattatt                                                                                                                <  2:547329‑M1/149‑28 (MQ=255)                         tGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgatgcagtacgacactatcattattga                                                                                                              >  1:391995‑M1/1‑120 (MQ=255)                         tGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgatgcagtacgacactatcattattga                                                                                                              >  1:616701‑M1/1‑120 (MQ=255)                         tGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgatgcagtacgacactatcattattga                                                                                                              <  2:326413‑M1/149‑30 (MQ=255)                         gcaagacggagCCGGGCGGTTACATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgatgcagtacgacactatcattattga                                                                                                              <  2:30607‑M1/144‑30 (MQ=255)                                             tGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgatgcagtacgacactatcattattgacgaagcgcacgaacgcagcc                                                                                          <  2:405493‑M1/149‑50 (MQ=255)                                                aTCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgatgcagtacgacactatcattattgacgaagcgcacgaacgcagcctga                                                                                       >  1:101709‑M1/1‑97 (MQ=255)                                                 tCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgatgcagtacgacactatcattattgacgaagggcacgaacgcagcctgaa                                                                                      >  1:295357‑M1/1‑96 (MQ=255)                                                     ttACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgatgcagtacgacactatcattattgacgaagcgcacgaacgcagcctgaatatc                                                                                  <  1:202788‑M1/149‑58 (MQ=255)                                                                                aaGTGATAACCTGGTGGTCATGCTGATCACCGACGGTATCCTGCAGGCGGAGATCCAGCAAGACCctgatgcagtactacactatcattattgacgaagcgcacgaacgcagcctgaatatcgattttttgctcggctatttgaaagag                                                       <  1:362663‑M1/149‑85 (MQ=255)                                                                                   tGAGAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgatgcagtacgacactatcattattgacgaagcgcacgaacgcagcctgaatatcgattttttgctcggctatttgaaagagttg                                                    <  1:603117‑M1/149‑88 (MQ=255)                                                                                          acGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgatgcagtacgacactatcattattgacgaagcgcacgaacgcagcctgaatatcgattttttgctcggctatttgaaagagttgctgccgc                                             >  2:89199‑M1/1‑55 (MQ=255)                                                                                           cGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGGCCctgatgcagtacgaggctatgattattgacgaagcgcacgaacgcggcctgaatatcgattttttgctcggctttttggaggagttgctgccgag                                            >  2:149427‑M1/1‑54 (MQ=255)                                                                                                gTCAAGAGGATGACGGACGGTATCCTGCTGCCGGAGATCCAGCAAGCCCctgatgcagtacgacactatcattattgacgaagcgcaggaacgcagcttgaatatcgattttttgctcggctatttgaaagagttgctgccgcggcgtc                                       <  2:493170‑M1/149‑101 (MQ=255)                                                                                                 tCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgatgcagtacgacactatcattattgacgaagcgcacgaacgcagcctgaatatcgattttttgctcggctatttgaaagagttgctgccgcggcgtcc                                      >  1:387724‑M1/1‑48 (MQ=255)                                                                                                         atgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCctgatgccgtacgacactatcattattgacgaagggcacgaacgcagcctgaatatcgattttttgctcggctatttgaaagagttgctgccgcggcgtcctgacctaa                              >  1:256743‑M1/1‑40 (MQ=255)                                                                                                                  ggTATCCTGCTGGCGGAGATCCAGCAAGACCgtgatgcagtaggacactatgattattgacgaggcgcacgaacgcagcctggatatcgattttttgctaggctttttgaaagagttggtgcagcggcgttctgaccttaaaatcatta                     >  2:28487‑M1/1‑31 (MQ=255)                                                                                                                            tgGCGGAGATCCAGCAAGACCctgatgcagtacgacactatcattattgacgaagcgcacgaacgcagcctgaatatcgattttttgctcggctatttgaaagagttgctgccgcggcgtcctgacctaaaaatcattatcacttccgc           >  1:398094‑M1/1‑21 (MQ=255)                                                                                                                                ggAGATCCAGCAAGACCctgatgcagtacgacactatcattattgacgaagcgcacgaacgcagcctgaatatcgattttttgctcggctatttgaaagagttgctgccgcggcgtcctgacctaaaaatcattatcacttccgcgact       >  2:474631‑M1/1‑17 (MQ=255)                                                                                                                                   gaTCCAGCAAGACCctgatgcagtacgacactatcattattgacgaagcgcacgaacgcagcctgaatatcgattttttgctcggctatttgaaagagttgctgccgcggcgtcctgacctaaaaatcattatcacttccgcgactatc    <  2:391995‑M1/149‑136 (MQ=255)                                                                                                                                     tCCAGCAAGACCctgatgcagtacgacactatcattattgtcgaagcgcacgaacgcagcctgaatatcgattttttgctcggctatttgaaagagttgctgccgcggcgtcctgacctaaaaatcattatcacttccgcgactatcga  <  2:231893‑M1/149‑138 (MQ=255)                                                                                                                                                |                                                                                                                                          TGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGACT  >  NC_000913_3_pae_tpiA/1483473‑1483753

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