Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,791,405 C→G 10.3% 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,4050CG10.3% 111.6 / 4.8 39G283G (GGC→GGG) PP_3353sulfatase‑modifying protein
Reads supporting (aligned to +/- strand):  ref base C (23/12);  new base G (2/2);  total (25/14)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TGGTTTTACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGAC  >  NC_002947/3791258‑3791548
                                                                                                                                                   |                                                                                                                                               
tGGTTTTACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAgtcc                                                                                                                                                >  7:219548/1‑149 (MQ=255)
   ttttACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggcccgg                                                                                                                                             <  4:2398/149‑1 (MQ=255)
   ttttACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggcccgg                                                                                                                                             <  8:248565/149‑1 (MQ=255)
    tttACCCGTACATCGCCCGAGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGAAGCCCGCGGCAATGACTGGGAGATGACACGCACGGGCTATCGGCCAGGACATGACGCAAAGCGCGACGCCAAAATCGACCCCTAAggcacggc                                                                                                                                            >  8:281418/1‑149 (MQ=255)
     ttACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggcccggcc                                                                                                                                           >  6:82340/1‑149 (MQ=255)
        ccGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACTCCAAGCTCGACCCCTCAGGCCCGGCCCt                                                                                                                                         <  1:18855/148‑1 (MQ=255)
         cGGTACATCGCCCGTGGGTTGCTTCCCGTCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggcccggcacaac                                                                                                                                       <  8:238876/149‑6 (MQ=255)
         cGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGa                                                                                                                                       <  5:82340/149‑1 (MQ=255)
            tACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGGCCGGGCCTGAggg                                                                                                                                    <  2:76341/149‑3 (MQ=255)
                  gCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCTAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggcccgg                                                                                                                                             >  3:263022/1‑134 (MQ=255)
                  gCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCTAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggcccgg                                                                                                                                             <  4:263022/134‑1 (MQ=255)
                   cccGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCACTCAGGGCCGGGCCTGAggggtcga                                                                                                                               >  3:177924/1‑140 (MQ=255)
                   cccGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCACTCAGGGCCGGGCCTGAggggtcga                                                                                                                               <  4:177924/147‑8 (MQ=255)
                                 cccGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGAcc                                                                                                               >  1:308773/1‑149 (MQ=255)
                                   cGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGAcccc                                                                                                             >  7:264369/1‑149 (MQ=255)
                                    ggCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCg                                                                                                            <  1:289246/149‑1 (MQ=255)
                                             cTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCgg                                                                                                   <  2:308773/149‑1 (MQ=255)
                                                 ggCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTgg                                                                                                <  2:116569/148‑1 (MQ=255)
                                                  gCCTGTTCGATGCCGGCGGCAATGTCTGGTAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCg                                                                                              <  7:281418/149‑1 (MQ=255)
                                                               cggcggCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGGGGCGGTAATCAAAggcg                                                                                 >  6:319151/1‑149 (MQ=255)
                                                                 gcggcAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAggcggc                                                                               >  3:77129/1‑149 (MQ=255)
                                                                    gcAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCg                                                                            >  3:201386/1‑149 (MQ=255)
                                                                    gcAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCAGTGGCGGTAATCAAAGGCGGCTCg                                                                            >  4:86962/1‑149 (MQ=255)
                                                                         gTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGccccg                                                                       >  6:200439/1‑148 (MQ=255)
                                                                          tCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCt                                                                       <  4:77129/148‑1 (MQ=255)
                                                                                      cTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGAc                                                          >  6:128341/1‑149 (MQ=255)
                                                                                             gggCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGc                                                    >  4:304/1‑148 (MQ=255)
                                                                                                      gccaggccaTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTAcc                                          >  5:262365/1‑149 (MQ=255)
                                                                                                         aggccaTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGcc                                       >  6:234579/1‑149 (MQ=255)
                                                                                                              aTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCg                                  >  8:64203/1‑149 (MQ=255)
                                                                                                                         cgcgACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCcagcc                       <  4:201386/149‑1 (MQ=255)
                                                                                                                            gACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGAACCGCTGCATgcg                                               <  8:44798/122‑1 (MQ=255)
                                                                                                                            gACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGAACCGCTGCATgcg                                               >  7:44798/1‑122 (MQ=255)
                                                                                                                              cGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCcagccgcggc                  >  7:237290/1‑149 (MQ=255)
                                                                                                                               gCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTcgg                                                               <  7:189027/103‑1 (MQ=255)
                                                                                                                               gCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTcgg                                                               >  8:189027/1‑103 (MQ=255)
                                                                                                                               gCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATTCGCTACCGCCCCTCGGCACGCcagccgcagcc                 >  4:9379/1‑149 (MQ=255)
                                                                                                                               gCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCcagccgcagcc                 >  7:172255/1‑149 (MQ=255)
                                                                                                                                ccAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCCCGCcagccgcagccg                >  7:89134/1‑149 (MQ=255)
                                                                                                                                 cAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCgg               >  8:291044/1‑149 (MQ=255)
                                                                                                                                 cAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCgg               >  5:216218/1‑149 (MQ=255)
                                                                                                                                  aaGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGgt              >  4:66916/1‑149 (MQ=255)
                                                                                                                                   aGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGgtg             <  7:64203/149‑1 (MQ=255)
                                                                                                                                        tgtcacTCAGGGCCGGGCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCAt        >  1:310698/6‑149 (MQ=255)
                                                                                                                                        gACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCAt        <  3:9379/149‑1 (MQ=255)
                                                                                                                                           cgtcatggcctggcccTGAGTGACAGCTTCGAACCGGCAGTCGCCGGAGTGGCTGTGGCGGTAATCAAAGGCGGGTCGCACATGTTTTCGGCGGACCGCTGCATGCGCTAGAGCCCCTCGGCACGCCAGCGGCAGCGGTTGTTCATga      <  8:206657/142‑1 (MQ=255)
                                                                                                                                           cccTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATgac     >  5:75510/1‑149 (MQ=255)
                                                                                                                                              tCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATgacgac  >  1:254868/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TGGTTTTACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGAC  >  NC_002947/3791258‑3791548

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