Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,638,525	A→G	E295G (GAA→GGA)	pitA →	phosphate transporter, low‑affinity; tellurite importer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,638,525	0	A	G	100.0%	200.5 / NA	57	E295G (GAA→GGA)	pitA	phosphate transporter, low‑affinity; tellurite importer
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/25);  total (32/25)

CGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGAAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGACCACCGATGTGGA  >  NC_000913/3638431‑3638621                                                                                               |                                                                                                 cgtGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTg                                                                                             >  1:145105/1‑100 (MQ=255)   tGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGgc                                                                                           <  1:474700/100‑1 (MQ=255)    gATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGgcg                                                                                          <  1:733487/100‑1 (MQ=255)     aTGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGgcgc                                                                                         <  2:50971/100‑1 (MQ=255)      tGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGcgca                                                                                        >  2:907140/1‑100 (MQ=255)       gCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGcgcaa                                                                                       >  2:467413/1‑100 (MQ=255)         cATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGcgcaacg                                                                                     >  2:691356/1‑100 (MQ=255)          aTCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGcgcaacgc                                                                                    <  1:676003/100‑1 (MQ=255)                aacGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTg                                                                              >  1:508686/1‑100 (MQ=255)                   gTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCgg                                                                           <  2:840008/100‑1 (MQ=255)                     cGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGTAg                                                                         <  2:491429/100‑1 (MQ=255)                          cTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCa                                                                    >  2:424364/1‑100 (MQ=255)                           ttACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCAc                                                                   <  2:836925/100‑1 (MQ=255)                           ttACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCAc                                                                   >  2:128025/1‑100 (MQ=255)                            tACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGATGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACt                                                                  <  2:834784/100‑1 (MQ=255)                              cTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGc                                                                <  1:267953/100‑1 (MQ=255)                               tttGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGcc                                                               >  1:1266731/1‑100 (MQ=255)                                ttGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCa                                                              >  2:1264112/1‑100 (MQ=255)                                ttGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCa                                                              >  1:575464/1‑100 (MQ=255)                                    gcagcaTCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCg                                                          <  1:760669/100‑1 (MQ=255)                                        caTCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGa                                                      <  2:824911/100‑1 (MQ=255)                                             tGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATAcc                                                 <  2:682106/100‑1 (MQ=255)                                               cgcTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCAt                                               >  1:796856/1‑100 (MQ=255)                                                 ctgctCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTa                                             >  1:1012265/1‑100 (MQ=255)                                                 ctgctCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTa                                             >  1:161279/1‑100 (MQ=255)                                                      cAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAAcgcg                                        >  2:211300/1‑100 (MQ=255)                                                       aaaCAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACgcgc                                       <  2:1130395/100‑1 (MQ=255)                                                        aaCAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACgcgc                                       >  2:597324/1‑99 (MQ=255)                                                            ggCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAAc                                  >  1:660399/1‑100 (MQ=255)                                                                  cGGTGCTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCATCCGcc                              <  1:895883/98‑1 (MQ=255)                                                                       cTGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAgg                       <  2:345062/100‑1 (MQ=255)                                                                        tGATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGt                      >  1:687489/1‑100 (MQ=255)                                                                         gATCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTa                     >  1:525097/1‑100 (MQ=255)                                                                          aTCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTAt                    >  1:447070/1‑100 (MQ=255)                                                                           tCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATg                   <  1:831926/100‑1 (MQ=255)                                                                           tCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATg                   >  1:212794/1‑100 (MQ=255)                                                                           tCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATg                   <  2:894147/100‑1 (MQ=255)                                                                           tCAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATg                   <  1:1188784/100‑1 (MQ=255)                                                                            cAGTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGt                  >  2:1185537/1‑100 (MQ=255)                                                                              gTTAGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTg                <  2:87090/100‑1 (MQ=255)                                                                                 aGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGacc             >  2:538489/1‑100 (MQ=255)                                                                                 aGTACCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGacc             <  1:907147/100‑1 (MQ=255)                                                                                    aCCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGaccacc          >  2:218688/1‑100 (MQ=255)                                                                                    aCCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGaccacc          <  2:781913/100‑1 (MQ=255)                                                                                    aCCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGaccacc          >  1:573928/1‑100 (MQ=255)                                                                                    aCCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGaccacc          <  1:34539/100‑1 (MQ=255)                                                                                    aCCGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGaccacc          >  1:158759/1‑100 (MQ=255)                                                                                     ccGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGACCACCg         >  2:176581/1‑100 (MQ=255)                                                                                     ccGGCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGACCACCg         >  1:183466/1‑100 (MQ=255)                                                                                     ccGGCTCCGGGAGCTCGCGCAACGCAACCGGCGGAGTGCGTCTTCCATCCGTCGAATTCCATTAACGCGCTCAGCCGTCTGAAATGTATGTTGaccacc          <  2:587399/99‑1 (MQ=255)                                                                                        gCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGACCACCGAtg      <  2:1255358/100‑1 (MQ=255)                                                                                        gCTCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGACCACCGAtg      >  1:1041044/1‑100 (MQ=255)                                                                                          tCCGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGACCACCGAtgtg    <  2:1266720/100‑1 (MQ=255)                                                                                           ccGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGACCACCGATGTgg   >  1:695358/1‑100 (MQ=255)                                                                                           ccGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGACCACCGATGTgg   >  1:557489/1‑100 (MQ=255)                                                                                            cGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGACCACCGATGTGGa  >  2:412014/1‑100 (MQ=255)                                                                                            cGGGAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGACCACCGATGTGGa  >  1:1180972/1‑100 (MQ=255)                                                                                               |                                                                                                 CGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGAAGCTGGCGCAACGCAACCTGCGGAGTTCCACTGCCATCCGTCGAATACCATTAACGCGCTCAACCGCCTGAAAGGTATGTTGACCACCGATGTGGA  >  NC_000913/3638431‑3638621

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