Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,554,507 A→C 7.0% V140G (GTC→GGC)  PP_4041 ← glycoside hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,554,5070AC 7.0% 170.1 / 4.5 57V140G (GTC→GGC) PP_4041glycoside hydrolase
Reads supporting (aligned to +/- strand):  ref base A (25/28);  new base C (2/2);  total (27/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 = 8.44e-01

ACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCAGGGG  >  NC_002947/4554369‑4554652
                                                                                                                                          |                                                                                                                                                 
aCCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAg                                                                                                                                         <  6:158667/149‑1 (MQ=255)
aCCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAg                                                                                                                                         <  1:240939/149‑1 (MQ=255)
  cGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGTCGCCGGcc                                                                                                                                              >  3:91301/1‑141 (MQ=255)
  cGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGTCGCCGGcc                                                                                                                                              <  4:91301/142‑2 (MQ=255)
        gTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTg                                                                                                                                 <  3:279467/149‑1 (MQ=255)
                         gCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGACAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCAcc                                                                                                                >  1:239359/1‑149 (MQ=255)
                             gcgtaaGTCCTGGGGGTGGGCCGGGGTGTAGGGAAAGCGGCAGGCGGTTTGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAgg                                                                                                            <  8:252592/145‑1 (MQ=255)
                             gCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAgg                                                                                                            >  8:30134/1‑149 (MQ=255)
                              cTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAggg                                                                                                           >  3:157157/1‑149 (MQ=255)
                                 aaGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCa                                                                                                         <  6:178607/148‑1 (MQ=255)
                                  aGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTgg                                                                                                                                >  5:54391/1‑124 (MQ=255)
                                  aGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTgg                                                                                                                                <  6:54391/124‑1 (MQ=255)
                                      cTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTcc                                                                                                                   <  4:132203/133‑1 (MQ=255)
                                      cTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTcc                                                                                                                   >  3:132203/1‑133 (MQ=255)
                                         gcgcTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGc                                                                                                <  7:30134/149‑1 (MQ=255)
                                           gcTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGa                                                                                                                                          <  6:226008/105‑1 (MQ=255)
                                           gcTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCg                                                                                              <  2:239359/149‑1 (MQ=255)
                                           gcTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGGGCGGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCg                                                                                              <  5:237709/149‑1 (MQ=255)
                                               gCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGt                                                                                          >  7:2816/1‑149 (MQ=255)
                                               gCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGt                                                                                          <  3:171762/149‑1 (MQ=255)
                                                cccGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTc                                                                                         >  1:231507/1‑149 (MQ=255)
                                                cccGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGGGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTc                                                                                         <  5:153569/149‑1 (MQ=255)
                                                 ccGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTc                                                                                         >  6:307987/1‑148 (MQ=255)
                                                      cTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAgc                                                                                    <  1:19270/148‑1 (MQ=255)
                                                       tGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAgcg                                                                                   >  7:50198/1‑148 (MQ=255)
                                                            cgAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCAt                                                                                <  5:169797/146‑1 (MQ=255)
                                                            cgAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCAt                                                                                >  6:169797/1‑146 (MQ=255)
                                                               aaGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAg                                                                          >  1:247888/1‑149 (MQ=255)
                                                                              tgttggtCGCGGGCTTGCGGGGGCCCGGTGCTTCCCAGGATCACCCGGTCCGGGCCGGCGACGGCACTGAGGGGCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAAt                                                            <  6:70578/145‑1 (MQ=255)
                                                                              tggtggTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAAt                                                            <  2:247888/148‑1 (MQ=255)
                                                                                gtggtCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATgg                                                          <  5:307987/148‑1 (MQ=255)
                                                                                    tCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTc                                                     >  1:175571/1‑149 (MQ=255)
                                                                                        gggCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCa                                                 <  4:157157/149‑1 (MQ=255)
                                                                                          gCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGcc                                                                                               <  2:327807/101‑1 (MQ=255)
                                                                                          gCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGcc                                                                                               <  2:322736/101‑1 (MQ=255)
                                                                                          gCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGcc                                                                                               >  1:322736/1‑101 (MQ=255)
                                                                                          gCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGcc                                                                                               >  1:327807/1‑101 (MQ=255)
                                                                                              gCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAg                                           >  6:7416/1‑149 (MQ=255)
                                                                                                        ggTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTcgacga                                 <  2:231507/149‑1 (MQ=255)
                                                                                                          tgctgcACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCTTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATg                               <  2:283573/149‑1 (MQ=255)
                                                                                                            ctgcACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATgc                              <  2:175571/148‑1 (MQ=255)
                                                                                                                   ggATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCaccacc                      >  2:178179/1‑149 (MQ=255)
                                                                                                                   ggATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGCGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGATAAGTTTTCGACGATGCGCaccacc                      <  8:50198/149‑1 (MQ=255)
                                                                                                                   ggATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGATCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGCTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCaccacc                      >  7:36892/1‑149 (MQ=255)
                                                                                                                    gATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGGTTCGCCAGACAAGCCTTCGACGATGCGCaccacc                      >  8:114301/1‑148 (MQ=255)
                                                                                                                         aCCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAAtt                 >  1:298373/1‑148 (MQ=255)
                                                                                                                           cGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGAc              >  7:93676/1‑149 (MQ=255)
                                                                                                                            cctCAGTGCCGTCGCCGGCCCTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAAtt                 >  7:158260/3‑145 (MQ=255)
                                                                                                                            cctCAGTGCCGTCGCCGGCCCTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAAtt                 <  8:158260/143‑1 (MQ=255)
                                                                                                                             gTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACtt            <  8:2816/149‑1 (MQ=255)
                                                                                                                                aGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGc         <  7:114301/149‑1 (MQ=255)
                                                                                                                                aGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGc         >  6:312378/1‑149 (MQ=255)
                                                                                                                                  ggCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGcc        >  8:120682/1‑148 (MQ=255)
                                                                                                                                   gCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGcc        >  6:343961/1‑147 (MQ=255)
                                                                                                                                       gCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCAgggg  >  7:44851/1‑149 (MQ=255)
                                                                                                                                         gACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGa               >  7:320013/1‑134 (MQ=255)
                                                                                                                                         gACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGa               <  8:320013/134‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
ACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCAGGGG  >  NC_002947/4554369‑4554652

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