Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,397,589 C→T 7.3% G9G (GGC→GGT PP_3002 → shikimate dehydrogenase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,397,5890CT 7.3% 151.3 / 2.0 55G9G (GGC→GGTPP_3002shikimate dehydrogenase‑like protein
Reads supporting (aligned to +/- strand):  ref base C (23/28);  new base T (2/2);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

CAAGTAGCCCCGGACCTTGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCAT  >  NC_002947/3397444‑3397738
                                                                                                                                                 |                                                                                                                                                     
cAAGTAGCCCCGGACCTTGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGccgac                                                                                                                                                   <  2:281791/150‑1 (MQ=255)
  aGTAGCCCCGGACCTTGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGccgaccg                                                                                                                                                 >  2:621072/1‑150 (MQ=255)
        cccGGACCTTGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAAc                                                                                                                                           <  1:670545/150‑1 (MQ=255)
        cccGGACCTTGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTAGGATTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAAc                                                                                                                                           <  1:965768/150‑1 (MQ=255)
                ttGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCaa                                                                                                                                   >  2:739184/1‑150 (MQ=255)
                ttGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCaa                                                                                                                                   >  2:670632/1‑150 (MQ=255)
                    ccAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTcc                                                                                                                               >  2:259620/1‑150 (MQ=255)
                           aaCGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCAATGAAGAAATCACTGCCTTGCAGGCAGATGGAGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCACCGCTGa                                                                                                                        >  1:203516/1‑150 (MQ=255)
                                cGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGCGATCGGCCGACCGATCAACCACACCGAGTCTCCGCTGATTCAt                                                                                                                   <  2:1007432/150‑1 (MQ=255)
                                      ggCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTg                                                                                                             <  2:343292/150‑1 (MQ=255)
                                        ccgccgGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGtt                                                                                                           <  1:904830/150‑1 (MQ=255)
                                          gccgGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCg                                                                                                         >  2:981149/1‑150 (MQ=255)
                                                    gcttCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCaa                                                                                                                                   <  1:238154/114‑1 (MQ=255)
                                                    gcttCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCaa                                                                                                                                   >  2:238154/1‑114 (MQ=255)
                                                        cAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCacac                                                                                                                                      >  2:598614/1‑107 (MQ=255)
                                                        cAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCacac                                                                                                                                      <  1:598614/107‑1 (MQ=255)
                                                            ccTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACcagc                                                                                       >  2:947575/1‑150 (MQ=255)
                                                              tGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACcagcag                                                                                     <  2:931309/150‑1 (MQ=255)
                                                                 agaaATCACTGCCTTGCAGGCAGATGGCTTGTTGTTTACCCACGTAGATGCACTATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTg                                                                                  <  2:203516/150‑1 (MQ=255)
                                                                     aTCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCCGTTCGCCCAGGCCAGCAACCAGCAGCTGGAAt                                                                              <  1:621072/150‑1 (MQ=255)
                                                                      tCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTg                                                                                                             >  1:492850/1‑118 (MQ=255)
                                                                      tCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTg                                                                                                             <  2:492850/118‑1 (MQ=255)
                                                                                    gCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAg                                                               <  1:1070795/150‑1 (MQ=255)
                                                                                         tGGCGTGTTGTTTACCCACGGAGATGCATGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGATGATTCATGGCCTGATCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCg                                                          >  2:583288/1‑150 (MQ=255)
                                                                                             gtgtTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTgg                                                      >  1:405718/1‑150 (MQ=255)
                                                                                              tgttgtTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGga                                                     <  1:947575/150‑1 (MQ=255)
                                                                                              tgttgtTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTAACCGCGGATTCATGTCCTGGTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGga                                                     <  1:981149/150‑1 (MQ=255)
                                                                                                  gttTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGgacgac                                                 >  1:1064823/1‑150 (MQ=255)
                                                                                                   tttACCCACGGAGATGCACGATGAGCGACCGCTCCCCAGTGAGCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACt                                                <  2:88080/150‑1 (MQ=255)
                                                                                                        ccACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCCGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGaa                                           <  1:815208/150‑1 (MQ=255)
                                                                                                        ccACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGaa                                           >  2:1033381/1‑150 (MQ=255)
                                                                                                          aCGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCa                                                                                                >  1:105721/1‑95 (MQ=255)
                                                                                                          aCGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCa                                                                                                <  2:105721/95‑1 (MQ=255)
                                                                                                                aTGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCa                                      <  2:1075140/147‑1 (MQ=255)
                                                                                                                aTGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCa                                      >  1:1075140/1‑147 (MQ=255)
                                                                                                                 tGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTa                                  >  1:181017/1‑150 (MQ=255)
                                                                                                                  gCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTAc                                 >  2:498595/1‑150 (MQ=255)
                                                                                                                    aCGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTg                               <  2:1089505/150‑1 (MQ=255)
                                                                                                                     cGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGc                              <  2:349371/150‑1 (MQ=255)
                                                                                                                      gATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCa                             >  1:822701/1‑150 (MQ=255)
                                                                                                                             gACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGt                      <  1:498595/150‑1 (MQ=255)
                                                                                                                               ccGCTACGCAGTGATCCGCCGTCCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAg                    <  2:255647/150‑1 (MQ=255)
                                                                                                                                cGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGc                   <  1:583288/150‑1 (MQ=255)
                                                                                                                                cGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGc                   <  1:1033381/150‑1 (MQ=255)
                                                                                                                                cGCTACGCAGTCATCGGCCGGCGGAGCAACCACACGAAGGCCCCGCGGATTGATGGCCTGTTCTCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTACCGGCCGCTCGTTGGACGACCCCTAAGCCCAGGTACTACAGATTCGTAGc                   >  2:618890/1‑150 (MQ=255)
                                                                                                                                  cTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGa                 >  1:323366/1‑150 (MQ=255)
                                                                                                                                  cTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGa                 >  1:786390/1‑150 (MQ=255)
                                                                                                                                        gttGATCGGTCGGCCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGaa                >  1:232203/3‑145 (MQ=255)
                                                                                                                                        gttGATCGGTCGGCCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGaa                <  2:232203/143‑1 (MQ=255)
                                                                                                                                        aGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACgg                                                                          >  2:429951/1‑87 (MQ=255)
                                                                                                                                        aGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACgg                                                                          <  1:429951/87‑1 (MQ=255)
                                                                                                                                        aGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGcc                                        >  2:352212/1‑121 (MQ=255)
                                                                                                                                        aGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGcc                                        <  1:352212/121‑1 (MQ=255)
                                                                                                                                         ttGATCGGTCGGCCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCg                                                          >  2:645252/2‑102 (MQ=255)
                                                                                                                                         ttGATCGGTCGGCCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCg                                                          <  1:645252/101‑1 (MQ=255)
                                                                                                                                            aTCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAg       <  1:739184/150‑1 (MQ=255)
                                                                                                                                              cGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGacgca     <  1:588175/150‑5 (MQ=255)
                                                                                                                                                 ccgaccgaTCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCAt  >  1:50174/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CAAGTAGCCCCGGACCTTGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCAT  >  NC_002947/3397444‑3397738

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: