Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,508,245	T→C	T116A (ACG→GCG)	php ←	phosphotriesterase homology protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,245	0	T	C	100.0%	172.3 / NA	49	T116A (ACG→GCG)	php	phosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/31);  total (18/31)

GCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCA  >  NC_000913/3508119‑3508390                                                                                                                               |                                                                                                                                                  gCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAAttt                                                                                                                            <  2:777187/150‑1 (MQ=255)   gTGGAGATCAGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCAt                                                                                                                         >  1:431581/1‑151 (MQ=255)     ggAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCg                                                                                                                       <  2:129699/151‑1 (MQ=255)       agaTCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGAc                                                                                                                     <  1:304629/151‑1 (MQ=255)         aTCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGAcc                                                                                                                    >  1:211563/1‑150 (MQ=255)         aTCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCa                                                                                                                   <  2:593839/151‑1 (MQ=255)            gggCGTCCGGTCTGGTTATGCGCCAGCGCCGCGGCAATAAATCCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGTTCTCCGCTATGATCCCGGCTTTCAGCTACGCGCCATCGATACCCTGTTCAATTTCATCGACCAtct                                                                                                                <  1:647848/151‑1 (MQ=255)            gggCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCAtct                                                                                                                >  1:629741/1‑151 (MQ=255)                  ccGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTgc                                                                                                           <  2:567512/150‑1 (MQ=255)                   cGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTgcgc                                                                                                         <  2:394639/151‑1 (MQ=255)                     gTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGcc                                                                                                        >  2:545070/1‑150 (MQ=255)                     gTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCa                                                                                                       <  1:682417/151‑1 (MQ=255)                               acgcCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCAc                                                                                             <  2:239966/150‑1 (MQ=255)                                  ccAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCt                                                                                          <  2:306374/151‑1 (MQ=255)                                   cagcgcagcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTg                                                                                         <  2:150910/151‑1 (MQ=255)                                    agcgcagcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGc                                                                                        >  1:487265/1‑151 (MQ=255)                                        cagcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGt                                                                                    <  1:695073/151‑1 (MQ=255)                                        cagcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGt                                                                                    <  1:770973/151‑1 (MQ=255)                                         agcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTc                                                                                   <  1:467564/151‑1 (MQ=255)                                         agcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTc                                                                                   <  1:545033/151‑1 (MQ=255)                                            ggCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCAAGTTCCTGCACGCTGCGGGTCg                                                                                  <  1:729259/149‑1 (MQ=255)                                              cAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCcac                                                                              >  1:200581/1‑151 (MQ=255)                                                      aCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTcc                                                                       <  1:497752/150‑1 (MQ=255)                                                         ttCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCggg                                                                    >  2:726767/1‑150 (MQ=255)                                                             ctTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGaaa                                                                 <  2:506635/149‑1 (MQ=255)                                                                ccAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAcg                                                             <  2:115836/150‑1 (MQ=255)                                                                  aGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGt                                                          >  1:158513/1‑151 (MQ=255)                                                                  aGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGt                                                          >  1:310732/1‑151 (MQ=255)                                                                    cggcgTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTcc                                                        >  2:837597/1‑151 (MQ=255)                                                                     ggcgTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCt                                                       <  2:848476/151‑1 (MQ=255)                                                                     ggcgTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCt                                                       >  2:431359/1‑151 (MQ=255)                                                                       cgTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTgg                                                     >  2:739434/1‑151 (MQ=255)                                                                               ttCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAAc                                              <  2:653783/150‑1 (MQ=255)                                                                                tCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGtaataa                                               <  1:17227/148‑1 (MQ=255)                                                                                 ccTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCg                                            <  2:30022/150‑1 (MQ=255)                                                                                    tCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTAc                                        >  1:247482/1‑151 (MQ=255)                                                                                     cGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACa                                       <  1:242387/151‑1 (MQ=255)                                                                                            gCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGccac                                 <  1:557713/150‑1 (MQ=255)                                                                                             ccGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGccacc                                >  1:204699/1‑150 (MQ=255)                                                                                                   tcCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTg                          >  2:169894/1‑150 (MQ=255)                                                                                                     cgcgATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATc                       <  1:275911/151‑1 (MQ=255)                                                                                                       cgATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATccc                     <  1:789384/151‑1 (MQ=255)                                                                                                            atCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGtct                >  1:183312/1‑151 (MQ=255)                                                                                                            atCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGtct                <  1:629593/151‑1 (MQ=255)                                                                                                             tCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGtctc               >  2:517276/1‑151 (MQ=255)                                                                                                                      tCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTa       <  1:111487/150‑1 (MQ=255)                                                                                                                       cAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCATGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTAc      <  2:377427/150‑1 (MQ=255)                                                                                                                          cTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACa     <  1:851314/148‑1 (MQ=255)                                                                                                                            ccGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCa  >  2:473512/1‑149 (MQ=255)                                                                                                                               |                                                                                                                                                  GCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCA  >  NC_000913/3508119‑3508390

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