Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,745,917	C→T	P372S (CCG→TCG)	yiaN →	2,3‑diketo‑L‑gulonate:Na(+) symporter ‑ membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,745,917	0	C	T	100.0%	73.2 / NA	24	P372S (CCG→TCG)	yiaN	2,3‑diketo‑L‑gulonate:Na(+) symporter ‑ membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/11);  total (13/11)

GCTGATTACCATTGCTGAACTGCCGATGATGGTTTCTGATTTACTGCAACCGTTGGTCGATTCACCGCGTCTGCTGTTTATCGTCATTATGGTGGCTATCCTGATTGTCGGCATGGTCATGGATTTAACGCCGACGGTATTAATTCTTACCCCGGTATTAATGCCTTTAGTTAAAGAGGCAGGAATCGATCCGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTCCCGATCTCATCATCCTCCCTTTGAAATGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTG  >  NC_000913/3745673‑3746162                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       gCTGATTACCATTGCTGAACTGCCGATGATGGTTTCTGATTTACTGCAACCGTTGGTCGATTCACCGCGTCTGCTGTTTATCGTCATTATGGGGGCTATCCTGATTGTCGGCATGGTCATGGATTTAACGCCGACGGTATTAATTCTTACCCCGGTATTAATGCCTTTAGTTAAAGAGGCAGGAATCGATCCGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTACATCGCCt                                                                                                                                                                                                                                                   >  2:396240/1‑250 (MQ=255)    gATTACCATTGCTGAACTGCCGATGATGGTTTCTGATTTACTGCAACCGTTGGTCGATTCACCGCGTCTGCTGTTTATCGTCATTATGGTGGCTATCCTGATTGTCGGCATGGTCATGGATTTAACGCCGACGGTATTAATTCTTACCCCGGTATTAATGCCTTTAGTTAAAGAGGCAGGAATCGATCCGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTACATCGCCTATc                                                                                                                                                                                                                                                <  2:290015/250‑1 (MQ=255)                 gAACTGCCGATGATGGTTTCTGATTTACTGCAACCGTTGGTCGATTCACCGCGTCTGCTTTTTATCGTCATTATGGTGGCTATCCTGATTGTCGGCATGGTCATGGATTTAACGCCGACGGTATTAATTCTTACCCCGGTATTAATGCCTTTAGTTAAAGAGGCAGGAATCGATCCGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTTCATCGCCTATCGGTAATGTACTTaa                                                                                                                                                                                                                                  <  2:260797/251‑1 (MQ=255)                                                       ggTCGATTCACCGCGTCTGCTGTTTATCGTCATTATGGTGGCTATCCTGATTGTCGGCATGGTCATGGATTTAACGCCGACGGTATTAATTCTTACCCCGGTATTAATGCCTTTAGTTAAAGAGGCAGGAATCGATCCGATTTATTTCGGTGTCATGTTTCTCATTAACTGCTCAATCGGATTAATTACATCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑CAAACTCAAATTCGatgatg                                                                                                                                                                                           >  1:133714/1‑251 (MQ=255)                                                           gATTCACCGCGTCTGCTGTTTATCGTCATTATGGTGGCTATCCTGATTGTCGGCATGGTCATGGATTTAACGCCGACGGTATTAATTCTTACCCCGGTATTAATGCCTTTAGTTAAAGAGGCAGGAATCGATCCGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTACATCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAATTCGATGATGCAGt                                                                                                                                                                                       <  1:1882/251‑1 (MQ=255)                                                                          ctgTTTATCGTCATTATGGTGGCTATCCTGATTGTCGGCATGGTCATGGATTTAACGCCGACGGTATTAATTCTTACCCCGGTATTAATGCCTTTCGTTAAAGAGGCAGGAATCGATCCGATTTATTTCGGTGCCATGTTTATCATTAACTGCTCAATAGGATTAATTACATCGCCTATCGGTAATGGACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAATTCGATGCTGCAGTCAGAGGCGTTTTccc                                                                                                                                                                        >  2:112483/1‑251 (MQ=255)                                                                                                             cGGCATGGTCATGGATTTAACGCCGACGGTATTAATTCTTACCCCGGTATTAATGCCTTTAGTTAAAGAGGCAGGAATCGATCCGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTACATCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTAt                                                                                                                                     >  2:94545/1‑251 (MQ=255)                                                                                                             cGGCATGGTCATGGATTTAACGCCGACGGTATTAATTCTTACCCCGGTATTAATGCCTTTAGTTAAAGAGGCAGGAATCGATCCGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTACATAGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGTTTTTATATTCATTATTAGTCGTAt                                                                                                                                     >  2:138182/1‑251 (MQ=255)                                                                                                                      cATGGATTTAACGCCGACGGTATTAATTCTTACCCCGGTATTAATGCCTTTAGTTAAAGAGTCAGGAATCGATCTGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTACATCGCGTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAGTTCGATGATGCAGTCAGAGTCGTGTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTa                                                                                                                            <  1:398275/251‑1 (MQ=255)                                                                                                                                          tATTAATTCTTACCCCGGTATTAATGCCTTTAGTTAAAGAGGAAGGAATCGATCCGATTTATTTCGGTGTAATGTTTATCTTTAACTGCTCAATCGACTTAATTACATCGCCTATTGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCTTTTTTAGTCGTATTTGTTTTTATTCCCGATCTCATCAtcctc                                                                                                        >  2:44672/1‑251 (MQ=255)                                                                                                                                                    tACCCCGGTATTAATGCCTTTAGTTAAAGAGTCAGGAATCCATCCGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTACATCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAAGTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCGGGTTTTATATTCATTATTAGTCGTATTTGTTGTTATTCCCGATCGCATCATCCTGCCTTTGaa                                                                                                <  2:261923/249‑1 (MQ=255)                                                                                                                                                          ggTATTAATGCCTTTAGTTAAAGAGGCAGGAATCGATCCGATTTATTTCGGTGTCATGTTCATCATTAACTGCTCAATCGGCTTAATTACATCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTCCCGATCTCATCATCCTCCCTTTGAAATGGatt                                                                                         >  1:351965/1‑250 (MQ=255)                                                                                                                                                                   gCCTTTAGTTAAAGAGGCAGGAATCGATCCGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTACATCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTCCCGATCTCATCATCCTCCCTTTGAAATGGATTAATTAAAAgg                                                                               <  1:113172/251‑1 (MQ=255)                                                                                                                                                                         agttacggagccAGGAATCGGTCCGATTTATTTCGGCGTCGTGTTTATCAATAACCGCTCAATCCGGTTACGTACATCGCCTATCGGTAATGTACTTAAGGTTATTTCATGG‑TGGCCAAAACTCAAATTCGCTGATGCAGTCATTGGCGTTTTACCTTACGTCCTGGTTTTGTGTTCATTCTTAGTCGCATTTGTTTTTATTCCCGATCTCATCATCCTCCCTTTGACATGGATTAATTAATAGGAAAatc                                                                         >  2:137870/12‑250 (MQ=255)                                                                                                                                                                                 agGCAGGAATCGATTCGATTTATTGCGGTGTCATGTTTATCATTAACTGCACAATCGGCTTAATTACATAGCTTATCGGTAATGTACTTAACGTTATTTCAGGGGTGTC‑AAAACTCCAATTCGATGATGCATTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTGTTATTCCCGATCTCACCATCCTCTCTTGGAAATGGATTAATTAAAAGGAAAATATTATGaaa                                                                 <  1:152185/251‑1 (MQ=255)                                                                                                                                                                                 agGCAGGAATCGATCCGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTACATCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAGATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTCCCGATCTCATCATCCTCCCTTTGAAATGGATTAATTAAAAGGAAAATATTATGaa                                                                  >  2:112572/1‑250 (MQ=255)                                                                                                                                                                                         aTCGATCCGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTACATCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTCCCGATCTCATCATCCTCCCTTTg                                                                                                  <  1:279046/210‑1 (MQ=255)                                                                                                                                                                                         aTCGATCCGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTACATCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTCCCGATCTCATCATCCTCCCTTTg                                                                                                  >  2:279046/1‑210 (MQ=255)                                                                                                                                                                                           cGATCCGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTACATCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTCCCGATCTCATCATCCTCCCTTTGAAATGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTg                                                       <  1:5126/251‑1 (MQ=255)                                                                                                                                                                                            gttCCGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATAGGCTTAATTACATCTCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAAATCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTCCCGATCTCATCATGCTCCCTTGGAAATGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGt                                                      <  1:359635/249‑1 (MQ=255)                                                                                                                                                                                               ccGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTACATCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTCCCGATCTCATCATCCTCCCTTTGAAATGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAAc                                                   >  2:137039/1‑251 (MQ=255)                                                                                                                                                                                                                          aTTAACTGCTCAATCGGCTTAATTACATCGCCTATCGGTAATGTACTTAACGTTATTTCAGGTGTGGC‑AAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTCCCGATCTCATCATCCTCCCTTTGAAATGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGctggccgggc                        <  2:96405/251‑1 (MQ=255)                                                                                                                                                                                                                               cTGCTCAATCGGCTTAATTACATCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTCCCGATCTCATCATCCTCCCTTTGAAATGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCAt                   >  2:324019/1‑251 (MQ=255)                                                                                                                                                                                                                                                ttACATCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTCCCGATCTCATCATCCTCCCTTTGAAATGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTg  >  1:161065/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       GCTGATTACCATTGCTGAACTGCCGATGATGGTTTCTGATTTACTGCAACCGTTGGTCGATTCACCGCGTCTGCTGTTTATCGTCATTATGGTGGCTATCCTGATTGTCGGCATGGTCATGGATTTAACGCCGACGGTATTAATTCTTACCCCGGTATTAATGCCTTTAGTTAAAGAGGCAGGAATCGATCCGATTTATTTCGGTGTCATGTTTATCATTAACTGCTCAATCGGCTTAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGC‑AAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTCCCGATCTCATCATCCTCCCTTTGAAATGGATTAATTAAAAGGAAAATATTATGAAATTACGCTCTGTAACCTACGCATTATTCATTGCTGGCCTGGCTGCATTCAGCACATCTTCTCTG  >  NC_000913/3745673‑3746162

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