Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,910,688	T→C	T280A (ACC→GCC)	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,688	0	T	C	100.0%	162.0 / NA	46	T280A (ACC→GCC)	pstS	phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/22);  total (24/22)

CTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGAT  >  NC_000913/3910555‑3910823                                                                                                                                      |                                                                                                                                        cTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGt                                                                                                                                   >  2:798985/1‑140 (MQ=255) cTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGt                                                                                                                                   >  1:469024/1‑140 (MQ=255) cTATCCGGCAGGCTAGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGt                                                                                                                                   >  2:356022/1‑140 (MQ=255)   aTCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTaa                                                                                                                                 <  2:469024/140‑1 (MQ=255)    tCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAAt                                                                                                                                <  2:386255/140‑1 (MQ=255)    tCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAAt                                                                                                                                >  2:673550/1‑140 (MQ=255)     ccGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATa                                                                                                                               <  1:412289/140‑1 (MQ=255)         cAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATa                                                                                                                               <  2:91236/136‑1 (MQ=255)         cAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATa                                                                                                                               >  1:91236/1‑136 (MQ=255)              tGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAACGAACGTGGCAGAGGTAATAGGCCATGCa                                                                                                                      >  2:775419/1‑140 (MQ=255)               ggCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCAt                                                                                                                     <  2:192848/140‑1 (MQ=255)               ggCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCAt                                                                                                                     <  2:41561/140‑1 (MQ=255)                gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAACGAACGTGGCAGAGGTAATAGGCCATGCATc                                                                                                                    <  1:775419/140‑1 (MQ=255)                 cGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCATTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCt                                                                                                                   <  1:280640/140‑1 (MQ=255)                          ggTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTc                                                                                                          >  2:58398/1‑140 (MQ=255)                           gTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCt                                                                                                         >  1:461357/1‑140 (MQ=255)                             cGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCtgg                                                                                                       <  2:461357/140‑1 (MQ=255)                              gTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggt                                                                                                      <  1:383747/140‑1 (MQ=255)                              gTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggt                                                                                                      <  1:62198/140‑1 (MQ=255)                              gTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggt                                                                                                      <  2:647102/140‑1 (MQ=255)                              attCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggt                                                                                                      <  2:708193/139‑1 (MQ=255)                               ttCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt                                                                                                     >  1:270573/1‑140 (MQ=255)                                tCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttg                                                                                                    >  2:617481/1‑140 (MQ=255)                                               cGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCg                                                                                     >  2:366923/1‑140 (MQ=255)                                                ggTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGa                                                                                    <  1:195939/140‑1 (MQ=255)                                                   tttGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAgg                                                                                 <  1:798985/140‑1 (MQ=255)                                                    ttGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGt                                                                                >  1:84772/1‑140 (MQ=255)                                                             cAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCa                                                                       <  2:591129/140‑1 (MQ=255)                                                                tCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTc                                                                    >  2:192538/1‑140 (MQ=255)                                                                tCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTc                                                                    >  1:133527/1‑140 (MQ=255)                                                                          tCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCtttt                                                          >  1:336493/1‑140 (MQ=255)                                                                                cTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCa                                                    >  1:285268/1‑140 (MQ=255)                                                                                        ccTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGaa                                            >  2:400142/1‑140 (MQ=255)                                                                                           tGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGtt                                         >  2:61838/1‑140 (MQ=255)                                                                                             ttCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGtttt                                       >  2:102564/1‑140 (MQ=255)                                                                                             ttCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGtttt                                       >  1:364217/1‑140 (MQ=255)                                                                                                       ttctGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTtcggtcgg                             >  1:640178/1‑140 (MQ=255)                                                                                                          tGATCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACt                          <  2:336493/140‑1 (MQ=255)                                                                                                            aTCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTaa                        >  1:429583/1‑140 (MQ=255)                                                                                                            aTCTTTGTGGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTaa                        <  1:356022/140‑1 (MQ=255)                                                                                                                 tgtgGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGt                   >  2:348801/1‑140 (MQ=255)                                                                                                                   tgGATCAGAATGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGttt                 <  2:178184/140‑1 (MQ=255)                                                                                                                           aaTGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAg         >  1:23035/1‑140 (MQ=255)                                                                                                                           aaTGAACGTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAg         <  2:270573/140‑1 (MQ=255)                                                                                                                                  gTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGAt  <  1:348801/140‑1 (MQ=255)                                                                                                                                  gTGGCAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGAt  <  1:617481/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                        CTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGAT  >  NC_000913/3910555‑3910823

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