Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,424,497 C→G 9.8% H32Q (CAC→CAG)  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,4970CG 9.8% 151.3 / 7.1 51H32Q (CAC→CAG) PP_4764histone deacetylase family protein
Reads supporting (aligned to +/- strand):  ref base C (27/19);  new base G (3/2);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

AGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGC  >  NC_002947/5424352‑5424644
                                                                                                                                                 |                                                                                                                                                    
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                                                                                                              >  4:55956/1‑143 (MQ=255)
                                          ccccATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCctgctg                                                                                                              <  3:55956/143‑1 (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)
                                                                                                                                        tgccGACCAGGTGGTCGTGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAgcgc           >  6:240378/4‑149 (MQ=255)
                                                                                                                                        gcACGACCACCTTGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGcgcc          <  6:152981/149‑1 (MQ=255)
                                                                                                                                            gACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGc                                   >  4:254995/1‑120 (MQ=255)
                                                                                                                                            gACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGc                                   <  3:254995/120‑1 (MQ=255)
                                                                                                                                                accTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGc  <  5:46582/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGC  >  NC_002947/5424352‑5424644

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