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_0009133,508,4510C.96.4% 226.2 / 1.8 56coding (140/879 nt)phpphosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base . (25/29);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

ATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAA  >  NC_000913/3508315‑3508583
                                                                                                                                        |                                                                                                                                    
aTGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGccc                                                                                                                                    <  1:1072946/139‑1 (MQ=255)
 tGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcccg                                                                                                                                   >  2:206971/1‑138 (MQ=255)
 tGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcccg                                                                                                                                   >  1:743266/1‑138 (MQ=255)
  gTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcccgg                                                                                                                                  <  1:921917/139‑3 (MQ=255)
    tCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCgggt                                                                                                                                >  2:357754/1‑138 (MQ=255)
             aaCGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTc                                                                                                                      >  1:505957/1‑139 (MQ=255)
                     cTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCAtct                                                                                                              <  2:1043265/139‑1 (MQ=255)
                       ggTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTcc                                                                                                            <  1:357754/139‑1 (MQ=255)
                              aCCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACCTCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     <  2:743266/139‑1 (MQ=255)
                              aCCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     >  1:439291/1‑139 (MQ=255)
                                  gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGccc                                                                                                                                    <  2:526909/105‑1 (MQ=255)
                                  gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGccc                                                                                                                                    >  1:526909/1‑105 (MQ=255)
                                  gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     <  2:386457/135‑1 (MQ=255)
                                  gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                                 <  1:781036/139‑1 (MQ=255)
                                  gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                                 >  2:781036/1‑139 (MQ=255)
                                  gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                                 <  2:439291/139‑1 (MQ=255)
                                  gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCAATGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     >  1:386457/1‑135 (MQ=255)
                                   tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  2:411900/139‑1 (MQ=255)
                                        gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  2:77292/134‑1 (MQ=255)
                                        gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                >  1:77292/1‑134 (MQ=255)
                                             cacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGAt                                                                                      >  2:614361/1‑139 (MQ=255)
                                                  tGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAgg                                                                                 >  1:608356/1‑139 (MQ=255)
                                                  tGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAgg                                                                                 >  1:819894/1‑139 (MQ=255)
                                                      ccccGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAggcggc                                                                             <  1:614361/139‑1 (MQ=255)
                                                      ccccGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAggcggc                                                                             <  1:341330/139‑1 (MQ=255)
                                                         cGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                >  1:2134/1‑117 (MQ=255)
                                                         cGTCTCGCGCATTACATCAAGCATAAATTGCCCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  2:2134/117‑1 (MQ=255)
                                                          gTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACAATACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAAt                                                                                                    <  1:136692/112‑1 (MQ=255)
                                                          gTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACAATACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAAt                                                                                                    >  2:136692/1‑112 (MQ=255)
                                                                cgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCAt                                                                                                                >  1:81563/1‑94 (MQ=255)
                                                                cgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCAt                                                                                                                <  2:81563/94‑1 (MQ=255)
                                                                            aGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGcc                                                       >  1:583308/1‑139 (MQ=255)
                                                                               aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTc                                                                                                                  >  2:238568/1‑77 (MQ=255)
                                                                               aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTc                                                                                                                  <  1:238568/77‑1 (MQ=255)
                                                                               aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     >  2:514326/1‑90 (MQ=255)
                                                                               aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     <  1:514326/90‑1 (MQ=255)
                                                                                 aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGt                                                                                                                    <  1:633072/73‑1 (MQ=255)
                                                                                 aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGt                                                                                                                    >  2:633072/1‑73 (MQ=255)
                                                                                 aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  <  1:761262/139‑1 (MQ=255)
                                                                                 aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  >  2:217302/1‑139 (MQ=255)
                                                                                  aaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                >  2:1021212/1‑92 (MQ=255)
                                                                                  aaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  1:1021212/92‑1 (MQ=255)
                                                                                   attgcgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgagaga                                                <  2:819894/139‑1 (MQ=255)
                                                                                   attgcgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgagaga                                                <  2:608356/139‑1 (MQ=255)
                                                                                       cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCaa                                            >  2:75038/1‑139 (MQ=255)
                                                                                            tgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAg                                                         >  1:591819/1‑121 (MQ=255)
                                                                                            tgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAg                                                         <  2:591819/121‑1 (MQ=255)
                                                                                                        aaCGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATgg                           >  2:744220/1‑139 (MQ=255)
                                                                                                           ggTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                        <  1:206971/139‑1 (MQ=255)
                                                                                                           ggTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                        <  2:583308/139‑1 (MQ=255)
                                                                                                           ggTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                        <  2:505957/139‑1 (MQ=255)
                                                                                                               ggTCATCTCAATCACATTACGCACGCCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCAGACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAgg                     >  1:679460/1‑139 (MQ=255)
                                                                                                               ggTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCaa                                            >  1:690022/1‑115 (MQ=255)
                                                                                                               ggTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCaa                                            <  2:690022/115‑1 (MQ=255)
                                                                                                               ggTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAggg                    >  1:851522/1‑139 (MQ=255)
                                                                                                                               ttACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     >  1:402192/1‑42 (MQ=255)
                                                                                                                               ttACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     <  2:402192/42‑1 (MQ=255)
                                                                                                                               ttACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACg                                                                     >  1:260839/1‑74 (MQ=255)
                                                                                                                               ttACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACg                                                                     <  2:260839/74‑1 (MQ=255)
                                                                                                                               ttACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgagaga                                                >  1:66415/1‑95 (MQ=255)
                                                                                                                               ttACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgagaga                                                <  2:66415/95‑1 (MQ=255)
                                                                                                                               ttACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATc    <  2:851522/139‑1 (MQ=255)
                                                                                                                                 acgcacgcCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaa  <  1:217302/139‑1 (MQ=255)
                                                                                                                                 acgcacgcCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaa  <  2:118734/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
ATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAA  >  NC_000913/3508315‑3508583

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: