Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,817,197	G→A	R103C (CGC→TGC)	yqaB ←	fructose‑1‑phosphate phosphatase YqaB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,817,197	0	G	A	100.0%	66.2 / NA	20	R103C (CGC→TGC)	yqaB	fructose‑1‑phosphate phosphatase YqaB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/11);  total (9/11)

CCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGC  >  NC_000913/2816954‑2817421                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                 ccACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGAcc                                                                                                                                                                                                                            <  1:60831/250‑1 (MQ=255)        gTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAAc                                                                                                                                                                                                                    >  2:57072/1‑251 (MQ=255)                  gctgcACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTcaccac                                                                                                                                                                                                          <  2:126694/251‑1 (MQ=255)                  gctgcACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTcacc                                                                                                                                                                                                            <  1:1791/249‑1 (MQ=255)                     gcACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCAGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACAt                                                                                                                                                                                                        <  2:168651/250‑1 (MQ=255)                                          aaaTCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGAc                                                                                                                                                                                  <  2:238863/251‑1 (MQ=255)                                                gCATCTTCAAAGACCACACACTGCGTCGGTTGAACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGAGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATc                                                                                                                                                                            >  1:207092/1‑251 (MQ=255)                                                     ttCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATAAACAAGAGGAAGCGGGTCGACGCTATCcagc                                                                                                                                                                        >  1:64661/1‑250 (MQ=255)                                                               acacacTGCGTCGGTAGCACGCCCATGCGCTGCGCGCACAAAAAAAATGTGTCTTGCTCGGGTTTATGGGGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGCAATCCCAGGTGCACCAGCAAAGCCTCAGCGATGGCGCTTACACTCCCCGTTCCGACAGCCATGGGGCAACGCCCATGCCACCTTTTCACCACATCCACAAGAGGAAGAGGTTCGACGCTAGCCAGCAGCATActtct                                                                                                                                                             <  2:166768/251‑1 (MQ=255)                                                                   acTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACt                                                                                                                                                         >  1:15413/1‑251 (MQ=255)                                                                             ttGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGttt                                                                                                                                               >  2:90331/1‑251 (MQ=255)                                                                                    cccATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCAcg                                                                                                                                        <  2:62130/251‑1 (MQ=255)                                                                                                  gcACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCg                                                                                                                          >  1:170389/1‑251 (MQ=255)                                                                                                                                          aCGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGc                                                                                  <  2:5293/251‑1 (MQ=255)                                                                                                                                                      gCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGt                                                                      >  1:70700/1‑251 (MQ=255)                                                                                                                                                                                tAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCg                                            >  2:182574/1‑251 (MQ=255)                                                                                                                                                                                             gacAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCAATTGGGCAACGACCATGCAAACTTTTCACCACATCAACAAGAGGAATCAGTTCGACGCTATCCAGCAGAATACTTATTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTAGGCGATCCATTAAGCGAAATCATCGACTGAATATCGt                                >  2:9513/3‑250 (MQ=255)                                                                                                                                                                                               cAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATGAAGCGCAATCATCGCCTGAATATCGTACt                             <  1:186756/251‑1 (MQ=255)                                                                                                                                                                                                                       ttCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATAc       <  1:220501/249‑1 (MQ=255)                                                                                                                                                                                                                          aCTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCGAATACTTcgc  <  1:207823/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                 CCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGC  >  NC_000913/2816954‑2817421

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