Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,424,545 C→G 7.9% 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,5450CG 7.9% 141.5 / 2.5 46H32Q (CAC→CAGPP_4764histone deacetylase family protein
Reads supporting (aligned to +/- strand):  ref base C (23/19);  new base G (2/2);  total (25/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 = 1.00e+00

GTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGC  >  NC_002947/5424401‑5424692
                                                                                                                                                |                                                                                                                                                    
gTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGt                                                                                                                                                 >  1:579305/1‑150 (MQ=255)
 tCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTc                                                                                                                                                >  2:428015/1‑150 (MQ=255)
     ccGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCa                                                                                                                                           >  2:47648/1‑150 (MQ=255)
     ccGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCa                                                                                                                                           <  2:644458/150‑1 (MQ=255)
          aTGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCggg                                                                                                                                      <  1:266320/150‑1 (MQ=255)
           tGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGc                                                                                                                                     <  1:424511/150‑1 (MQ=255)
           tGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGc                                                                                                                                     >  1:47254/1‑150 (MQ=255)
           tGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGc                                                                                                                                     <  2:659194/150‑1 (MQ=255)
           tGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCCGCTGCCGACCAGGTGGTCGTGCAGCAGGc                                                                                                                                     <  1:330139/150‑1 (MQ=255)
                    cAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGa                                                                                                                           >  1:750634/1‑150 (MQ=255)
                            gCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACACGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGcc                                                                                                                    >  2:107913/1‑150 (MQ=255)
                            gCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGcc                                                                                                                    >  1:109622/1‑150 (MQ=255)
                               cTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCctg                                                                                                                 <  2:47254/150‑1 (MQ=255)
                               cTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCctg                                                                                                                 <  2:109622/150‑1 (MQ=255)
                                tGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGATGATCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCctgc                                                                                                                <  1:650390/150‑1 (MQ=255)
                                 gATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCctgct                                                                                                               <  2:82197/150‑1 (MQ=255)
                                             atccatGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACAt                                                                                                   <  1:142444/150‑1 (MQ=255)
                                                         gCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGCCAAGATCCGCCTGCCGCCAAGCTAGCTATGCTTCATCGGTCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTCCCCCAACGa                                                                                       >  1:61425/1‑150 (MQ=255)
                                                           cGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACa                                                                                     >  2:272461/1‑150 (MQ=255)
                                                             cTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCg                                                                                                                            >  1:804374/1‑109 (MQ=255)
                                                             cTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCg                                                                                                                            <  2:804374/109‑1 (MQ=255)
                                                              tGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGa                                                                                                      >  1:666790/1‑130 (MQ=255)
                                                              tGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGa                                                                                                      <  2:666790/130‑1 (MQ=255)
                                                                aTCTACCACGAAGATTTTATCCCTGAGGGCCTGTCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTc                                                                                <  1:189534/150‑1 (MQ=255)
                                                                     ccACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCCGCTGCCGACCAGGTGGTCgtgc                                                                                                                                            >  1:250048/1‑82 (MQ=37)
                                                                     ccACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCCGCTGCCGACCAGGTGGTCgtgc                                                                                                                                            <  2:250048/85‑4 (MQ=37)
                                                                        cGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGc                                                                        >  2:142125/1‑150 (MQ=255)
                                                                            gATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCAt                                                                    <  1:272461/150‑1 (MQ=255)
                                                                              ttACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGa                                                                  <  1:428015/150‑1 (MQ=255)
                                                                                  aGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGc                                                              >  1:364931/1‑150 (MQ=255)
                                                                                          gTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCTACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACAt                                                      <  2:750634/150‑1 (MQ=255)
                                                                                              cccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCAGGGGCCAGTCCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAg                                                  >  2:23112/1‑150 (MQ=255)
                                                                                               ccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgc                                                 >  2:232590/1‑150 (MQ=255)
                                                                                                 gCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgcgc                                               >  1:704381/1‑150 (MQ=255)
                                                                                                 gCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgcgc                                               >  1:715398/1‑150 (MQ=255)
                                                                                                                ccGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGcc                                                                                                                    <  2:94714/66‑1 (MQ=255)
                                                                                                                ccGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGcc                                                                                                                    >  1:94714/1‑66 (MQ=255)
                                                                                                                  gATGGACAAGTTCCACCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAAc                                                                                         <  2:724334/91‑1 (MQ=255)
                                                                                                                  gATGGACAAGTTCCACCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAAc                                                                                         >  1:724334/1‑91 (MQ=255)
                                                                                                                       aCAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgc                         >  2:131788/1‑150 (MQ=255)
                                                                                                                         aaGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATGTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgcgc                       >  1:274249/1‑150 (MQ=255)
                                                                                                                         aaGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGc                                                                        >  1:842685/1‑101 (MQ=255)
                                                                                                                         aaGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGc                                                                        <  2:842685/101‑1 (MQ=255)
                                                                                                                            ttCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTc                           >  2:482826/1‑143 (MQ=255)
                                                                                                                            ttCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTc                           <  1:482826/143‑1 (MQ=255)
                                                                                                                                       gcACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCCGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGcgccg         >  1:539247/1‑150 (MQ=255)
                                                                                                                                              caccTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGc  <  2:274249/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGC  >  NC_002947/5424401‑5424692

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