Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,911,170	A→G	V119A (GTG→GCG)	pstS ←	phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,911,170	0	A	G	100.0%	120.3 / NA	34	V119A (GTG→GCG)	pstS	phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/9);  total (25/9)

ATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAAC  >  NC_000913/3911040‑3911303                                                                                                                                   |                                                                                                                                      aTCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGc                                                                                                                              >  2:864996/1‑140 (MQ=255)      cgcgGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGAc                                                                                                                         <  1:317338/140‑1 (MQ=255)        cggcgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACtt                                                                                                                       <  1:676038/140‑1 (MQ=255)           cgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAg                                                                                                                    >  1:11627/1‑140 (MQ=255)             tactacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGcc                                                                                                                  <  2:192444/140‑1 (MQ=255)             tactacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGcc                                                                                                                  <  1:1044167/140‑1 (MQ=255)               ctacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCt                                                                                                                <  2:550830/140‑1 (MQ=255)                     cAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATg                                                                                                          >  2:665103/1‑140 (MQ=255)                     cAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATg                                                                                                          >  1:597954/1‑140 (MQ=255)                        tGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGc                                                                                                                   >  2:714890/1‑128 (MQ=255)                        tGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGc                                                                                                                   <  1:714890/128‑1 (MQ=255)                          ttttGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAAc                                                                                                     >  2:80636/1‑140 (MQ=255)                          ttttGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAAc                                                                                                     >  1:46126/1‑140 (MQ=255)                                     gatttccAGCCCGGGTTTCGGTTTGGCGATGGCTTCATCACCCCCCTTCTTGATTTTCCCCGGGTGGGTGTCGCGGGGGGTTTTCCAATCAGGGGCCATTTCGCTAGCCTCCCGCCCGGAATTGTTAACCGCCAGcacca                                                                                          <  2:302783/134‑1 (MQ=255)                                         ttCAGACCCGGATTCATTTTGGCGATGGTTTCATCACCCCACTTCTTGATTTTGCCCGGGTAGATGCCGCCGGGGTTTTTCCCATCCAGCCCCAGTTCGCCAGCCTTCAGCCCTGAAATGTTAACCGCCAGCACCAcgcc                                                                                      <  2:853516/140‑1 (MQ=255)                                                  ggATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCAc                                                                             >  2:707875/1‑140 (MQ=255)                                                   gATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACg                                                                            >  2:599115/1‑140 (MQ=255)                                                       cattttGGCGATGGCTTCATCACCCCTCTTCTTGATTTTGCCCGGGTAGATGTCGCCGAGGTTTTTCCCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGATTGTTAACCGCCAGCACCACGCCGCCAATCACGGTCg                                                                        <  2:46126/137‑1 (MQ=255)                                                       cAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGCCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCg                                                                        <  1:557342/140‑1 (MQ=255)                                                          tttGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGa                                                                     >  1:142377/1‑140 (MQ=255)                                                            gggCGATGGCTTCATCACCCCACTTCTTGATTTTGCCCGGGTAGATGTCGCCGGGGTTTTTCCCACCCAGCGCCAGTTCGCCAGACTTCAGCCCGGAAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAAc                                                                   <  2:785664/139‑1 (MQ=255)                                                                  tGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAAc                                                             >  2:257485/1‑140 (MQ=255)                                                                    gCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCAc                                                                             <  1:513035/122‑1 (MQ=255)                                                                    gCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCAc                                                                             >  2:513035/1‑122 (MQ=255)                                                                                 ccttcttGATTTTGCCCGGATAGTGGTCGCGGGGGGTTTTACCACCCACCGCCAGTTCGCCATACTTCAGCCCTGGATTGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGCCCTTCCTGAGcc                                              <  2:142377/139‑1 (MQ=255)                                                                                  cttcttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCa                                             >  1:172109/1‑140 (MQ=255)                                                                                   ttcttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAg                                            >  2:495607/1‑140 (MQ=255)                                                                                             tGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCa                                  >  1:983470/1‑140 (MQ=255)                                                                                               cccAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGa                                >  2:16992/1‑140 (MQ=255)                                                                                                   ggTAGATGTCGCCGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAgc                            >  2:871066/1‑140 (MQ=255)                                                                                                    gTAGATGTCGCCGGGGTTTTTCCCATCCAGCGCCAGTCCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAgcg                           <  2:597954/140‑1 (MQ=255)                                                                                                    gTAGATGTCGCCGGGGTTTTTCCCATCCAGCGCAAGTCCGCAAGCCTTCGGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAgcg                           <  2:915098/140‑1 (MQ=255)                                                                                                              ccGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGNCTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCaga                 >  1:254190/1‑140 (MQ=255)                                                                                                              ccGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCaga                 >  2:154643/1‑140 (MQ=255)                                                                                                              ccGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCaga                 >  2:423175/1‑140 (MQ=255)                                                                                                              ccGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCaga                 >  1:915098/1‑140 (MQ=255)                                                                                                              ccGAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCaga                 >  2:77533/1‑140 (MQ=255)                                                                                                                gAGGGTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAgg               <  1:80636/140‑1 (MQ=255)                                                                                                                  gggTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCa             >  2:1056020/1‑140 (MQ=255)                                                                                                                  gggTTTTACCATCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCa             >  2:605577/1‑140 (MQ=255)                                                                                                                            aTCCAGCGCCATTCCGCCAGACTTCACCCCTGAAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCaa   <  2:172109/140‑1 (MQ=255)                                                                                                                            aTCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCaa   <  1:495607/140‑1 (MQ=255)                                                                                                                            aTCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCaa   <  2:983470/140‑1 (MQ=255)                                                                                                                             tCCAGCGCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAAc  >  2:390897/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                      ATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAAC  >  NC_000913/3911040‑3911303

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