Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,318,758 G→A 7.7% T93T (ACC→ACT)  treSA ← trehalose synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,318,7580GA 7.7% 148.7 / 2.6 52T93T (ACC→ACT) treSAtrehalose synthase A
Reads supporting (aligned to +/- strand):  ref base G (25/23);  new base A (2/2);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

TGCGCAGCGGCCATGCGGCTCATCTGGATCAGTTCCTGCTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGC  >  NC_002947/3318612‑3318903
                                                                                                                                                  |                                                                                                                                                 
tgCGCAGCGGCCATGCGGCTCATCTGGATCAGTTCCTGCTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGgt                                                                                                                                                 >  5:336844/1‑149 (MQ=255)
 gcgcAGCGGCCATGCGGCTCATCTGGATCAGTTCCTGCTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGgtg                                                                                                                                                >  4:242956/1‑149 (MQ=255)
         gCCATGCGGCTCATCTGGATCAGTTCCTGCTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGcc                                                                                                                                        <  2:262728/149‑1 (MQ=255)
             tGCGGCTCATCTGGATCAGTTCCTGCTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGt                                                                                                                                    <  4:180152/149‑1 (MQ=255)
                 gCTCATCTGGATCAGTTCCTGCTCGCTGCCGTTCAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGc                                                                                                                                <  4:97120/149‑1 (MQ=255)
                     aTCTGGATCAGTTCCTGCTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGa                                                                                                                             >  5:160514/1‑148 (MQ=255)
                          gATCAGTTCCTGCTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTcc                                                                                                                        <  6:75775/148‑1 (MQ=255)
                          gATCAGTTCCTGCTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTcc                                                                                                                        >  5:268540/1‑148 (MQ=255)
                          gATCAGTTCCTGCTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGACGGTGTGCAGGCCCTGTACGCCGATCTcc                                                                                                                        >  3:294816/1‑148 (MQ=255)
                           aTCAGTTCCTGCTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGa                                                                                                                      >  3:24829/1‑149 (MQ=255)
                           aTCAGTTCCTGCTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGa                                                                                                                      >  3:47216/1‑149 (MQ=255)
                            tCAGTTCCTGCTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGa                                                                                                                      >  2:35130/1‑148 (MQ=255)
                                   cTGCTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCt                                                                                                              >  7:95799/1‑149 (MQ=255)
                                    tGCTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGCCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCtt                                                                                                             >  6:193760/1‑149 (MQ=255)
                                      cTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTc                                                                                                            <  5:202722/148‑1 (MQ=255)
                                           tGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCt                                                                                                              >  5:133087/1‑141 (MQ=255)
                                           tGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCt                                                                                                              <  6:133087/141‑1 (MQ=255)
                                                    gggcgTCCATGTCGAAGCTGATGCGGTCGAAGTGGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTgg                                                                                             <  8:145017/145‑1 (MQ=255)
                                                        gTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCa                                                                                         <  4:24829/149‑1 (MQ=255)
                                                        gTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCCACAGGCGCTTCCACAACGCTTCGTGGGCCa                                                                                         <  1:35130/149‑1 (MQ=255)
                                                           gATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCca                                                                                                          >  4:204408/1‑129 (MQ=255)
                                                           gATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCca                                                                                                          <  3:204408/129‑1 (MQ=255)
                                                           gATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAgg                                                                                       <  6:268540/148‑1 (MQ=255)
                                                            aTGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAggg                                                                                      >  5:370029/1‑148 (MQ=255)
                                                            aTGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGgc                                                                                     <  6:336844/149‑1 (MQ=255)
                                                                 gAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGgcacc                                                                                  >  1:359508/1‑147 (MQ=255)
                                                                            ggTCGAAGTTGCCGCCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCaa                                                                     <  5:193760/149‑1 (MQ=255)
                                                                                gAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCt                                                                 >  4:326581/1‑149 (MQ=255)
                                                                                  aGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCg                                                               >  6:319834/1‑149 (MQ=255)
                                                                                   gTTGCCGTCCACGCTGGGGGTGAGTGCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCgg                                                              <  4:347536/149‑1 (MQ=255)
                                                                                    ttGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCggg                                                             <  3:242956/149‑1 (MQ=255)
                                                                                    ttGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCggg                                                             <  5:319834/149‑1 (MQ=255)
                                                                                         gTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGT‑GGCCAGGGCACCGAGCACCGAGCAACACTCGGGCGCgat                                                       >  6:12796/1‑149 (MQ=255)
                                                                                            cACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCgat                                                       <  1:372505/147‑1 (MQ=255)
                                                                                              cGCTGGGGGTGAGTTCGTGGCCGCGGATGCCACCGGAGAGTTTGATCCGGTCGGTGTGGAGGCCCTGTACGCCTATCTCCGATTGGCGCTTCCACATCGCTTCGTTGGCCATGGAACCGTGCACCGAGGAAGCCTCGTGCGGGattatg                                                   <  4:294816/149‑1 (MQ=255)
                                                                                              cGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGatgatg                                                   <  4:47216/149‑1 (MQ=255)
                                                                                              cGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGatgatg                                                   <  6:342583/149‑1 (MQ=255)
                                                                                                      gTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCg                                             >  1:216655/1‑147 (MQ=255)
                                                                                                      gTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCg                                             <  2:216655/147‑1 (MQ=255)
                                                                                                       tGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCgggg                                          >  2:352687/1‑149 (MQ=255)
                                                                                                         aGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTa                                        <  8:241348/149‑1 (MQ=255)
                                                                                                              gcgGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCg                                   >  5:360223/1‑149 (MQ=255)
                                                                                                                 gCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCa                                <  6:370029/149‑1 (MQ=255)
                                                                                                                            cACCGGACAGTTTGATCGGGCCAGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCgatgat                                                    <  5:162018/118‑1 (MQ=255)
                                                                                                                            cACCGGACAGTTTGATCGGGCCAGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCgatgat                                                    >  6:168020/1‑118 (MQ=255)
                                                                                                                            cACCGGACAGTTTGATCGGGCCAGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCgatgat                                                    >  6:162018/1‑118 (MQ=255)
                                                                                                                            cACCGGACAGTTTGATCGGGCCAGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCgatgat                                                    <  5:168020/118‑1 (MQ=255)
                                                                                                                                  aCAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGAtt               >  7:45338/1‑149 (MQ=255)
                                                                                                                                  aCAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGAtt               >  3:90516/1‑149 (MQ=255)
                                                                                                                                      tttGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCgg           >  1:290989/1‑149 (MQ=255)
                                                                                                                                              ggCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACAGTCAGCCACACCGAGGCGATTTCGGTGGCgcgg   >  8:143696/1‑149 (MQ=255)
                                                                                                                                               gCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCgcggc  >  7:366202/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGCGCAGCGGCCATGCGGCTCATCTGGATCAGTTCCTGCTCGCTGCCGTACAGTGGGTCGATGTCGAAGCTGATGCGGTCGAAGTTGCCGTCCACGCTGGGGGTGAGTTCGCGGCCGCGGATGCCACCGGACAGTTTGATCGGGCCGGTGTGCAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGC  >  NC_002947/3318612‑3318903

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: ‑