Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,811,302	T→C	G123G (GGT→GGC)	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,302	0	T	C	100.0%	166.7 / NA	48	G123G (GGT→GGC)	tcyP	cystine/sulfocysteine:cation symporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/24);  total (24/24)

GATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGTGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTATGTTGGTAAAGTCTCTGATCTGAGCGTTCCGCAGCTGGTCTT  >  NC_000913/1811167‑1811436                                                                                                                                        |                                                                                                                                       gATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGc                                                                                                                                    <  2:421358/140‑1 (MQ=255)  aTCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCt                                                                                                                                   <  1:698555/140‑1 (MQ=255)   tCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCt                                                                                                                                   >  2:428448/1‑139 (MQ=255)   tCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCt                                                                                                                                   <  1:428448/139‑1 (MQ=255)    cGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTgg                                                                                                                                 >  1:274908/1‑140 (MQ=255)            cGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAAc                                                                                                                         <  2:709733/140‑1 (MQ=255)             gTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAAcc                                                                                                                        <  2:259263/140‑1 (MQ=255)              ttAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCt                                                                                                                       >  1:83902/1‑140 (MQ=255)                aGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACAACCCTGATTCCGCCCCTGGCCGGCGGGCTGGCAACCACCCTgt                                                                                                                     >  2:913200/1‑140 (MQ=255)                        ccTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTgtttggtttg                                                                                                             <  2:846556/140‑1 (MQ=255)                            tATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTgtttggttgacgac                                                                                                         >  2:36287/1‑134 (MQ=255)                                cTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGa                                                                                                     <  1:12195/140‑1 (MQ=255)                                cTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGa                                                                                                     <  2:4175/140‑1 (MQ=255)                                cTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGa                                                                                                     <  1:427264/140‑1 (MQ=255)                                   aGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAgg                                                                                                  >  2:595141/1‑140 (MQ=255)                                    gcgcGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTGCCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGt                                                                                                 <  1:501218/140‑1 (MQ=255)                                    gcgcGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGt                                                                                                 <  2:83902/140‑1 (MQ=255)                                     cgcgGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGa                                                                                                     <  1:332791/135‑1 (MQ=255)                                     cgcgGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGa                                                                                                     >  2:332791/1‑135 (MQ=255)                                     cgcgGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTc                                                                                                >  2:154827/1‑140 (MQ=255)                                     cgcgGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTc                                                                                                <  2:321948/140‑1 (MQ=255)                                        ggTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTc                                                                                                <  2:444980/137‑1 (MQ=255)                                        ggTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTc                                                                                                >  1:444980/1‑137 (MQ=255)                                        ggTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTgg                                                                                             >  1:963287/1‑140 (MQ=255)                                        ggTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTgg                                                                                             >  2:598784/1‑140 (MQ=255)                                        ggTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTgg                                                                                             >  2:671664/1‑140 (MQ=255)                                        ggTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTgg                                                                                             >  1:735849/1‑140 (MQ=255)                                            gCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCa                                                                                         <  2:592787/140‑1 (MQ=255)                                                   tGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGggtggt                                                                                  >  2:826416/1‑140 (MQ=255)                                                    gCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGgtggtg                                                                                 >  1:82074/1‑140 (MQ=255)                                                         aCGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGaa                                                                            >  2:13812/1‑140 (MQ=255)                                                               ctcAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCa                                                                      >  2:963420/1‑140 (MQ=255)                                                                       ggCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTTATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGaa                                                              <  2:328440/140‑1 (MQ=255)                                                                       ggCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGaa                                                              >  2:580539/1‑140 (MQ=255)                                                                               cAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTg                                                      <  1:598784/140‑1 (MQ=255)                                                                                       tGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAAc                                              <  1:580539/140‑1 (MQ=255)                                                                                         aCCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTa                                            <  1:472034/140‑1 (MQ=255)                                                                                              cGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTATGTTg                                       >  1:232071/1‑140 (MQ=255)                                                                                                gTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTATGTTGGt                                     <  2:232071/140‑1 (MQ=255)                                                                                                        ttGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTATGTTGGTAAAGtctc                             <  1:595141/140‑1 (MQ=255)                                                                                                        ttGTTTACCACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTATGTTGGTAAAGtctc                             >  2:688245/1‑140 (MQ=255)                                                                                                                cACGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTATGTTGGTAAAGTCtctgatctga                     >  2:734505/1‑140 (MQ=255)                                                                                                                  cGCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTATGTTGGTAAAGTCTCTGATCTGAGc                   >  2:239814/1‑140 (MQ=255)                                                                                                                   gCTGATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTATGTTGGTAAAGTCTCTGATCTGAGCg                  <  1:13812/140‑1 (MQ=255)                                                                                                                      gATTGCGGCGCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTATGTTGGTAAAGTCTCTGATCTGAGCGTTc               <  2:963287/140‑1 (MQ=255)                                                                                                                            ggcgCTGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTATGTTGGTAAAGTCTCTGATCTGAGCGTTCCGCAGc         <  2:82074/140‑1 (MQ=255)                                                                                                                                 tGGTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTATGTTGGTAAAGTCTCTGATCTGAGCGTTCCGCAGCTGGTc    >  1:300436/1‑140 (MQ=255)                                                                                                                                   gTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTATGTTGGTAAAGTCTCTGATCTGAGCGTTCCGCAGCTGGTc    <  2:767315/138‑1 (MQ=255)                                                                                                                                   gTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTATGTTGGTAAAGTCTCTGATCTGAGCGTTCCGCAGCTGGTc    >  1:767315/1‑138 (MQ=255)                                                                                                                                   gTCGGCGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTATGTTGGTAAAGTCTCTGATCTGAGCGTTCCGCAGCTGGTCtt  >  1:119188/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       GATCGTTATGCCGTTAGTCTTCGCCTCTATTCTGAGCGCGGTTGCCCGTCTGCATAACGCATCTCAGTTAGGCAAAATCAGTTTTCTGACCATCGGTACGCTTTTGTTTACCACGCTGATTGCGGCGCTGGTCGGTGTGCTGGTCACCAACCTGTTTGGTTTGACGGCTGAAGGTCTGGTTCAGGGTGGTGCAGAAACTGCACGTCTGAACGCCATTGAAAGTAACTATGTTGGTAAAGTCTCTGATCTGAGCGTTCCGCAGCTGGTCTT  >  NC_000913/1811167‑1811436

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