Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,414,699	T→C	intergenic (+27/‑48)	ackA → / → pta	acetate kinase A and propionate kinase 2/phosphate acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,414,699	0	T	C	100.0%	140.4 / NA	39	intergenic (+27/‑48)	ackA/pta	acetate kinase A and propionate kinase 2/phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/21);  total (18/21)

CAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCG  >  NC_000913/2414570‑2414786                                                                                                                                  |                                                                                        cAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTg                                                                                >  1:180349/1‑139 (MQ=255)    aTCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTaa                                                                             >  1:361837/1‑139 (MQ=255)                       gaagGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTaa                                                          <  2:180349/139‑1 (MQ=255)                       gaagGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTaa                                                          <  1:420077/139‑1 (MQ=255)                           gTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGaa                                                      >  2:64452/1‑139 (MQ=255)                           gTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGaa                                                      <  1:64452/139‑1 (MQ=255)                                       cGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATcggc                                                               <  2:567689/118‑1 (MQ=255)                                       cGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATcggc                                                               >  1:567689/1‑118 (MQ=255)                                        ggtggtTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCG‑CAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCgt                                        <  1:746946/139‑1 (MQ=255)                                            gtTATCCCAACCAACGGAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTc                                     <  1:857611/139‑1 (MQ=255)                                            gtTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTc                                     >  2:426600/1‑139 (MQ=255)                                            gtTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTc                                     <  1:305181/139‑1 (MQ=255)                                            gtTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTc                                     <  2:1111752/139‑1 (MQ=255)                                               aTCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCg                                  <  1:426600/139‑1 (MQ=255)                                               aTCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCg                                  >  1:22798/1‑139 (MQ=255)                                                    aaccaacGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATcggc                                                               <  1:580331/105‑1 (MQ=255)                                                    aaccaacGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATcggc                                                               >  2:580331/1‑105 (MQ=255)                                                     accaacGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATcggc                                                               >  1:613899/1‑104 (MQ=255)                                                     accaacGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATcggc                                                               <  2:613899/104‑1 (MQ=255)                                                      ccaacGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGtattatt                           <  2:149438/139‑1 (MQ=255)                                                      ccaacGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGtattatt                           >  1:1020023/1‑139 (MQ=255)                                                      ccaacGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGtattatt                           >  1:149438/1‑139 (MQ=255)                                                          cGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCg                                  <  2:838039/128‑1 (MQ=255)                                                          cGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCg                                  >  1:838039/1‑128 (MQ=255)                                                              gaaCTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGAt                   <  2:22798/139‑1 (MQ=255)                                                                         cgcgCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATcggc                                                               >  1:464861/1‑84 (MQ=255)                                                                         cgcgCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATcggc                                                               <  2:464861/84‑1 (MQ=255)                                                                           cgcAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCt                                                                                      >  2:1043862/1‑59 (MQ=255)                                                                           cgcAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCt                                                                                      <  1:1043862/59‑1 (MQ=255)                                                                                           ccTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTc                                     <  1:343101/92‑1 (MQ=255)                                                                                           ccTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTc                                     >  2:343101/1‑92 (MQ=255)                                                                                                                  ccAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGaaa                                                     <  2:849684/53‑1 (MQ=255)                                                                                                                  ccAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGaaa                                                     >  1:849684/1‑53 (MQ=255)                                                                                                                  ccAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGc                       >  2:695631/1‑83 (MQ=255)                                                                                                                  ccAGCTCAGATGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGc                       <  1:695631/83‑1 (MQ=255)                                                                                                                          gctGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGAt                   <  2:1005037/79‑1 (MQ=255)                                                                                                                          gctGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGAt                   >  1:1005037/1‑79 (MQ=255)                                                                                                                            tGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCg  <  1:882979/94‑1 (MQ=255)                                                                                                                            tGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCg  >  2:882979/1‑94 (MQ=255)                                                                                                                                  |                                                                                        CAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCG  >  NC_000913/2414570‑2414786

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑