Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,139,628 A→T 7.7% R165R (CGT→CGA)  PP_5481 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,139,6280AT 7.7% 150.7 / 2.9 52R165R (CGT→CGA) PP_5481hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/26);  new base T (2/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

GGCATTCTAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCAC  >  NC_002947/2139484‑2139766
                                                                                                                                                |                                                                                                                                          
ggCATTCTAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCgacgac                                                                                                                                        >  2:43808/1‑149 (MQ=255)
    ttCTAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCAGTGGCGTCGTCGcc                                                                                                                                     >  2:215685/1‑148 (MQ=255)
     tCTAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCAc                                                                                                                                   >  7:149742/1‑149 (MQ=255)
          gggTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCg                                                                                                                              >  8:102974/1‑149 (MQ=255)
          gggTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCg                                                                                                                              <  2:260240/149‑1 (MQ=255)
          gggTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCAGTGGCGTCGTCGCCAATGCCg                                                                                                                              <  1:215685/149‑1 (MQ=255)
          gggTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCAGTGGCGTCGTCGCCAATGCCg                                                                                                                              <  5:76759/149‑1 (MQ=255)
              gtCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGccc                                                                                                                          <  4:43216/149‑1 (MQ=255)
                  gATATCAGCATTGGCACGTTGATCACCTTCATGCTGGTGGGCTCTGCGGTGCTCTACGTTGTTTGCTTTCGCTTGCTGGTCTGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGTcccaccc                                                                                                                      <  4:131860/149‑1 (MQ=255)
                  gATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGcccaccc                                                                                                                      <  1:43808/149‑1 (MQ=255)
                        aGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAgc                                                                                                                <  2:173197/149‑1 (MQ=255)
                         gCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAgcg                                                                                                               <  3:3192/149‑1 (MQ=255)
                            ttGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCGCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCa                                                                                                            <  1:161152/149‑1 (MQ=255)
                                 aCGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAg                                                                                                           <  2:282777/145‑1 (MQ=255)
                                 aCGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAg                                                                                                           >  1:282777/1‑145 (MQ=255)
                                    gctaTCACCTTCATGCTGGTTGCCACTGGTGTGCTATACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGACTGAAGCCTTGCACCGGCATTGGCGACGACGCCAATGCCGGCCCACCCACCAGCGCCAGTAGAACa                                                                                                    <  7:103893/146‑1 (MQ=255)
                                           ccTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCACGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACAtggtg                                                                                               >  3:311138/1‑147 (MQ=255)
                                            cTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACAtggtg                                                                                               >  1:314464/1‑146 (MQ=255)
                                             ttCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGAAGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGgag                                                                                            >  3:207966/1‑148 (MQ=255)
                                                 tGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAg                                                                                         >  7:83664/1‑147 (MQ=255)
                                                      gtggCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGc                                                                                    >  2:109306/1‑147 (MQ=255)
                                                       tggCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCACGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAg                                                                                 <  4:311138/149‑1 (MQ=255)
                                                            aCTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGc                                                                            >  6:222156/1‑149 (MQ=255)
                                                             cTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGcc                                                                           <  2:314464/149‑1 (MQ=255)
                                                             cTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGcc                                                                           <  8:149742/149‑1 (MQ=255)
                                                              tGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGcc                                                                           >  2:276994/1‑148 (MQ=255)
                                                                  gTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGca                                                                      >  3:153804/1‑149 (MQ=255)
                                                                        tACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGACCAGCACCAcgc                                                                >  1:123672/1‑149 (MQ=255)
                                                                        tACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGCCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCAcgc                                                                >  2:40576/1‑149 (MQ=255)
                                                                            ttgttTGCGGTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGca                                                            <  4:11941/149‑1 (MQ=255)
                                                                              gttTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGccc                                                                          <  6:222696/133‑1 (MQ=255)
                                                                              gttTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGAGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGccc                                                                          >  5:222696/1‑133 (MQ=255)
                                                                                    gTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACaccacc                                                     >  5:166592/1‑148 (MQ=255)
                                                                                     ttCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCAccga                                                   <  7:102974/149‑1 (MQ=255)
                                                                                            gCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCAcacca                                                       >  7:122814/1‑138 (MQ=255)
                                                                                            gCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCAcacca                                                       <  8:122814/138‑1 (MQ=255)
                                                                                             cTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCAccgagccgag                                             >  3:21445/1‑147 (MQ=255)
                                                                                             cTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCa                                           <  8:83664/149‑1 (MQ=255)
                                                                                                   aGCGGCCGCTGCATGCGGGACTGAAGCCTTGCACCTGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAgg                                     <  4:207966/149‑1 (MQ=255)
                                                                                                      ggCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGtt                                  <  5:222156/149‑1 (MQ=255)
                                                                                                        ccGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGc                                >  5:234848/1‑149 (MQ=255)
                                                                                                          gcgGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCca                              >  4:106213/1‑149 (MQ=255)
                                                                                                            ggCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCcaca                            <  6:166592/149‑1 (MQ=255)
                                                                                                                 gCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCaa                                                                                  >  2:259475/1‑90 (MQ=255)
                                                                                                                 gCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCaa                                                                                  <  1:259475/90‑1 (MQ=255)
                                                                                                                   gggCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAg                                      >  2:256047/1‑132 (MQ=255)
                                                                                                                   gggCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAg                                      <  1:256047/132‑1 (MQ=255)
                                                                                                                          aaGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTg               <  6:234848/148‑1 (MQ=255)
                                                                                                                            gCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTc            <  1:276994/149‑1 (MQ=255)
                                                                                                                                   aCCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGt                        >  2:33562/1‑130 (MQ=255)
                                                                                                                                   aCCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGt                        <  1:33562/130‑1 (MQ=255)
                                                                                                                                   aCCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATc     <  1:166403/149‑1 (MQ=255)
                                                                                                                                      ggCAGTGGCGTCGTCGCCAATGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCAc  >  2:309998/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                          
GGCATTCTAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCAC  >  NC_002947/2139484‑2139766

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