Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,069,043	G→A	S472S (AGC→AGT)	pitB ←	phosphate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,069,043	0	G	A	100.0%	114.6 / NA	32	S472S (AGC→AGT)	pitB	phosphate transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/16);  total (16/16)

GCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGAC  >  REL606/3068970‑3069114                                                                          |                                                                        gCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATACt                                                                       <  1:774895/76‑1 (MQ=255)      gCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATACTGGTTa                                                                  <  2:729883/76‑1 (MQ=255)            cAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATACTGGTTAACGttt                                                            >  1:686176/1‑76 (MQ=255)              gTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATACTGGTTACCGTTTTa                                                          >  1:933454/1‑76 (MQ=255)                     ccaccaGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTa                                                   >  1:75336/1‑76 (MQ=255)                       accaGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAAc                                                 >  2:748256/1‑76 (MQ=255)                           gaaaaaaaAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCAc                                             <  1:363917/76‑1 (MQ=255)                            aaagaaaAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCaac                                            >  1:1142652/1‑74 (MQ=255)                            aaagaaaAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCAcc                                            >  2:889688/1‑76 (MQ=255)                            aaagaaaAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCAcc                                            >  1:1142664/1‑76 (MQ=255)                             aaaaaaaaTTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCAccg                                           <  2:310307/73‑1 (MQ=255)                                      tGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTccacc                                  <  1:295971/76‑1 (MQ=255)                                       gccgccGGTAAAGTAAATACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTccacca                                 <  1:397594/76‑1 (MQ=255)                                       gccgccGGTAAAGTAAATACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTccacca                                 <  1:1238990/76‑1 (MQ=255)                                       cccgcctGTAAAGTAAATACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTccacca                                 <  2:413504/75‑1 (MQ=255)                                         cgccgGTAAAGTAAATACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTATAACCCACCGCCGTccaccacc                               >  2:886147/1‑76 (MQ=255)                                         cgccgGTAAAGTAAAAACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTATAACCCACCGCCGTcccccacc                               >  2:886161/1‑76 (MQ=255)                                            cgGTAAAGTAAATACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATc                            <  2:1227270/76‑1 (MQ=255)                                               taaagtaaaTACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTc                         >  1:204070/1‑76 (MQ=255)                                                 aagtaaaTACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTccc                       >  2:505579/1‑76 (MQ=255)                                                    taaaTACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGc                    >  1:281890/1‑76 (MQ=255)                                                     aaaTACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCa                   >  1:702668/1‑76 (MQ=255)                                                     aaaTACCCACGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCa                   <  2:519491/76‑1 (MQ=255)                                                              cGCCATCAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGAc          <  2:426880/76‑1 (MQ=255)                                                               gCCATCAGGATACTGGTTACCGTTTTACGCTATAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACg         <  1:886659/76‑1 (MQ=255)                                                               gCCATCAGGATACTGGTTACCGTTTTACGCTATAACCCACCGCCGTCCACCAACATTGTCCCTGCAACTGCAGACg         <  1:886673/76‑1 (MQ=255)                                                                ccATCAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACga        >  2:132487/1‑76 (MQ=255)                                                                ccATCAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACga        >  1:458777/1‑76 (MQ=255)                                                                 cATCAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACgag       <  2:310790/76‑1 (MQ=255)                                                                 cATCAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACgag       <  1:74136/76‑1 (MQ=255)                                                                    cAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAgg    <  1:1092982/76‑1 (MQ=255)                                                                    cAGGATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAgg    <  1:1142234/76‑1 (MQ=255)                                                                      ggATACTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGAc  >  2:499005/1‑76 (MQ=255)                                                                          |                                                                        GCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGAC  >  REL606/3068970‑3069114

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