Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,791,448 C→G 9.8% G283G (GGC→GGG)  PP_3353 → sulfatase‑modifying protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,791,4480CG 9.8% 149.1 / 7.8 51G283G (GGC→GGG) PP_3353sulfatase‑modifying protein
Reads supporting (aligned to +/- strand):  ref base C (26/20);  new base G (2/3);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

TTTTACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGAC  >  NC_002947/3791304‑3791591
                                                                                                                                                |                                                                                                                                               
ttttACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggcccggc                                                                                                                                            <  1:790377/150‑1 (MQ=255)
      cGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGGCCGGGCCTGAg                                                                                                                                      <  1:5127/150‑1 (MQ=255)
       ggTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGGCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGt                                                                                                                                     >  2:461441/1‑150 (MQ=255)
        gTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTg                                                                                                                                    >  1:551723/1‑150 (MQ=255)
            aTCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAg                                                                                                                                >  1:666913/1‑150 (MQ=255)
            aTCGCCCGTGGATTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCCCGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAg                                                                                                                                >  2:122356/1‑150 (MQ=255)
               gCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCtt                                                                                                                             >  2:682178/1‑150 (MQ=255)
                 ccGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTcc                                                                                                                           >  2:212499/1‑149 (MQ=255)
                 ccGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTcc                                                                                                                           >  2:212119/1‑149 (MQ=255)
                  cGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGa                                                                                                                          >  2:1066207/1‑150 (MQ=255)
                    tGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGAcc                                                                                                                        <  2:449807/150‑1 (MQ=255)
                       gTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCCCAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCgg                                                                                                                     >  1:852503/1‑150 (MQ=255)
                            ttCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCa                                                                                                                   <  1:757641/147‑1 (MQ=255)
                            ttCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCa                                                                                                                   >  2:757641/1‑147 (MQ=255)
                                    cAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGt                                                                                                        <  1:1066207/150‑1 (MQ=255)
                                                ccTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCgtgccggtgccgct                                                                                            <  1:18163/150‑3 (MQ=255)
                                                ccTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGt                                                                                            >  1:826899/1‑150 (MQ=255)
                                                    ttCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACATCTTCTACCCTTCAGACCCCGGCGTGCCGGTGGCGGTAATc                                                                                        >  2:1030248/1‑150 (MQ=255)
                                                               cggcAATGTCTGGGGGCTGTCTCGCACGGGCGAGCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCTGTGGCGGTAATCAAAGGCGGCTc                                                                             <  2:423841/150‑1 (MQ=255)
                                                                     tGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCAATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCt                                                                       <  2:67867/150‑1 (MQ=255)
                                                                       tCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGa                                                                                                                                       >  1:957504/1‑84 (MQ=255)
                                                                       tCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGa                                                                                                                                       <  2:957504/84‑1 (MQ=255)
                                                                            gAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGCAATCAAAGGCGGCTCGCACCTGTGTTcg                                                                >  1:734690/1‑150 (MQ=255)
                                                                              gCTGACTCGCACGGGCTAGCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCACAGGCGGCTCGCACCTGTGTTcggc                                                              <  1:386720/150‑1 (MQ=255)
                                                                                  aCTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGAc                                                          <  1:368447/150‑1 (MQ=255)
                                                                                      gCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCt                                                      <  1:259570/150‑1 (MQ=255)
                                                                                          gggCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCt                                                                                                                                         <  1:451896/63‑1 (MQ=255)
                                                                                          ggacTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCt                                                                                                                                         >  2:451896/4‑63 (MQ=255)
                                                                                             cTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCGCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATgcg                                               >  2:214637/1‑150 (MQ=255)
                                                                                              tATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATgcgc                                              >  1:417554/1‑150 (MQ=255)
                                                                                              tATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGCAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATgcgc                                              <  2:734690/150‑1 (MQ=255)
                                                                                                  ggccaggccaTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAggcggc                                                                               <  2:915054/113‑1 (MQ=255)
                                                                                                  ggccaggccaTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAggcggc                                                                               >  1:915054/1‑113 (MQ=255)
                                                                                                           aTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTAcc                                          >  2:242146/1‑141 (MQ=255)
                                                                                                           aTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTAcc                                          <  1:242146/141‑1 (MQ=255)
                                                                                                           aTGACGCACAGCGCGACGCCAAGCTCGACCCCTAAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCgg                                 >  1:533201/1‑150 (MQ=255)
                                                                                                             gACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATAAAAGGCGGCTCGCACCTGTGTTCCGCGGACCTCGTCATGCGCTACCGCCCCTCGGCa                               <  2:852503/150‑1 (MQ=255)
                                                                                                              aCGCACAGCGCGACGCCAAGCTCGACCCCTCAGGGCCGGGCCTGAgg                                                                                                                                     >  2:72677/1‑46 (MQ=25)
                                                                                                              aCGCACAGCGCGACGCCAAGCTCGACCCCTCAGGGCCGGGCCTGAgg                                                                                                                                     <  1:72677/47‑2 (MQ=25)
                                                                                                              aCGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCGTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCAc                              <  1:1030248/150‑1 (MQ=255)
                                                                                                              aCGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCAc                              >  2:1068643/1‑150 (MQ=255)
                                                                                                                  acaGCGCGACGCCAAGCTCGACCCCTCATGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCca                          >  2:581489/1‑150 (MQ=255)
                                                                                                                  acaGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCca                          >  2:1131368/1‑150 (MQ=255)
                                                                                                                  acaGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCca                          <  2:826899/150‑1 (MQ=255)
                                                                                                                            gCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCcagccgcagccg                >  1:434867/1‑150 (MQ=255)
                                                                                                                              cAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGgt              <  1:461441/150‑1 (MQ=255)
                                                                                                                                  cTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTc          <  1:732606/150‑1 (MQ=255)
                                                                                                                                       tcacTCAGGGCCGGGCCTGAGGGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGCGTTCCGCCCAGAGCAGCACGCACTAGCTCCCCTCGGCACGCCAGCGGCCGCAGGTGTTAATgac     >  2:1030566/4‑150 (MQ=255)
                                                                                                                                       ccccTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATgac     <  1:1089702/150‑1 (MQ=255)
                                                                                                                                         ccTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTAACGCGCCTCGGCACGCCAGCCGCAGCCGGTGTTCATgacga   >  1:337222/1‑150 (MQ=255)
                                                                                                                                         acTCAGGGCCGGGCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATgacga   <  2:39872/149‑1 (MQ=255)
                                                                                                                                          cTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATgacgac  >  1:1081339/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TTTTACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGAC  >  NC_002947/3791304‑3791591

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑