Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,211,029 G→A 11.5% D307D (GAC→GAT)  PP_4592 ← metal dependent hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,211,0290GA11.5% 112.9 / 3.0 43D307D (GAC→GAT) PP_4592metal dependent hydrolase
Reads supporting (aligned to +/- strand):  ref base G (20/18);  new base A (3/2);  total (23/20)
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

CAAGATCTGGTCATGGCATTCGCCATCGGCGCACAGCTGCGGGTCCTTGAGGTTGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGC  >  NC_002947/5210882‑5211172
                                                                                                                                                   |                                                                                                                                               
caAGATCTGGTCATGGCATTCGCCATCGGCGCACAGCTGCGGGTCCTTGAGGTTGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGt                                                                                                                                                >  4:191846/1‑149 (MQ=255)
        ggTCATGGCATTCGCCATCGGCGCACAGCTGCGGGTCCTTGAGGTTGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCTGGACATCCtgcagc                                                                                                                                        >  7:123782/1‑149 (MQ=255)
             tGGCATTCGCCATCGGCGCACAGCTGCGGGTCCTTGAGGTTGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTc                                                                                                                                   <  7:179892/149‑1 (MQ=255)
               gCATTCGCCATCGGCGCACAGCTGCGGGTCCTTGAGGTTGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCTGGACATCCtgcagcgtgccg                                                                                                                                  <  8:123782/148‑7 (MQ=255)
               gCATTCGCCATCGGCGCACAGCTGCGGGTCCTTGAGGTTGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTcgc                                                                                                                                 <  4:108909/149‑1 (MQ=255)
                cATTCGCCATCGGCGCACAGCTGCGGGTCCTTGAGGATGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTcgcg                                                                                                                                >  1:234091/1‑149 (MQ=255)
                            gcgcACAGCTGCGGGTCCTTGAGGTTGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTAc                                                                                                                    >  3:93021/1‑149 (MQ=255)
                                        gggTCCTTGAGGTTGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTa                                                                                                        <  3:191846/149‑1 (MQ=255)
                                          gTCCTTGAGGTTGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCTGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGt                                                                                                      >  6:253120/1‑149 (MQ=255)
                                            ccTTGAGGTTGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTcc                                                                                                    >  2:172587/1‑149 (MQ=255)
                                              ttGAGGTTGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCTGGACATCCtgcagcgt                                                                                                                                      >  6:188469/1‑111 (MQ=255)
                                              ttGAGGTTGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCTGGACATCCtgcagcgt                                                                                                                                      <  5:188469/113‑3 (MQ=255)
                                                gAGGTTGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCtt                                                                                                  <  4:93021/147‑1 (MQ=255)
                                                 aGGTTGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTgc                                                                                                <  2:234091/148‑1 (MQ=255)
                                                               gATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAg                                                                                                       <  6:35863/127‑1 (MQ=255)
                                                               gATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAg                                                                                                       >  5:35863/1‑127 (MQ=255)
                                                                 tGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCgg                                                                                                                >  1:189691/1‑116 (MQ=255)
                                                                 tGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCgg                                                                                                                <  2:189691/116‑1 (MQ=255)
                                                                          cgcgATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCgg                                                                                                                <  6:80260/107‑1 (MQ=255)
                                                                          cgcgATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCgg                                                                                                                >  5:80260/1‑107 (MQ=255)
                                                                            cgATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCTGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCAc                                                                    <  5:253120/149‑1 (MQ=255)
                                                                                       gCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAAccgccg                                                         >  6:245729/1‑149 (MQ=255)
                                                                                          tttGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGc                                                                                       <  8:57572/116‑1 (MQ=255)
                                                                                          tttGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGc                                                                                       >  7:57572/1‑116 (MQ=255)
                                                                                                 cccGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTAcc                                               <  8:107159/149‑1 (MQ=255)
                                                                                                         gTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGt                                       >  1:303668/1‑149 (MQ=255)
                                                                                                         gTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGt                                       >  3:71510/1‑149 (MQ=255)
                                                                                                             tCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCAc                                   >  7:299404/1‑149 (MQ=255)
                                                                                                                 gCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCg                               >  4:31388/1‑149 (MQ=255)
                                                                                                                    aCTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGaa                            <  2:303668/149‑1 (MQ=255)
                                                                                                                      tGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGaaaa                          >  7:303862/1‑149 (MQ=255)
                                                                                                                         aCGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGc                       <  8:299404/149‑1 (MQ=255)
                                                                                                                           ggCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCca                     <  4:71510/149‑1 (MQ=255)
                                                                                                                             ccAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTGATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTg                                      <  4:145595/130‑1 (MQ=255)
                                                                                                                                   cGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGccc             >  8:97614/1‑149 (MQ=255)
                                                                                                                                    ggCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCcag                    <  7:177783/141‑1 (MQ=255)
                                                                                                                                    ggCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCcag                    >  8:177783/1‑141 (MQ=255)
                                                                                                                                       aCGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCt         <  1:172587/149‑1 (MQ=255)
                                                                                                                                         gctgctGGACATCCTGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCACCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTcc       >  4:271472/1‑149 (MQ=255)
                                                                                                                                             cAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGa                                         <  3:181573/111‑1 (MQ=255)
                                                                                                                                             cAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGa                                         >  4:181573/1‑111 (MQ=255)
                                                                                                                                             cAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAAACGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACCCCAGCACGCCCTGCTCCCtcg   >  4:339275/1‑147 (MQ=255)
                                                                                                                                               ggATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGc  >  8:249588/1‑148 (MQ=255)
                                                                                                                                                   gTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGc  <  3:91581/144‑1 (MQ=255)
                                                                                                                                                   gTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGc  >  4:91581/1‑144 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CAAGATCTGGTCATGGCATTCGCCATCGGCGCACAGCTGCGGGTCCTTGAGGTTGCTGACGTCGATGAACTGCACGCGATCGCAGGTGCCTTTGCAGCCCGACTGGTTGTCACGCCACTGCACGGCCAGGCCGGCACGCTGCAGGATGTCCAGCAGCCCTTCGCGGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCCGAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGC  >  NC_002947/5210882‑5211172

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