Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,424,484 T→C 8.8% L28P (CTG→CCG)  PP_4764 → histone deacetylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,424,4840TC 8.8% 121.5 / 2.4 46L28P (CTG→CCG) PP_4764histone deacetylase family protein
Reads supporting (aligned to +/- strand):  ref base T (25/16);  new base C (2/2);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTACTCCCATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGAC  >  NC_002947/5424337‑5424629
                                                                                                                                                   |                                                                                                                                                  
aGGTACTCCCATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCctg                                                                                                                                                   <  7:45590/149‑1 (MQ=255)
  gTACTCCCATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCctgct                                                                                                                                                 <  1:20922/149‑1 (MQ=255)
     cTCCCATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGGCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCtgctgc                                                                                                                                               <  7:28640/148‑1 (MQ=255)
        ccATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCAGATCTACCACGAAGATTACAGCCCGCAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGa                                                                                                                                           >  5:78089/1‑149 (MQ=255)
          aTTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCACGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGacc                                                                                                                                         >  8:303150/1‑149 (MQ=255)
               aGTCACGCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACGACGTAGATTACAGCCCGGATTTCACCGCTGAACACCGCTTCCCGATGTACAAGTTCCGCCTTCTGCACGACCACCt                                                                                                                                     <  1:298636/148‑1 (MQ=255)
                        tCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGc                                                                                                                          >  7:240705/1‑149 (MQ=255)
                           tGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCggg                                                                                                                       >  7:72528/1‑149 (MQ=255)
                           tGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCCGCTGC‑CGACCAGGTGGTCGTGCAGCagg                                                                                                                       <  8:55527/149‑1 (MQ=255)
                                aGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGTAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGac                                                                                                                  <  1:19705/149‑1 (MQ=255)
                                aGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGac                                                                                                                  <  2:63643/149‑1 (MQ=255)
                                 gATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGccc                                                                                                                 >  8:178284/1‑149 (MQ=255)
                                     aGCACTGGCCCTGATTGCGACCCCATCCATGACGTTGCCGCTGAGCTCCCACGAAGATTACAGCCCGGAGTGCCCCGCTGAACACCGCTTCCCGATGGACACGTTACGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCg                                                                                                             >  6:329099/1‑149 (MQ=255)
                                                  aTTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCctgct                                                                                                <  7:303150/149‑1 (MQ=255)
                                                         ccccATCCATGCGGTTGCCGCTGATCTACCACGACGATTACCGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCGTGCTGCGCCCg                                                                                         <  6:78089/149‑1 (MQ=255)
                                                         ccccATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCctgctg                                                                                               <  3:55956/143‑1 (MQ=255)
                                                         ccccATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCctgctg                                                                                               >  4:55956/1‑143 (MQ=255)
                                                         ccccATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCg                                                                                         <  8:72528/149‑1 (MQ=255)
                                                                 atGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTACACGCTGAACACCTCTTCCCTATGCACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACAAGGCCCTCCTGCGCCCGGACATCTg                                                                                 >  8:101698/1‑149 (MQ=255)
                                                                                   tACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGCCCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGCCATCCTCGc                                                               >  2:318328/1‑149 (MQ=255)
                                                                                   tACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGc                                                               >  6:134375/1‑149 (MQ=255)
                                                                                    aCCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGcc                                                              >  5:262328/1‑149 (MQ=255)
                                                                                               tACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGATGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGa                                                   >  5:305482/1‑149 (MQ=255)
                                                                                                  aGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCg                                                >  4:227231/1‑149 (MQ=255)
                                                                                                       ggAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGc                                            >  7:262493/1‑148 (MQ=255)
                                                                                                           ttCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACAt                                       >  6:297255/1‑149 (MQ=255)
                                                                                                             ccccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATc                                      <  3:227231/148‑1 (MQ=255)
                                                                                                              cccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCCACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATc                                      >  7:231575/1‑147 (MQ=255)
                                                                                                               ccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGc                                            >  3:234015/1‑140 (MQ=255)
                                                                                                               ccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGc                                            <  4:234015/140‑1 (MQ=255)
                                                                                                                cGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgc                                  <  1:318328/149‑1 (MQ=255)
                                                                                                                             ttCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCg                                     >  4:283132/1‑133 (MQ=255)
                                                                                                                             ttCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCg                                     <  3:283132/133‑1 (MQ=255)
                                                                                                                                ccGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGa                  <  8:240705/149‑1 (MQ=255)
                                                                                                                                 cGATGGACAAGTTCCGCCAGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAg                 >  2:264438/1‑149 (MQ=255)
                                                                                                                                      gACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTc            >  1:186147/1‑149 (MQ=255)
                                                                                                                                        cAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGaa                        >  2:8309/1‑135 (MQ=255)
                                                                                                                                        cAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGaa                        <  1:8309/135‑1 (MQ=255)
                                                                                                                                        cAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgc          <  7:178284/149‑1 (MQ=255)
                                                                                                                                         aaGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGCCTGACCACCGACCACGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCCCCCTGGCCCATGACCGCAGCTACATCCAGCCCGACATGAACGGCGAGCTGTcgcg         >  8:51732/1‑149 (MQ=255)
                                                                                                                                           gTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgcgcg       <  5:297255/149‑1 (MQ=255)
                                                                                                                                           gTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATCAACGGCGAGCTGTcgcgc        >  4:306935/1‑148 (MQ=255)
                                                                                                                                           gTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGCGCTGACCACCGACCAGGACCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCAAGACCCCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgagc        >  6:46582/1‑148 (MQ=255)
                                                                                                                                            ttCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgcgc        >  6:277998/1‑147 (MQ=255)
                                                                                                                                              ccGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAg     >  2:187361/1‑148 (MQ=255)
                                                                                                                                               cGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAgg    <  5:277998/148‑1 (MQ=255)
                                                                                                                                                gCCCGCTGC‑CGACCAGGTGGTCGTGCAGCAGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCt               <  6:69501/136‑1 (MQ=38)
                                                                                                                                                gCCCGCTGC‑CGACCAGGTGGTCGTGCAGCAGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCt               >  5:69501/1‑136 (MQ=38)
                                                                                                                                                gCCCGCTGC‑CGACCAGGTGGTCGTGCAGCAGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGAc  >  3:149048/1‑149 (MQ=255)
                                                                                                                                                  ctgctgCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGAc  >  3:247893/1‑147 (MQ=255)
                                                                                                                                                   tgctgcACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgc          <  6:306892/138‑1 (MQ=255)
                                                                                                                                                   tgctgcACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgc          >  5:306892/1‑138 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGGTACTCCCATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGAC  >  NC_002947/5424337‑5424629

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: ‑