Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,032,032	G→A	S222S (TCG→TCA)	uacT →	uric acid permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,032,032	0	G	A	100.0%	111.0 / NA	35	S222S (TCG→TCA)	uacT	uric acid permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/19);  total (16/19)

GCAATATGGTAATCCCGTTTATTTAGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCGTGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGCC  >  NC_000913/3031897‑3032152                                                                                                                                        |                                                                                                                         gCAATAGGGTAATCCCGTTTATTTAGGTATCTCTTTTACCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATgg                                                                                                                       <  2:671133/139‑1 (MQ=255)           aaTCCCGTTTATTTAGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGatgaatg                                                                                                             <  1:408162/139‑1 (MQ=255)             tCCCGTTTATTTAGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGAtgaatgaa                                                                                                           <  1:351734/139‑1 (MQ=255)               ccGTTTATTTAGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGt                                                                                                         >  2:499645/1‑139 (MQ=255)                    tatttaGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATg                                                                                                                        >  1:915348/1‑119 (MQ=255)                    tatttaGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATg                                                                                                                        <  2:915348/119‑1 (MQ=255)                            tATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCggg                                                                                            <  1:859231/139‑1 (MQ=255)                              tctcCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATgg                                                                                                                       >  1:32381/1‑110 (MQ=255)                              tctcCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATgg                                                                                                                       <  2:32381/110‑1 (MQ=255)                               ctcCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTa                                                                                         >  1:412048/1‑139 (MQ=255)                                tcCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTAc                                                                                        >  1:608808/1‑139 (MQ=255)                                 ccTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACa                                                                                       <  2:412048/139‑1 (MQ=255)                                           cTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCa                                                                             <  2:345789/139‑1 (MQ=255)                                                tttttATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTGCATGGtt                                                                        <  2:608808/139‑1 (MQ=255)                                                                 aCTCGCTATGCGAAAGGATTTATGTCCAACGTCGGCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGa                                                       <  1:152866/139‑1 (MQ=255)                                                                     gCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGa                                                       >  1:923/1‑135 (MQ=255)                                                                     gCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGa                                                       <  2:923/135‑1 (MQ=255)                                                                     gCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTc                                                   >  2:255946/1‑139 (MQ=255)                                                                        aTGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGtt                                                >  1:913360/1‑139 (MQ=255)                                                                        aTGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGtt                                                >  1:439373/1‑139 (MQ=255)                                                                          gCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTg                                              >  2:820952/1‑139 (MQ=255)                                                                                                  gCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGttt                      <  1:820952/139‑1 (MQ=255)                                                                                                              gggATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCggg                                                                                            >  1:810469/1‑57 (MQ=255)                                                                                                              gggATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCggg                                                                                            <  2:810469/57‑1 (MQ=255)                                                                                                              gggATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCa          >  1:204423/1‑139 (MQ=255)                                                                                                               ggATTGTATTTGGCTTTTTACTTTCATTGATGATGAATGAAGTCAATTTATCCGGGCTACATGATTCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCAt         <  2:809965/139‑1 (MQ=255)                                                                                                                   tGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCggg                                                                                            <  2:291095/52‑1 (MQ=255)                                                                                                                   tGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCggg                                                                                            >  1:291095/1‑52 (MQ=255)                                                                                                                   tGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATcc                           >  2:386016/1‑117 (MQ=255)                                                                                                                   tGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATcc                           <  1:386016/117‑1 (MQ=255)                                                                                                                     tATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCg                                                                                              >  2:332503/1‑48 (MQ=255)                                                                                                                     tATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCg                                                                                              <  1:332503/48‑1 (MQ=255)                                                                                                                      aTTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGcc  <  1:499645/139‑1 (MQ=255)                                                                                                                          ggCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGcc  <  2:378711/135‑1 (MQ=255)                                                                                                                          ggCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGcc  >  1:378711/1‑135 (MQ=255)                                                                                                                                        |                                                                                                                         GCAATATGGTAATCCCGTTTATTTAGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCGTGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGCC  >  NC_000913/3031897‑3032152

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