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

GTTTATTTAGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCGTGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGCCGTGTTAATCAT  >  NC_000913/3031913‑3032163                                                                                                                        |                                                                                                                                    gTTTATTTAGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCa                                                                                                                  >  2:372584/1‑139 (MQ=255)  tttatttaGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCaa                                                                                                                 >  2:446916/1‑139 (MQ=255)       ttaGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTAt                                                                                                            <  1:372584/139‑1 (MQ=255)            tATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCggg                                                                                                       <  2:472992/139‑1 (MQ=255)                   tttGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACatg                                                                                                <  1:446916/139‑1 (MQ=255)                    ttGCCGTCTTAATTTTTATCTTGCTCATTACTCACTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACatga                                                                                               >  1:801945/1‑139 (MQ=255)                            ttAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCAt                                                                                       >  1:38536/1‑139 (MQ=255)                                          gCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGTGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACatg                                                                                                <  2:417929/116‑1 (MQ=255)                                          gCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACatg                                                                                                >  1:417929/1‑116 (MQ=255)                                           cTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGgat                                                                                                                                <  1:229678/83‑1 (MQ=255)                                           cTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGgat                                                                                                                                >  2:229678/1‑83 (MQ=255)                                                 aCTCACTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGa                                                                  <  2:801945/139‑1 (MQ=255)                                                     gCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGAtgaatgaa                                                                                                                      <  1:23293/83‑1 (MQ=255)                                                     gCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTa                                                                                                    <  1:884353/101‑1 (MQ=255)                                                     gCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTa                                                                                                    >  2:884353/1‑101 (MQ=255)                                                     gCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTCTTTACTTTCATGGATGAtgaatgaa                                                                                                                      >  2:23293/1‑83 (MQ=255)                                                        aTGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGtt                                                           >  1:511942/1‑139 (MQ=255)                                                        aTGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGtt                                                           >  2:928127/1‑139 (MQ=255)                                                             aaaGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCggg                                                                                                       >  1:830738/1‑90 (MQ=255)                                                             aaaGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCggg                                                                                                       <  2:830738/90‑1 (MQ=255)                                                                   tttATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGa                                                >  1:501698/1‑139 (MQ=255)                                                                     tATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATgg                                                                                                                                  >  1:625739/1‑55 (MQ=255)                                                                             aCGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGAt                                        >  2:375490/1‑137 (MQ=255)                                                                             aCGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGAt                                        <  1:375490/137‑1 (MQ=255)                                                                              cGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGtt                                                           <  1:185563/117‑1 (MQ=255)                                                                              cGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGtt                                                           >  2:185563/1‑117 (MQ=255)                                                                                               ggATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCAt                    <  2:501698/139‑1 (MQ=255)                                                                                                gATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCg                                                                                                         <  2:624955/53‑1 (MQ=255)                                                                                                gATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCg                                                                                                         >  1:624955/1‑53 (MQ=255)                                                                                                gATTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATg                   >  1:952724/1‑139 (MQ=255)                                                                                                 aTTGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGa                  <  1:928127/139‑1 (MQ=255)                                                                                                   tGTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACt                >  1:1003207/1‑139 (MQ=255)                                                                                                    gTATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGTCTg               <  1:521630/139‑1 (MQ=255)                                                                                                     tATTTGGCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGc              >  1:955380/1‑139 (MQ=255)                                                                                                          ggCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATc                                                                                                           >  1:725332/1‑41 (MQ=255)                                                                                                          ggCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATc                                                                                                           <  2:725332/41‑1 (MQ=255)                                                                                                          ggCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCg                                                                   >  2:799070/1‑81 (MQ=255)                                                                                                          ggCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCg                                                                   <  1:799070/81‑1 (MQ=255)                                                                                                           gCTTTTTACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGCCGTGtt        >  1:500345/1‑139 (MQ=255)                                                                                                               tttACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCa                                                                                        >  1:361786/1‑55 (MQ=255)                                                                                                               tttACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCa                                                                                        <  2:361786/55‑1 (MQ=255)                                                                                                               tttACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGCCGTGTTAatc    >  1:76445/1‑139 (MQ=255)                                                                                                                 tACTTTCATGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGCCGTGTTAatcat  <  2:952724/139‑1 (MQ=255)                                                                                                                        |                                                                                                                                    GTTTATTTAGGTATCTCCTTTGCCGTCTTAATTTTTATCTTGCTCATTACTCGCTATGCGAAAGGATTTATGTCCAACGTCGCCGTATTACTGGGGATTGTATTTGGCTTTTTACTTTCGTGGATGATGAATGAAGTCAATTTATCCGGGCTACATGATGCTTCATGGTTTGCGATTGTTACGCCGATGTCGTTTGGTATGCCGATTTTCGATCCCGTTTCCATTCTGACCATGACTGCCGTGTTAATCAT  >  NC_000913/3031913‑3032163

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