Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,508,454	Δ1 bp	coding (137/879 nt)	php ←	phosphotriesterase homology protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,451	0	C	.	100.0%	290.0 / NA	69	coding (140/879 nt)	php	phosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (31/38);  total (31/38)

ATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATC  >  NC_000913/3508315‑3508581                                                                                                                                         |                                                                                                                                   aTGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGccc                                                                                                                                  <  1:472022/139‑1 (MQ=255)                        ggTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCAt                                                                                                              <  2:462200/135‑1 (MQ=255)                        ggTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCAt                                                                                                              >  1:462200/1‑135 (MQ=255)                               aCCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                   <  2:1236195/139‑1 (MQ=255)                               aCCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                   >  2:1167715/1‑139 (MQ=255)                                 cGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                 <  1:900709/139‑1 (MQ=255)                                    tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                              <  2:1307245/139‑1 (MQ=255)                                    tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                              <  2:677934/139‑1 (MQ=255)                                         gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                         >  1:1238352/1‑139 (MQ=255)                                         gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                         <  1:1167715/139‑1 (MQ=255)                                         gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                         >  1:337293/1‑139 (MQ=255)                                                       ccccGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                 <  1:1172210/117‑1 (MQ=255)                                                       ccccGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                 >  2:1172210/1‑117 (MQ=255)                                                       ccccGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                         >  1:1381489/1‑125 (MQ=255)                                                       ccccGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                         <  2:1381489/125‑1 (MQ=255)                                                          cGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGt                                                                        <  2:1238352/139‑1 (MQ=255)                                                            tctcGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                 >  1:503242/1‑112 (MQ=255)                                                            tctcGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                 <  2:503242/112‑1 (MQ=255)                                                            tctcGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACg                                                                                             <  1:433678/116‑1 (MQ=255)                                                            tctcGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACg                                                                                             >  2:433678/1‑116 (MQ=255)                                                                 cgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAg                                                                                                                       >  1:144111/1‑85 (MQ=255)                                                                 cgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAg                                                                                                                       <  2:144111/85‑1 (MQ=255)                                                                 cgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                               >  1:521211/1‑109 (MQ=255)                                                                 cgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                               <  2:521211/109‑1 (MQ=255)                                                                              gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                   <  2:1424033/92‑1 (MQ=255)                                                                              gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                   >  1:1424033/1‑92 (MQ=255)                                                                              gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCg                                                    >  2:635953/1‑139 (MQ=255)                                                                                aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGga                                                  >  2:1202800/1‑139 (MQ=255)                                                                                    attgcgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgagaga                                              <  2:337293/139‑1 (MQ=255)                                                                                         gcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAat                                         <  1:635953/139‑1 (MQ=255)                                                                                             tgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                   <  2:529892/77‑1 (MQ=255)                                                                                             tgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                   >  1:529892/1‑77 (MQ=255)                                                                                             tgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGc                                     >  1:1370161/1‑139 (MQ=255)                                                                                               cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAg                                   <  2:1370161/139‑1 (MQ=255)                                                                                                gcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGa                                  <  1:1420814/139‑1 (MQ=255)                                                                                                     aTGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                              >  1:695239/1‑74 (MQ=255)                                                                                                     aTGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                              <  2:695239/74‑1 (MQ=255)                                                                                                     aTGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACg                                                                   <  2:359792/101‑1 (MQ=255)                                                                                                     aTGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACg                                                                   >  1:359792/1‑101 (MQ=255)                                                                                                     aTGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTc                             <  1:1202800/139‑1 (MQ=255)                                                                                                      tGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCa                            <  1:100249/139‑1 (MQ=255)                                                                                                            ggTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                      <  2:417476/139‑1 (MQ=255)                                                                                                            ggTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                      <  2:816709/139‑1 (MQ=255)                                                                                                              ttGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                 >  1:135691/1‑62 (MQ=255)                                                                                                              ttGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                 <  2:135691/62‑1 (MQ=255)                                                                                                                ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCAtct                                                                                                            >  1:606949/1‑49 (MQ=39)                                                                                                                ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCAtct                                                                                                            <  2:606949/49‑1 (MQ=39)                                                                                                                ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                 >  1:1002539/1‑60 (MQ=255)                                                                                                                ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                 <  2:1002539/60‑1 (MQ=255)                                                                                                                ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAgg                                                                               >  1:296632/1‑78 (MQ=255)                                                                                                                ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAgg                                                                               <  2:296632/78‑1 (MQ=255)                                                                                                                ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                             <  2:61594/112‑1 (MQ=255)                                                                                                                ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                             >  1:61594/1‑112 (MQ=255)                                                                                                                ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAggg                  >  2:771496/1‑139 (MQ=255)                                                                                                                 gTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                   >  1:636795/1‑57 (MQ=255)                                                                                                                 gTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                   <  2:636795/57‑1 (MQ=255)                                                                                                                 gTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAg                                                       <  1:1042715/101‑1 (MQ=255)                                                                                                                 gTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAg                                                       >  2:1042715/1‑101 (MQ=255)                                                                                                                 gTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                      >  2:1126039/1‑134 (MQ=255)                                                                                                                 gTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                      <  1:1126039/134‑1 (MQ=255)                                                                                                                     tctcAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTaa             <  1:771496/139‑1 (MQ=255)                                                                                                                        cAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTgg                                                                                                       <  2:502675/46‑1 (MQ=39)                                                                                                                        cAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTgg                                                                                                       >  1:502675/1‑46 (MQ=39)                                                                                                                                ttACGCACG‑CCCGTGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGga                                                  <  2:1112339/91‑1 (MQ=255)                                                                                                                                ttACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGga                                                  >  1:1112339/1‑91 (MQ=255)                                                                                                                                ttACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATc  >  1:63789/1‑139 (MQ=255)                                                                                                                                ttACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGCTTTTAAAGCCGgaga                                                <  2:955928/93‑1 (MQ=255)                                                                                                                                ttACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGACAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                >  1:955928/1‑93 (MQ=255)                                                                                                                                 tACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                   >  2:531303/1‑41 (MQ=255)                                                                                                                                 tACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                   <  1:531303/41‑1 (MQ=255)                                                                                                                                   cgcacg‑cCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                              <  1:198843/44‑1 (MQ=255)                                                                                                                                   cgcacg‑cCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                              >  2:198843/1‑44 (MQ=255)                                                                                                                                         |                                                                                                                                   ATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATC  >  NC_000913/3508315‑3508581

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