Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,909,930	C→T	V157I (GTT→ATT)	pstC ←	phosphate ABC transporter membrane subunit PstC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,909,930	0	C	T	100.0%	212.7 / NA	58	V157I (GTT→ATT)	pstC	phosphate ABC transporter membrane subunit PstC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/33);  total (25/33)

AATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAG  >  NC_000913/3909853‑3910005                                                                              |                                                                            aataatCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCaataat                                                                          >  2:762550/1‑81 (MQ=255)  ataatCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATc                                                                         <  1:762484/81‑1 (MQ=255)   taatCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCg                                                                        <  2:762485/81‑1 (MQ=255)    aatCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCgg                                                                       >  2:762544/1‑81 (MQ=255)    aatCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCgg                                                                       >  1:762509/1‑81 (MQ=255)     atCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCggg                                                                      <  2:762486/81‑1 (MQ=255)     atCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCgg                                                                       <  1:762487/80‑1 (MQ=255)      tCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGa                                                                     >  2:762547/1‑81 (MQ=255)        atgatgGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGata                                                                   <  1:762488/81‑1 (MQ=255)         tgatgGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGatat                                                                  <  1:762489/81‑1 (MQ=255)           atgGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATAtt                                                                 <  1:762490/80‑1 (MQ=255)           atgGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATAtt                                                                 <  1:762492/80‑1 (MQ=255)           atgGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTc                                                                >  1:762560/1‑81 (MQ=255)           atgGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTc                                                                <  2:762491/81‑1 (MQ=255)              gCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGAc                                                             <  2:762493/81‑1 (MQ=255)               cGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACa                                                            <  1:762494/81‑1 (MQ=255)                gAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACAt                                                           >  2:762563/1‑81 (MQ=255)                  ggATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGa                                                         >  2:762558/1‑81 (MQ=255)                   gATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGat                                                        <  2:762495/81‑1 (MQ=255)                       aCGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATAtt                                                     <  1:762496/80‑1 (MQ=255)                       aCGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTg                                                    >  1:762557/1‑81 (MQ=255)                       aCGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTg                                                    >  1:762582/1‑81 (MQ=255)                        cGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGc                                                   <  2:762497/81‑1 (MQ=255)                         gCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGcc                                                  >  1:762553/1‑81 (MQ=255)                          ccTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCg                                                 <  2:762498/81‑1 (MQ=255)                          ccTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCg                                                 <  1:762500/81‑1 (MQ=255)                          ccTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCg                                                 <  1:762499/81‑1 (MQ=255)                           cTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGa                                                >  2:762559/1‑81 (MQ=255)                             gCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGAcc                                              >  1:762552/1‑81 (MQ=255)                              ccGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGAcc                                              <  1:762502/80‑1 (MQ=255)                              ccGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGAc                                               <  2:762501/79‑1 (MQ=255)                              ccGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCg                                             >  1:762571/1‑81 (MQ=255)                                  gAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTc                                         >  1:762586/1‑81 (MQ=255)                                  gAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTc                                         >  2:762580/1‑81 (MQ=255)                                      ataccgataccAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTg                                      >  1:762567/1‑80 (MQ=255)                                      ataccgataccAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGa                                     <  2:762503/81‑1 (MQ=255)                                           gataccAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGt                                 >  1:762592/1‑80 (MQ=255)                                           gataccAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTa                                <  1:762504/81‑1 (MQ=255)                                            ataccAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTa                                >  2:762561/1‑80 (MQ=255)                                               ccAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACg                            <  1:762505/81‑1 (MQ=255)                                                cAAATGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACgg                           >  2:762564/1‑81 (MQ=255)                                                 aaaTGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGc                          <  1:762506/81‑1 (MQ=255)                                                  aaTGCGGGGCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCg                         <  1:762507/81‑1 (MQ=255)                                                   aTGCGGGGCCAGAGAACAGCGCGCCAATAACCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCg                         >  1:762569/1‑80 (MQ=255)                                                       ggggCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAg                    >  2:762556/1‑81 (MQ=255)                                                        gggCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAgc                   <  2:762508/81‑1 (MQ=255)                                                         ggCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAgcg                  >  1:762575/1‑81 (MQ=255)                                                         ggCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAgcg                  <  2:762509/81‑1 (MQ=255)                                                         ggCCAGAGAACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAgcg                  <  2:762510/81‑1 (MQ=255)                                                              gagaACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTGAACGGCGAACAGCGGCGCa             <  2:762511/81‑1 (MQ=255)                                                                gaACAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCCCCTGAAAGTAAACGGCGAACAGCGGCGCa             >  2:762593/1‑79 (MQ=255)                                                                 aaCAGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAg          <  2:762513/81‑1 (MQ=255)                                                                 aaCAGCGCGCCAATAATCGGGATATTCGACAAGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAg          <  2:762512/81‑1 (MQ=255)                                                                    aGCGCGCCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGAt        >  1:762588/1‑80 (MQ=255)                                                                     gcgcgcCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATaa      <  1:762514/81‑1 (MQ=255)                                                                       gcgcCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAAc    <  1:762515/81‑1 (MQ=255)                                                                        cgcCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATaaa     >  2:762584/1‑79 (MQ=255)                                                                         gcCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAg  <  2:762516/81‑1 (MQ=255)                                                                         gcCAATAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAg  >  2:762577/1‑81 (MQ=255)                                                                              |                                                                            AATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAG  >  NC_000913/3909853‑3910005

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