Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,372,139	T→G	S146R (AGT→CGT)	nanT ←	N‑acetylneuraminate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,372,139	0	T	G	100.0%	95.5 / NA	27	S146R (AGT→CGT)	nanT	N‑acetylneuraminate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (11/16);  total (11/16)

GCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCTGAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGT  >  NC_000913/3372015‑3372252                                                                                                                             |                                                                                                                  gCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGCCCGGAACCACCAGGCTATAGACCTGAGCGGCAACGACGGCGCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGt                                                                                                           <  2:107013/133‑1 (MQ=255) gCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCTGAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGt                                                                                                           <  2:438644/133‑1 (MQ=255)  cAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCTGAGCGGCAACGACGGCCCCCACAGAGATGCCTGAAATCATAAAACCACGGGCTTTGtt                                                                                                          >  1:205755/1‑133 (MQ=255)            tGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCTGAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGt                                                                                                >  1:535287/1‑133 (MQ=255)                    tGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCTGAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCa                                                                                        <  2:488554/133‑1 (MQ=255)                       cGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCTGAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCt                                                                                     <  2:412764/133‑1 (MQ=255)                         aTAAAGAACAACGCAGGCCAGCCCCAGAACGGAACCACCAGGCTATAGACCTGAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCttt                                                                                   <  2:205755/133‑1 (MQ=255)                                    cgcACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCTGAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAg                                                                        >  2:380136/1‑133 (MQ=255)                                        cGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCTGAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTgg                                                                    <  2:535287/133‑1 (MQ=255)                                         gCCAGCCCCAGACCGGAACCACCAGGCTATAGACCTGAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGgc                                                                   >  1:97847/1‑133 (MQ=255)                                                   gACCGGAACCACCAGGCTATAGACCTGAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCGTTCCATGACATAAGTGGTGCTGGAACCg                                                         >  2:170326/1‑133 (MQ=255)                                                                           cttaGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCg                                 <  2:445882/130‑1 (MQ=255)                                                                           cTGAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCg                                 >  1:488275/1‑133 (MQ=255)                                                                           cTGAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCg                                 <  1:341475/133‑1 (MQ=255)                                                                            tGAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGa                                >  2:516353/1‑133 (MQ=255)                                                                              aGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATg                              >  1:102259/1‑133 (MQ=255)                                                                                cGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTTACCCGCCAACCCCCTGCCGATGAc                            >  2:255274/1‑133 (MQ=255)                                                                                     acgtcgGCCCCCCCAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCTCCAGCCCCATGCCGATGACCAGAc                       <  2:82500/133‑1 (MQ=255)                                                                                       gacgGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGa                     <  2:105741/133‑1 (MQ=255)                                                                                          ggCCCCCACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCg                  <  2:97847/133‑1 (MQ=255)                                                                                            cccccACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGAt                <  1:516353/133‑1 (MQ=255)                                                                                            cccccACACAGAAGCTTGCAATTAAAAAACTACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGAt                <  1:170326/133‑1 (MQ=255)                                                                                               ccACAGAGAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATaaa             <  1:477325/133‑1 (MQ=255)                                                                                                   agagAAGCTTGAAATCAATAAACCACGGGCCTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATg         <  1:255274/133‑1 (MQ=255)                                                                                                   agagAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATg         <  1:495737/133‑1 (MQ=255)                                                                                                   agagAAGCCTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATg         <  1:380136/133‑1 (MQ=255)                                                                                                          ccTGAAATCAAAAAACCACGGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGt  >  2:19001/1‑133 (MQ=255)                                                                                                          ccTGAAATCAAAAAACCACGGGCTTTGTCACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGGGATGt  >  2:323625/1‑133 (MQ=255)                                                                                                                             |                                                                                                                  GCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGACCTGAGCGGCAACGACGGCCCCCACAGAGAAGCCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGT  >  NC_000913/3372015‑3372252

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