Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,507,942	A→G	Y217H (TAC→CAC)	php ←	phosphotriesterase homology protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,507,942	0	A	G	100.0%	226.1 / NA	66	Y217H (TAC→CAC)	php	phosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/34);  total (32/34)

TAAATAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGC  >  NC_000913/3507808‑3508071                                                                                                                                       |                                                                                                                                  tAAATAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTaa                                                                                                                               <  1:1113385/139‑1 (MQ=255) tAAATAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTaa                                                                                                                               <  1:146160/139‑1 (MQ=255)          aTAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTg                                                                                                                      >  2:633150/1‑139 (MQ=255)                aTAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCTCATGAAGCATCGCAATACGCTTTTCATCCTGGTGGTAACTGTTCTTGCCGATg                                                                                                                <  1:1202450/139‑1 (MQ=255)                aTAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTaa                                                                                                                               >  2:680735/1‑124 (MQ=255)                aTAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTaa                                                                                                                               <  1:680735/124‑1 (MQ=255)                 tAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATgg                                                                                                               >  2:1251875/1‑139 (MQ=255)                 tAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATgg                                                                                                               >  1:63986/1‑139 (MQ=255)                  aaCCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGgt                                                                                                              >  2:111167/1‑139 (MQ=255)                        cGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                                        >  2:969382/1‑139 (MQ=255)                           tGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                                        <  2:244437/136‑1 (MQ=255)                           tGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGaa                                                                                                        >  1:244437/1‑136 (MQ=255)                            ggCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGtt                                                                                                                          <  2:1230417/117‑1 (MQ=255)                            ggCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGtt                                                                                                                          >  1:1230417/1‑117 (MQ=255)                                 ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                               <  1:1251875/139‑1 (MQ=255)                                 ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                               <  2:1104112/139‑1 (MQ=255)                                 ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                               <  2:61557/139‑1 (MQ=255)                                 ttAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                               <  1:588531/139‑1 (MQ=255)                                  tAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCTGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                               <  2:514092/138‑1 (MQ=255)                                  tAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                               >  1:514092/1‑138 (MQ=255)                                          aGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCgag                                                                                      >  1:1203892/1‑139 (MQ=255)                                                gcgTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgatcgat                                                                                <  1:633150/139‑1 (MQ=255)                                                gcgTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgatcgat                                                                                <  1:969382/139‑1 (MQ=255)                                                 cgTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgatcgatc                                                                               <  2:1203892/139‑1 (MQ=255)                                                  gTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTGGCCGATGGTGTCGAACTGCACGTATGCGCCGAGATCGatca                                                                              <  1:15405/139‑1 (MQ=255)                                                         ccATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTGTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCaa                                                                       >  1:942395/1‑139 (MQ=255)                                                             cGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAAtg                                                                   >  1:295492/1‑139 (MQ=255)                                                              gACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGGGCCGAGATCGATCATCTTCAAAAtgt                                                                  >  2:1465256/1‑139 (MQ=255)                                                                    aTGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGc                                                                                                                     <  2:1539029/82‑1 (MQ=255)                                                                    aTGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGc                                                                                                                     >  1:1539029/1‑82 (MQ=255)                                                                      gACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAgg                                                          >  1:53533/1‑139 (MQ=255)                                                                             ttCAGCAACCCATGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                               >  2:1454706/1‑95 (MQ=255)                                                                             ttCAGCAACCCATGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                               <  1:1454706/95‑1 (MQ=255)                                                                             ttCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                               <  1:239594/95‑1 (MQ=255)                                                                             ttCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                               >  2:239594/1‑95 (MQ=255)                                                                                       cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAgg                                                          >  1:81742/1‑122 (MQ=255)                                                                                       cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAgg                                                          >  1:142881/1‑122 (MQ=255)                                                                                       cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAgg                                                          <  2:142881/122‑1 (MQ=255)                                                                                       cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAgg                                                          <  2:81742/122‑1 (MQ=255)                                                                                       cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGtt                                                        <  2:148773/124‑1 (MQ=255)                                                                                       cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGtt                                                        >  1:148773/1‑124 (MQ=255)                                                                                       cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCa                                         >  1:1292711/1‑139 (MQ=255)                                                                                       cACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCa                                         <  2:942395/139‑1 (MQ=255)                                                                                         cGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGt                                       <  2:63986/139‑1 (MQ=255)                                                                                               cgTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                                          >  2:1494532/1‑66 (MQ=255)                                                                                               cgTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCg                                                                                                          <  1:1494532/66‑1 (MQ=255)                                                                                               cgTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCaa                                 >  1:615204/1‑139 (MQ=255)                                                                                               cgTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGACGATGGTGTCGAACTGCACGTACGCGCCGAGAGCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCaa                                 >  2:1027448/1‑139 (MQ=255)                                                                                                      cATGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGAc                          >  2:278625/1‑139 (MQ=255)                                                                                                        tGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                                  <  1:718152/49‑1 (MQ=255)                                                                                                        tGAAGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGa                                                                                                                  >  2:718152/1‑49 (MQ=255)                                                                                                          aaGCATCGCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTTTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGAcgcgc                      <  1:86854/139‑1 (MQ=255)                                                                                                                 gCAATACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGAt               <  2:295492/139‑1 (MQ=255)                                                                                                                     tACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCg                                                             <  2:500888/89‑1 (MQ=255)                                                                                                                     tACGCTTTTCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCg                                                             >  1:500888/1‑89 (MQ=255)                                                                                                                           tttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgat                                                                                    >  2:1089075/1‑60 (MQ=255)                                                                                                                           tttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAgat                                                                                    <  1:1089075/60‑1 (MQ=255)                                                                                                                           tttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCaaaa                                                                     <  2:1418316/75‑1 (MQ=255)                                                                                                                           tttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCaaaa                                                                     >  1:1418316/1‑75 (MQ=255)                                                                                                                            ttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                               >  2:53107/1‑48 (MQ=255)                                                                                                                            ttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTa                                                                                               <  1:53107/48‑1 (MQ=255)                                                                                                                            ttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGtt                                                        >  2:698549/1‑87 (MQ=255)                                                                                                                            ttCGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGtt                                                        <  1:698549/87‑1 (MQ=255)                                                                                                                              cGTCCGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGc  <  2:574334/139‑1 (MQ=255)                                                                                                                                 ccGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACgcgc                                                                                          >  1:174987/1‑48 (MQ=255)                                                                                                                                 ccGGGTGGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACgcgc                                                                                          <  2:174987/48‑1 (MQ=255)                                                                                                                                       |                                                                                                                                  TAAATAGTCATAGCCATAACCACCGTTGGCTTTTAAATGGGAGCGGCGCGTAATATCCATCGACAGCATGACGCGGTTCAGCAACCCACGGTCGCGTAGCGCATGAAGCATCGCAATACGCTTTTCGTCCGGGTAGTAACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGC  >  NC_000913/3507808‑3508071

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