Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,811,208	T→C	V92A (GTT→GCT)	tcyP →	cystine/sulfocysteine:cation symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,811,208	0	T	C	100.0%	260.6 / NA	66	V92A (GTT→GCT)	tcyP	cystine/sulfocysteine:cation symporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/33);  total (33/33)

GCAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCT  >  NC_000913/1811161‑1811253                                                |                                              gcAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTg                                             >  1:1987681/1‑50 (MQ=255)  cAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGc                                            >  1:2597120/1‑50 (MQ=255)  cAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGc                                            >  1:1730820/1‑50 (MQ=255)   aaaTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGcc                                           >  1:575846/1‑50 (MQ=255)   aaaTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGcc                                           >  1:1330065/1‑50 (MQ=255)   aaaTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGcc                                           <  1:2297446/50‑1 (MQ=255)    aaTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGccc                                          <  1:3516258/50‑1 (MQ=255)     aTGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCg                                         >  1:803936/1‑50 (MQ=255)      tGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGt                                        <  1:2803050/50‑1 (MQ=255)      tGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGt                                        >  1:2800228/1‑50 (MQ=255)       gATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTc                                       >  1:2587363/1‑50 (MQ=255)       gATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTc                                       >  1:910464/1‑50 (MQ=255)        aTCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCt                                      >  1:3871915/1‑50 (MQ=255)        aTCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCt                                      <  1:3671084/50‑1 (MQ=255)         tCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTg                                     <  1:850547/50‑1 (MQ=255)          cGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTg                                     >  1:518981/1‑49 (MQ=255)          cGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTg                                     >  1:720166/1‑49 (MQ=255)          cGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTg                                     >  1:1507459/1‑49 (MQ=255)           gTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGc                                    <  1:1532594/49‑1 (MQ=255)            ttATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCa                                   >  1:1916708/1‑49 (MQ=255)            ttATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCa                                   >  1:2608978/1‑49 (MQ=255)            ttATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCa                                   >  1:3029112/1‑49 (MQ=255)            ttATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCAt                                  <  1:2575719/50‑1 (MQ=255)             tATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCAt                                  >  1:3099884/1‑49 (MQ=255)               tGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAAc                               <  1:2839499/50‑1 (MQ=255)               tGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAAc                               <  1:408496/50‑1 (MQ=255)               tGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAAc                               <  1:3941780/50‑1 (MQ=255)               tGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAAc                               <  1:1503955/50‑1 (MQ=255)                 ccGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGc                             <  1:32406/50‑1 (MQ=255)                  cGTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCa                            >  1:305788/1‑50 (MQ=255)                   gTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAt                           >  1:3089293/1‑50 (MQ=255)                   gTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAt                           <  1:789924/50‑1 (MQ=255)                   gTTAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAt                           >  1:2757851/1‑50 (MQ=255)                    ttAGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAtc                          >  1:403761/1‑50 (MQ=255)                      aGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAtct                         >  1:1789629/1‑49 (MQ=255)                      aGTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAtct                         <  1:397764/49‑1 (MQ=255)                       gTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCAtctc                        <  1:3622542/49‑1 (MQ=255)                       gTCTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCa                       <  1:891544/50‑1 (MQ=255)                         cTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGt                     >  1:3450067/1‑50 (MQ=255)                         cTTCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGt                     >  1:1605946/1‑50 (MQ=255)                          ttCGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGt                     >  1:3432715/1‑49 (MQ=255)                            cGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGt                     <  1:1359887/47‑1 (MQ=255)                            cGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTa                   <  1:1823085/49‑1 (MQ=255)                            cGCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTa                   <  1:2946053/49‑1 (MQ=255)                             gCCTCTATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAgg                 >  1:2693537/1‑50 (MQ=255)                               ctctATTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCa               <  1:1316553/50‑1 (MQ=255)                                   aTTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAAt           >  1:1034732/1‑50 (MQ=255)                                   aTTCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAAt           >  1:131883/1‑50 (MQ=255)                                    ttCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCaaa             <  1:2078002/47‑1 (MQ=255)                                    ttCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCaaa             <  1:3305430/47‑1 (MQ=255)                                    ttCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATc          <  1:2533348/50‑1 (MQ=255)                                    ttCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATc          <  1:764918/50‑1 (MQ=255)                                    ttCTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATc          <  1:3765815/50‑1 (MQ=255)                                      cTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCa         <  1:1270182/49‑1 (MQ=255)                                      cTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCa         <  1:1605633/49‑1 (MQ=255)                                      cTGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCa         <  1:161778/49‑1 (MQ=255)                                       tGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGt       >  1:945623/1‑50 (MQ=255)                                       tGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGt       >  1:3781418/1‑50 (MQ=255)                                       tGAGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGt       >  1:3678989/1‑50 (MQ=255)                                         aGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAg        <  1:3015230/47‑1 (MQ=255)                                         aGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGtt      <  1:2263806/49‑1 (MQ=255)                                         aGCGCGGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGt       <  1:41305/48‑1 (MQ=255)                                           cgcgGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTc   <  1:3027452/50‑1 (MQ=255)                                           cgcgGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTc   <  1:1295370/50‑1 (MQ=255)                                           cgcgGCTGCCCGTCTGCATAACGCATCTCAGTTAGACAAAATCAGTTTTc   >  1:2142685/1‑50 (MQ=255)                                            gcgGCTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCt  >  1:3535482/1‑50 (MQ=255)                                                |                                              GCAAATGATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCT  >  NC_000913/1811161‑1811253

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