Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,317,213 C→G 9.2% P608P (CCG→CCC treSA ← trehalose synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,317,2130CG 9.2% 164.1 / 6.7 54P608P (CCG→CCCtreSAtrehalose synthase A
Reads supporting (aligned to +/- strand):  ref base C (22/27);  new base G (3/2);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

CGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCGGGCAC  >  NC_002947/3317069‑3317358
                                                                                                                                                |                                                                                                                                                  
cGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgc                                                                                                                                                <  6:113631/149‑1 (MQ=255)
  tctcGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgctg                                                                                                                                              <  4:76275/149‑1 (MQ=255)
                  gACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCTTCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCa                                                                                                                              >  4:111655/1‑149 (MQ=255)
                  gACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTAcc                                                                                                                               <  7:209274/148‑1 (MQ=255)
                     caccGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCa                                                                                                                              >  8:17831/1‑146 (MQ=255)
                     caccGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCa                                                                                                                              <  7:17831/146‑1 (MQ=255)
                      accGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAgg                                                                                                                                   <  1:301522/140‑1 (MQ=255)
                      accGGCCCAGGTGTGAACCAGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATCCACCATCACCAGCAGCCCCCGGCTGCTAACCTCAgg                                                                                                                                   >  2:301522/1‑140 (MQ=255)
                             cAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAg                                                                                                                                    >  4:113775/1‑132 (MQ=255)
                             cAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAg                                                                                                                                    <  3:113775/132‑1 (MQ=255)
                             cAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACgg                                                                                                                   <  2:181840/149‑1 (MQ=255)
                                   tgAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCtggtg                                                                                                             <  7:106087/149‑1 (MQ=255)
                                   tgAACCCGGTCAACAGCAGTTCCTAGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCtggt                                                                                                              <  5:302862/148‑1 (MQ=255)
                                      aCCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCAcc                                                                                                                                        >  6:165617/1‑119 (MQ=255)
                                      aCCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCAcc                                                                                                                                        <  5:165617/119‑1 (MQ=255)
                                          ggTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGa                                                                                                         <  7:332845/146‑1 (MQ=255)
                                          ggTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCAACATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGa                                                                                                         >  8:332845/1‑146 (MQ=255)
                                                cagcagTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAgg                                                                                                >  1:321432/1‑149 (MQ=255)
                                                 agcagTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGc                                                                                               <  3:111655/149‑1 (MQ=255)
                                                   cagTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCtg                                                                                                                >  3:91689/1‑130 (MQ=255)
                                                   cagTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCtg                                                                                                                <  4:91689/130‑1 (MQ=255)
                                                     gTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATg                                                                                                       >  7:199704/1‑137 (MQ=255)
                                                     gTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATg                                                                                                       <  8:205067/137‑1 (MQ=255)
                                                     gTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATg                                                                                                       >  7:205067/1‑137 (MQ=255)
                                                     gTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATg                                                                                                       <  8:199704/137‑1 (MQ=255)
                                                      ttCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGgc                                                                                          >  3:43532/1‑149 (MQ=255)
                                                          tCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTg                                                                                            <  5:78235/143‑1 (MQ=255)
                                                          tCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTg                                                                                            >  6:78235/1‑143 (MQ=255)
                                                          tCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGACTGGTGGCGATGCCGTAGGCCTGGCGca                                                                                       >  8:202274/1‑148 (MQ=255)
                                                                 tCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCcagca                                                                              <  2:215677/149‑1 (MQ=255)
                                                                      tCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCAGCAGTttc                                                                         <  4:43532/149‑1 (MQ=255)
                                                                          ggCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttctt                                                                       >  7:334105/1‑148 (MQ=255)
                                                                           gCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttcttc                                                                     <  2:321432/148‑1 (MQ=255)
                                                                                        gcAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCaaa                                                       >  5:284644/1‑149 (MQ=255)
                                                                                         cAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAAt                                                      <  8:334105/149‑1 (MQ=255)
                                                                                              aTCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAAt                                                  <  8:113914/148‑1 (MQ=255)
                                                                                                       ccGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGACTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTc                                        <  7:202274/149‑1 (MQ=255)
                                                                                                             ccGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCAcagc                                                                                   <  6:70662/101‑1 (MQ=255)
                                                                                                             ccGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCAcagc                                                                                   >  5:70662/1‑101 (MQ=255)
                                                                                                                       gCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAgg                        <  6:284644/149‑1 (MQ=255)
                                                                                                                             gCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCg                  >  4:182916/1‑149 (MQ=255)
                                                                                                                               aCCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAg                                                            >  8:202621/1‑105 (MQ=255)
                                                                                                                               aCCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTc                                        >  3:37896/1‑125 (MQ=255)
                                                                                                                               aCCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTc                                        <  4:37896/125‑1 (MQ=255)
                                                                                                                               aCCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGt                 <  3:182916/148‑1 (MQ=255)
                                                                                                                                 cATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGta                >  7:239306/1‑147 (MQ=255)
                                                                                                                                     gtgagcagcCCGGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCCACCTTGCAGGCAAATGAATCACtctc                                           <  1:4616/114‑1 (MQ=255)
                                                                                                                                     gtgagcagcCCGGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACtctc                                           >  2:4616/4‑116 (MQ=255)
                                                                                                                                     aCCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAgcgc          >  4:207615/1‑149 (MQ=255)
                                                                                                                                      tgagcagcCCGGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTc                            >  3:15294/3‑130 (MQ=255)
                                                                                                                                      tgagcagcCCGGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGcc         >  6:255256/3‑149 (MQ=255)
                                                                                                                                      tcagcagCCCGGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTc                            <  4:15294/129‑1 (MQ=255)
                                                                                                                                       cagcagCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttcttc                                                                     >  5:354210/1‑88 (MQ=255)
                                                                                                                                       cagcagCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttcttc                                                                     <  6:354210/88‑1 (MQ=255)
                                                                                                                                             ccccGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCGGGCAc  >  4:250646/1‑149 (MQ=255)
                                                                                                                                             ccccGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCGGGCAc  >  2:154720/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCGGGCAC  >  NC_002947/3317069‑3317358

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: