Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,424,504:1 +T 9.6% coding (103/915 nt) 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,5041.T 9.6% 187.9 / 16.2 52coding (103/915 nt)PP_4764histone deacetylase family protein
Reads supporting (aligned to +/- strand):  ref base . (27/20);  new base T (3/2);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

TCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGCCTGC  >  NC_002947/5424361‑5424648
                                                                                                                                                |                                                                                                                                                
tCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGc                                                                                                                                             >  7:240705/1‑149 (MQ=255)
   tGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCggg                                                                                                                                          >  7:72528/1‑149 (MQ=255)
   tGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCCGCTGCCGACCAGGTGGTCGTGCAGCagg                                                                                                                                          <  8:55527/149‑1 (MQ=255)
        aGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGTAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGac                                                                                                                                     <  1:19705/149‑1 (MQ=255)
        aGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGac                                                                                                                                     <  2:63643/149‑1 (MQ=255)
         gATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGccc                                                                                                                                    >  8:178284/1‑149 (MQ=255)
             aGCACTGGCCCTGATTGCGACCCCATCCATGACGTTGCCGCTGAGCTCCCACGAAGATTACAGCCCGGAGTGCCCCGCTGAACACCGCTTCCCGATGGACACGTTACGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCg                                                                                                                                >  6:329099/1‑149 (MQ=255)
                          aTTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCctgct                                                                                                                   <  7:303150/149‑1 (MQ=255)
                                 ccccATCCATGCGGTTGCCGCTGATCTACCACGACGATTACCGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCGTGCTGCGCCCg                                                                                                            <  6:78089/149‑1 (MQ=255)
                                 ccccATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCctgctg                                                                                                                  >  4:55956/1‑143 (MQ=255)
                                 ccccATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCctgctg                                                                                                                  <  3:55956/143‑1 (MQ=255)
                                 ccccATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCg                                                                                                            <  8:72528/149‑1 (MQ=255)
                                         atGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTACACGCTGAACACCTCTTCCCTATGCACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACAAGGCCCTCCTGCGCCCGGACATCTg                                                                                                    >  8:101698/1‑149 (MQ=255)
                                                           tACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGCCCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGCCATCCTCGc                                                                                  >  2:318328/1‑149 (MQ=255)
                                                           tACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGc                                                                                  >  6:134375/1‑149 (MQ=255)
                                                            aCCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGcc                                                                                 >  5:262328/1‑149 (MQ=255)
                                                                       tACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGATGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGa                                                                      >  5:305482/1‑149 (MQ=255)
                                                                          aGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCg                                                                   >  4:227231/1‑149 (MQ=255)
                                                                               ggAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGc                                                               >  7:262493/1‑148 (MQ=255)
                                                                                   ttCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACAt                                                          >  6:297255/1‑149 (MQ=255)
                                                                                     ccccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATc                                                         <  3:227231/148‑1 (MQ=255)
                                                                                      cccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCCACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATc                                                         >  7:231575/1‑147 (MQ=255)
                                                                                       ccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGc                                                               >  3:234015/1‑140 (MQ=255)
                                                                                       ccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGc                                                               <  4:234015/140‑1 (MQ=255)
                                                                                        cGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgc                                                     <  1:318328/149‑1 (MQ=255)
                                                                                                     ttCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCg                                                        >  4:283132/1‑133 (MQ=255)
                                                                                                     ttCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCg                                                        <  3:283132/133‑1 (MQ=255)
                                                                                                        ccGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGa                                     <  8:240705/149‑1 (MQ=255)
                                                                                                         cGATGGACAAGTTCCGCCAGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAg                                    >  2:264438/1‑149 (MQ=255)
                                                                                                              gACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTc                               >  1:186147/1‑149 (MQ=255)
                                                                                                                cAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGaa                                           <  1:8309/135‑1 (MQ=255)
                                                                                                                cAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGaa                                           >  2:8309/1‑135 (MQ=255)
                                                                                                                cAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgc                             <  7:178284/149‑1 (MQ=255)
                                                                                                                 aaGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGCCTGACCACCGACCACGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCCCCCTGGCCCATGACCGCAGCTACATCCAGCCCGACATGAACGGCGAGCTGTcgcg                            >  8:51732/1‑149 (MQ=255)
                                                                                                                   gTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgcgcg                          <  5:297255/149‑1 (MQ=255)
                                                                                                                   gTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATCAACGGCGAGCTGTcgcgc                           >  4:306935/1‑148 (MQ=255)
                                                                                                                   gTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGCGCTGACCACCGACCAGGACCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCAAGACCCCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgagc                           >  6:46582/1‑148 (MQ=255)
                                                                                                                    ttCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgcgc                           >  6:277998/1‑147 (MQ=255)
                                                                                                                      ccGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAg                        >  2:187361/1‑148 (MQ=255)
                                                                                                                       cGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAgg                       <  5:277998/148‑1 (MQ=255)
                                                                                                                        gCCCGCTGCCGACCAGGTGGTCGTGCAGCAGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCt                                  <  6:69501/136‑1 (MQ=38)
                                                                                                                        gCCCGCTGCCGACCAGGTGGTCGTGCAGCAGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCt                                  >  5:69501/1‑136 (MQ=38)
                                                                                                                        gCCCGCTGCCGACCAGGTGGTCGTGCAGCAGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGAc                     >  3:149048/1‑149 (MQ=255)
                                                                                                                          ctgctgCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGAc                     >  3:247893/1‑147 (MQ=255)
                                                                                                                           tgctgcACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgc                             >  5:306892/1‑138 (MQ=255)
                                                                                                                           tgctgcACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgc                             <  6:306892/138‑1 (MQ=255)
                                                                                                                               tgccGACCAGGTGGTCGTGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAgcgc               >  6:240378/4‑149 (MQ=255)
                                                                                                                               gcACGACCACCTTGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGcgcc              <  6:152981/149‑1 (MQ=255)
                                                                                                                                   gACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGc                                       >  4:254995/1‑120 (MQ=255)
                                                                                                                                   gACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGc                                       <  3:254995/120‑1 (MQ=255)
                                                                                                                                       accTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGc      <  5:46582/149‑1 (MQ=255)
                                                                                                                                            gTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGCCTGc  <  5:329099/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGCCTGC  >  NC_002947/5424361‑5424648

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: