Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,139,622 T→G 7.4% Q167H (CAA→CAC)  PP_5481 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,139,6220TG 7.4% 164.1 / 4.5 54Q167H (CAA→CAC) PP_5481hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/26);  new base G (2/2);  total (26/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 = 9.88e-01

TTCCGCTACGGCATTCTAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCAC  >  NC_002947/2139475‑2139766
                                                                                                                                                   |                                                                                                                                                
ttCCGCTACGGCATTCTAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCAtt                                                                                                                                                 >  2:166403/1‑149 (MQ=255)
    gCTACGGCATTCTAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGc                                                                                                                                              >  2:173314/1‑148 (MQ=255)
         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                                                                                                           >  1:282777/1‑145 (MQ=255)
                                          aCGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAg                                                                                                           <  2:282777/145‑1 (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                                                       <  8:122814/138‑1 (MQ=255)
                                                                                                     gCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCAcacca                                                       >  7:122814/1‑138 (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                                                                                  <  1:259475/90‑1 (MQ=255)
                                                                                                                          gCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCaa                                                                                  >  2:259475/1‑90 (MQ=255)
                                                                                                                            gggCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAg                                      <  1:256047/132‑1 (MQ=255)
                                                                                                                            gggCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAg                                      >  2:256047/1‑132 (MQ=255)
                                                                                                                                   aaGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTg               <  6:234848/148‑1 (MQ=255)
                                                                                                                                     gCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTc            <  1:276994/149‑1 (MQ=255)
                                                                                                                                            aCCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGt                        <  1:33562/130‑1 (MQ=255)
                                                                                                                                            aCCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGt                        >  2:33562/1‑130 (MQ=255)
                                                                                                                                            aCCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATc     <  1:166403/149‑1 (MQ=255)
                                                                                                                                               ggCAGTGGCGTCGTCGCCAATGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCAc  >  2:309998/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TTCCGCTACGGCATTCTAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCAC  >  NC_002947/2139475‑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: ‑