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%	122.5 / NA	38	S222S (TCG→TCA)	uacT	uric acid permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/19);  total (19/19)

TCCCGTTTATTTAGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCGTGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGCCGTGTTAATCAT  >  NC_000913/3031909‑3032163                                                                                                                            |                                                                                                                                    tCCCGTTTATTTAGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGAtgaatgaa                                                                                                                      >  2:359644/1‑139 (MQ=255)  cccGTTTATTTAGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAg                                                                                                                     <  2:37212/139‑1 (MQ=255)   ccGTTTATTTAGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGt                                                                                                                    >  2:298584/1‑139 (MQ=255)     gTTTATTTAGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTGTTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCa                                                                                                                  >  1:25557/1‑139 (MQ=255)      tttatttaGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCaa                                                                                                                 >  1:219695/1‑139 (MQ=255)            taGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATc                                                                                                           <  1:400642/139‑1 (MQ=255)            taGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATc                                                                                                           <  2:393718/139‑1 (MQ=255)              ggTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCg                                                                                                         <  2:52174/139‑1 (MQ=255)               gTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGatgaatg                                                                                                                        >  1:450311/1‑123 (MQ=255)               gTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGatgaatg                                                                                                                        <  2:450311/123‑1 (MQ=255)                  tctcCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCt                                                                                                     >  2:538320/1‑139 (MQ=255)                    tcCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTAc                                                                                                   >  2:229387/1‑139 (MQ=255)                       tttGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACatg                                                                                                <  2:678381/139‑1 (MQ=255)                          gCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCggg                                                                                                       <  2:442281/129‑1 (MQ=255)                          gCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCggg                                                                                                       >  1:442281/1‑129 (MQ=255)                               cTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCa                                                                                        <  2:568685/139‑1 (MQ=255)                                ttAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCAt                                                                                       <  2:351779/139‑1 (MQ=255)                                             tGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGt                                                                          >  2:169892/1‑139 (MQ=255)                                               cTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTTCATGATGCTTCATGGTTTGCGATTGTTa                                                                        <  2:495285/139‑1 (MQ=255)                                               cTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTa                                                                        >  2:384320/1‑139 (MQ=255)                                                     aCTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGa                                                                  <  2:585752/139‑1 (MQ=255)                                                            aTGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGatgaat                                                                                                                         >  1:650817/1‑77 (MQ=255)                                                            aTGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGatgaat                                                                                                                         <  2:650817/77‑1 (MQ=255)                                                              gCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTg                                                         >  2:709177/1‑139 (MQ=255)                                                                       tttATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGa                                                >  2:719432/1‑139 (MQ=255)                                                                               cAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGAt                                        >  1:7499/1‑139 (MQ=255)                                                                                  cGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATccc                                     >  2:277862/1‑139 (MQ=255)                                                                                                   ggATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTATACGATGCTTCATGGTTTGCGATTGTTTCGCCGATGTCGTTTGATATGCTGATTTTCGATCCCGTTTCCATTCTGACAAt                    <  1:298584/139‑1 (MQ=255)                                                                                                    gATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTCGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCGTTCTGACCATg                   >  2:521515/1‑139 (MQ=255)                                                                                                     aTTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGa                  <  1:373760/139‑1 (MQ=255)                                                                                                     aTTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGa                  <  1:777994/139‑1 (MQ=255)                                                                                                      ttGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCggg                                                                                                       <  1:748635/53‑1 (MQ=255)                                                                                                      ttGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCggg                                                                                                       >  2:748635/1‑53 (MQ=255)                                                                                                          aTTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGcc             <  2:283549/139‑1 (MQ=255)                                                                                                               gCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGCCGTGtt        >  2:251320/1‑139 (MQ=255)                                                                                                                   tttACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGCCGTGTTAatc    >  1:374134/1‑139 (MQ=255)                                                                                                                    ttACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGCCGTGTTAatca   <  2:363427/139‑1 (MQ=255)                                                                                                                     tACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGCCGTGTTAatcat  <  2:7499/139‑1 (MQ=255)                                                                                                                            |                                                                                                                                    TCCCGTTTATTTAGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCGTGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGCCGTGTTAATCAT  >  NC_000913/3031909‑3032163

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