Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,910,707	A→T	D273E (GAT→GAA)	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,910,707	0	A	T	100.0%	219.9 / NA	62	D273E (GAT→GAA)	pstS	phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (38/24);  total (38/24)

ACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAA  >  NC_000913/3910610‑3910801                                                                                                  |                                                                                               aCGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTc                                                                                              <  1:1360908/100‑1 (MQ=255)  ngCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCt                                                                                             >  2:479499/2‑100 (MQ=255)  cGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCt                                                                                             >  2:554742/1‑100 (MQ=255)  cGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCt                                                                                             >  2:1171418/1‑100 (MQ=255)  cGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCTTc                                                                                              >  2:111745/1‑99 (MQ=255)    cccAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCttcttc                                                                                           >  2:980857/1‑100 (MQ=255)    cccAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCttcttc                                                                                           <  1:915549/100‑1 (MQ=255)     ccAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCg                                                                                          <  1:1533847/100‑1 (MQ=255)           gaagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTc                                                                                    >  2:393813/1‑100 (MQ=255)                aTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCtggtt                                                                               >  1:1317469/1‑100 (MQ=255)                aTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCtggtt                                                                               >  2:1511566/1‑100 (MQ=255)                 tttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCtggttg                                                                              <  1:980865/100‑1 (MQ=255)                  ttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCtggttgg                                                                             <  1:723249/100‑1 (MQ=255)                   tCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCtggttggt                                                                            <  1:115491/100‑1 (MQ=255)                    cAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTc                                                                           >  2:1539524/1‑100 (MQ=255)                     aGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCa                                                                          >  1:293127/1‑100 (MQ=255)                     aGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCa                                                                          <  1:709058/100‑1 (MQ=255)                          ttCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATcc                                                                     <  1:616613/100‑1 (MQ=255)                           tCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCt                                                                    >  1:1175917/1‑100 (MQ=255)                              gtgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATACGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAg                                                                 >  2:795574/1‑100 (MQ=255)                               tgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGc                                                                <  2:1317460/100‑1 (MQ=255)                                  cTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGaa                                                             >  2:55075/1‑100 (MQ=255)                                  cTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTACGCCTTTCAGGTTGGTCAGATCCTGAGCGaa                                                             >  2:1061775/1‑100 (MQ=255)                                   ttGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAg                                                            <  1:1192404/100‑1 (MQ=255)                                   ttGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAg                                                            <  2:715109/100‑1 (MQ=255)                                        cTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGtttt                                                       >  2:69145/1‑100 (MQ=255)                                          ggTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTTGTCAGATCCTGAGCGAAGGTTTTGc                                                     >  2:869121/1‑100 (MQ=255)                                          ggTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATACTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGc                                                     >  2:839574/1‑100 (MQ=255)                                           gTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCt                                                    >  2:158329/1‑100 (MQ=255)                                             ttcttcTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTcn                                                  <  1:607715/100‑2 (MQ=255)                                               cttctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAg                                                >  1:1201375/1‑100 (MQ=255)                                                ttctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGt                                               >  1:141420/1‑100 (MQ=255)                                                 tctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTc                                              >  1:1508310/1‑100 (MQ=255)                                                   tgttCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTg                                            >  2:205628/4‑100 (MQ=255)                                                    gttCTTTGTGGATCAGAATGAACGTGGTAGAGGTGATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCGAGTCTGc                                           <  2:825985/98‑1 (MQ=255)                                                    gttCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGc                                           <  2:1257642/98‑1 (MQ=255)                                                          tgtgGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCttt                                     >  2:1386634/1‑100 (MQ=255)                                                            tgGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTtg                                   >  1:1170024/1‑100 (MQ=255)                                                              naTCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTAAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTtgct                                 >  2:1495500/2‑100 (MQ=255)                                                              gATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTtgct                                 >  2:375890/1‑100 (MQ=255)                                                              gATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTtgct                                 >  1:1345109/1‑100 (MQ=255)                                                               aTCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTtgctg                                <  2:880918/100‑1 (MQ=255)                                                               aTCAGAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTtgctg                                <  2:1307197/100‑1 (MQ=255)                                                                   gAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCAtt                            >  1:1430465/1‑100 (MQ=255)                                                                   gAATGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCAtt                            <  1:280354/100‑1 (MQ=255)                                                                    atctaaCGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTa                           <  1:1205720/96‑1 (MQ=255)                                                                    aaTGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTa                           <  1:1290063/100‑1 (MQ=255)                                                                      tGAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGc                         >  2:1039986/1‑100 (MQ=255)                                                                       gAACGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCg                        <  2:1047414/100‑1 (MQ=255)                                                                        aaCGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTaagat                       >  2:1205713/1‑96 (MQ=255)                                                                         aCGTGGTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGAGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCAtt                            >  2:1453605/1‑94 (MQ=255)                                                                             ggTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGttt                  >  2:336648/1‑100 (MQ=255)                                                                              gTAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGtttt                 >  2:1081668/1‑100 (MQ=255)                                                                               tAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTttc                >  2:1476616/1‑100 (MQ=255)                                                                               tAGAGGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTttc                >  2:1526909/1‑100 (MQ=255)                                                                                 gagGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTttctt              <  2:1170017/100‑1 (MQ=255)                                                                                 gagGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTttctt              >  2:357239/1‑100 (MQ=255)                                                                                  agGTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTttcttc             <  1:55075/100‑1 (MQ=255)                                                                                   ggTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTtcg            >  2:1292587/1‑100 (MQ=255)                                                                                   ggTAATAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTtcg            >  2:1420528/1‑100 (MQ=255)                                                                                      aaTAGGCCATGCTTCATCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTtcggtc         <  1:1547839/100‑1 (MQ=255)                                                                                        tAGGCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTtcggtcgg       <  1:799723/100‑1 (MQ=255)                                                                                          ggCCATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGAc     <  1:718467/100‑1 (MQ=255)                                                                                             cATGCTTCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTaa  >  2:1485274/1‑100 (MQ=255)                                                                                                  |                                                                                               ACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAA  >  NC_000913/3910610‑3910801

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: ‑