Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,637,751	A→C	N37T (AAC→ACC)	pitA →	metal phosphate:H(+) symporter PitA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,637,751	0	A	C	100.0%	206.8 / NA	57	N37T (AAC→ACC)	pitA	metal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (32/25);  total (32/25)

CATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTATGCCATTGTGCATATGCTGCCGA  >  NC_000913/3637615‑3637885                                                                                                                                         |                                                                                                                                       cATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGc                                                                                                                                     >  2:17950/1‑140 (MQ=255)       cTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGc                                                                                                                               <  2:74957/140‑1 (MQ=255)         cAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCaa                                                                                                                             <  1:112559/140‑1 (MQ=255)            aaTGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCg                                                                                                                          <  1:17950/140‑1 (MQ=255)            aaTGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCg                                                                                                                          <  1:53000/140‑1 (MQ=255)                 cGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATc                                                                                                                     >  2:881528/1‑140 (MQ=255)                 cGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATc                                                                                                                     >  2:808599/1‑140 (MQ=255)                  gTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCt                                                                                                                    >  2:285797/1‑140 (MQ=255)                      cGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATAc                                                                                                                >  2:127625/1‑140 (MQ=255)                          cTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCcgc                                                                                                            <  1:127625/140‑1 (MQ=255)                            aTGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCcgcgc                                                                                                          <  1:412939/140‑1 (MQ=255)                                     ttgtttgCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTtc                                                                                                 >  1:708565/1‑140 (MQ=255)                                      tgtttgCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTtct                                                                                                <  1:931083/140‑1 (MQ=255)                                       gtttgCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTtctc                                                                                               >  2:688385/1‑140 (MQ=255)                                        tttgCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCa                                                                                              >  2:170706/1‑140 (MQ=255)                                         ttgCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAg                                                                                             <  1:881528/140‑1 (MQ=255)                                            cTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTc                                                                                          >  2:811431/1‑140 (MQ=255)                                                    gATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGt                                                                                  >  1:283746/1‑140 (MQ=255)                                                    gATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGt                                                                                  >  2:230408/1‑140 (MQ=255)                                                              ccGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGgcggcg                                                                        <  1:774495/140‑1 (MQ=255)                                                                   cTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTAtt                                                                   >  1:887914/1‑140 (MQ=255)                                                                              gcttgcACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTgg                                                        <  2:293990/140‑1 (MQ=255)                                                                                ttgcACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGgt                                                      >  1:686847/1‑140 (MQ=255)                                                                                ttgcACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGgt                                                      >  1:49654/1‑140 (MQ=255)                                                                                ttgcACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGgt                                                      >  1:620016/1‑140 (MQ=255)                                                                                ttgcACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCAGTGGTTATGGCGGCGGTATTCAACTTTTTgggg                                                      >  1:820177/1‑139 (MQ=255)                                                                                 tgcACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGgtg                                                     >  1:351544/1‑140 (MQ=255)                                                                                 tgcACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGgtg                                                     >  1:555641/1‑140 (MQ=255)                                                                                       ggCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGc                                               <  1:515380/140‑1 (MQ=255)                                                                                       ggCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGc                                               <  1:797185/140‑1 (MQ=255)                                                                                        gCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCt                                              >  2:881269/1‑140 (MQ=255)                                                                                           tttGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTggg                                           <  1:170706/140‑1 (MQ=255)                                                                                                gCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTc                                      <  2:686847/140‑1 (MQ=255)                                                                                                gCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTc                                      <  2:620016/140‑1 (MQ=255)                                                                                                  tGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTg                                    >  2:523936/1‑140 (MQ=255)                                                                                                    ttCTACGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCTACTTTTTGGGTGTTTTTCTGGGGGGTCTGAg                                  >  2:210848/1‑140 (MQ=255)                                                                                                        aCGAAGCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGtt                              <  2:283746/140‑1 (MQ=255)                                                                                                             gCCATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTa                         <  2:49654/140‑1 (MQ=255)                                                                                                              ccATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGc                            <  1:746721/136‑1 (MQ=255)                                                                                                              ccATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGc                            >  2:746721/1‑136 (MQ=255)                                                                                                              ccATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTa                         >  1:960396/1‑139 (MQ=255)                                                                                                              ccATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTa                         <  2:960396/139‑1 (MQ=255)                                                                                                              ccATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTAt                        >  1:734485/1‑140 (MQ=255)                                                                                                              ccATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTAt                        >  1:406879/1‑140 (MQ=255)                                                                                                               cATCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTATg                       <  2:708565/140‑1 (MQ=255)                                                                                                                aTCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTAAGTGTTGCCTa                         >  2:833313/1‑137 (MQ=255)                                                                                                                aTCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCACGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTa                         <  1:833313/137‑1 (MQ=255)                                                                                                                 tCAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTATGcc                     >  1:716419/1‑140 (MQ=255)                                                                                                                  cAATGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTATGCCa                    <  2:406879/140‑1 (MQ=255)                                                                                                                   aaTGGTTTCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTATGCCAt                   <  2:887914/140‑1 (MQ=255)                                                                                                                        tttCCATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGGGGTCTGAGTGTTGCCTATGCCATTGTGc              >  1:509876/1‑140 (MQ=255)                                                                                                                           ccATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTATGCCATTGTGCata           >  2:95637/1‑140 (MQ=255)                                                                                                                            cATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTATg                       >  2:90644/1‑127 (MQ=255)                                                                                                                            cATGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTATg                       <  1:90644/127‑1 (MQ=255)                                                                                                                             aTGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTATGCCATTGTGCATAtg         <  2:555641/140‑1 (MQ=255)                                                                                                                             aTGACACAGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTATGCCATTGTGCATAtg         <  2:351544/140‑1 (MQ=255)                                                                                                                                 cacaGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTATGCCATTGTGCATAtgctgc     >  2:541450/1‑140 (MQ=255)                                                                                                                                    aGCCACCGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTATGCCATTGTGCATATGCTGCCGa  >  2:837389/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                       CATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATACCCGCGCGATGCGTTCTCAGCTCGCCGTGGTTATGGCGGCGGTATTCAACTTTTTGGGTGTTTTGCTGGGTGGTCTGAGTGTTGCCTATGCCATTGTGCATATGCTGCCGA  >  NC_000913/3637615‑3637885

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