Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,424,552:1 +T 12.1% 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,5521.T12.1% 113.2 / 12.0 33coding (103/915 nt)PP_4764histone deacetylase family protein
Reads supporting (aligned to +/- strand):  ref base . (20/9);  new base T (2/2);  total (22/11)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

CAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGCCTGCCCT  >  NC_002947/5424415‑5424699
                                                                                                                                          |                                                                                                                                                   
cAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGa                                                                                                                                         >  2:385325/1‑150 (MQ=255)
                ccTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCct                                                                                                                         >  2:468375/1‑150 (MQ=255)
                         gACCCCATCAATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGccc                                                                                                                >  1:609949/1‑150 (MQ=255)
                                        gTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTCACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCaa                                                                                                 <  1:344873/150‑1 (MQ=255)
                                             cGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACa                                                                                            >  2:474386/1‑150 (MQ=255)
                                              gCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCgg                                                                                                                                                    <  1:101833/93‑1 (MQ=255)
                                              gCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCgg                                                                                                                                                    >  2:101833/1‑93 (MQ=255)
                                                 gATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGaccac                                                                                                                                     >  1:20261/1‑105 (MQ=255)
                                                 gATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGaccac                                                                                                                                     <  2:20261/105‑1 (MQ=255)
                                                   tCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCg                                                                                      <  1:385325/150‑1 (MQ=255)
                                                         aCGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTgg                                                                                >  2:482937/1‑150 (MQ=255)
                                                              gATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCAt                                                                           >  1:301851/1‑150 (MQ=255)
                                                                       ccGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGc                                                                  >  1:301919/1‑150 (MQ=255)
                                                                                cccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAg                                                         >  1:421301/1‑150 (MQ=255)
                                                                                cccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAg                                                         >  2:237001/1‑150 (MQ=255)
                                                                                cccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAg                                                         >  2:287997/1‑150 (MQ=255)
                                                                                  cGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCAGGTGGTCGTGCAGcag                                                                                                                                              >  1:352705/1‑62 (MQ=25)
                                                                                  cGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCAGGTGGTCGTGCAGcag                                                                                                                                              <  2:352705/64‑3 (MQ=25)
                                                                                            cGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAAc                                             >  1:540467/1‑150 (MQ=255)
                                                                                                tCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCg                                         >  2:460261/1‑150 (MQ=255)
                                                                                                 cccGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGa                                        <  1:482937/150‑1 (MQ=255)
                                                                                                    gATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCt                                     >  2:570460/1‑150 (MQ=255)
                                                                                                           aaGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgcgc                              >  2:575583/1‑150 (MQ=255)
                                                                                                             gTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGa                            <  2:301919/150‑1 (MQ=255)
                                                                                                                  gCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCa                                                                            <  2:547419/97‑1 (MQ=255)
                                                                                                                  gCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCa                                                                            >  1:547419/1‑97 (MQ=255)
                                                                                                                  gCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGAcc                       >  1:220439/1‑150 (MQ=255)
                                                                                                                     tgctgcACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAgc                    <  1:287997/150‑1 (MQ=255)
                                                                                                                     tgctgcACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAgc                    <  2:553526/150‑1 (MQ=255)
                                                                                                                               ccaccTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTTGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATTAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTccg          >  2:234789/1‑148 (MQ=255)
                                                                                                                                caccTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGc         <  1:221018/150‑1 (MQ=255)
                                                                                                                                  ggtGGTCGTGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGc                                                                                                      <  2:429769/54‑1 (MQ=255)
                                                                                                                                  ggtGGTCGTGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGc                                                                                                      >  1:429769/3‑56 (MQ=255)
                                                                                                                                    tGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGCCTGc     >  2:415410/1‑150 (MQ=255)
                                                                                                                                       tCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGCCTGCCCt  >  1:276357/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
CAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGCCTGCCCT  >  NC_002947/5424415‑5424699

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