Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,793,145 G→A 7.4% V69V (GTG→GTA PP_2445 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,793,1450GA 7.4% 156.3 / 2.5 55V69V (GTG→GTAPP_2445membrane protein
Reads supporting (aligned to +/- strand):  ref base G (23/27);  new base A (2/2);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

ATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCCTGCTGTGGAAGGCGACCAAGGAAATTCACGAAAACGTCGACCCGCAT  >  NC_002947/2793005‑2793286
                                                                                                                                            |                                                                                                                                             
aTGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGAt                                                                                                                                       >  2:300341/1‑149 (MQ=255)
   gAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGt                                                                                                                                    >  7:78219/1‑149 (MQ=255)
      tttGTAGTGGGCATCGAAAACCGTATTTTCTTTTCTTTCTTTACAAAAAATTTCTTGGTGGAGTTTTGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCa                                                                                                                                 <  1:259847/148‑1 (MQ=255)
            cTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTgac                                                                                                                           <  3:105901/149‑1 (MQ=255)
            cTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTga                                                                                                                            <  3:232795/148‑1 (MQ=255)
             tGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTgacc                                                                                                                          >  7:200990/1‑149 (MQ=255)
                 cATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTgaccgacc                                                                                                                      >  5:41389/1‑149 (MQ=255)
                       cAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGgtgt                                                                                                                >  2:194911/1‑149 (MQ=255)
                          ccTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTAGCCTGGATCGTCCAGTTGACCGACCCGGTGTTc                                                                                                              >  2:268030/1‑148 (MQ=255)
                           cTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCg                                                                                                             >  8:243038/1‑148 (MQ=255)
                              aTCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAg                                                                                                          <  2:179625/148‑1 (MQ=255)
                              aTCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTAGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGaa                                                                                                           <  1:268030/147‑1 (MQ=255)
                                   cATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTgac                                                                                                                           <  8:178778/126‑1 (MQ=255)
                                   cATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTgac                                                                                                                           >  7:178778/1‑126 (MQ=255)
                                   cATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTc                                                                                                     <  1:300341/148‑1 (MQ=255)
                                    aTCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCgg                                                                                                   >  1:155983/1‑149 (MQ=255)
                                           tCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTAGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTc                                                                                                     <  6:1442/140‑1 (MQ=255)
                                           tCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTAGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTc                                                                                                     >  5:1442/1‑140 (MQ=255)
                                             cTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCtt                                                                                          <  1:194911/149‑1 (MQ=255)
                                                  caacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCt                                                                                                                                           >  7:80811/1‑95 (MQ=255)
                                                  caacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCt                                                                                                                                           <  8:80811/95‑1 (MQ=255)
                                                  caacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGt                                                                                     >  1:247159/1‑149 (MQ=255)
                                                  caacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGt                                                                                     <  6:41389/149‑1 (MQ=255)
                                                  caacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGt                                                                                     >  8:242975/1‑149 (MQ=255)
                                                  caacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGt                                                                                     >  6:137692/1‑149 (MQ=255)
                                                   aacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATc                                                                                                                                      <  3:230751/99‑1 (MQ=255)
                                                   aacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATc                                                                                                                                      >  4:230751/1‑99 (MQ=255)
                                                   aacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTg                                                                                    >  5:298299/1‑149 (MQ=255)
                                                    acaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTg                                                                                    >  6:88606/1‑148 (MQ=255)
                                                            ccGGTGGAGTACCGCTCCAAGGCGCGGCGTATTGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGAtgt                                                                           >  3:199189/1‑149 (MQ=255)
                                                            ccGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGAtg                                                                            <  6:298299/148‑1 (MQ=255)
                                                                             cAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCggtg                                                          <  2:247159/149‑1 (MQ=255)
                                                                             cAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCggtg                                                          <  8:104780/149‑1 (MQ=255)
                                                                             cAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATACTGATTGCCggt                                                           <  8:78219/148‑1 (MQ=255)
                                                                              aaGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGAt                                                                             <  8:174542/129‑1 (MQ=255)
                                                                              aaGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGAt                                                                             >  7:174542/1‑129 (MQ=255)
                                                                               aGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCggtggt                                                        <  7:243038/149‑1 (MQ=255)
                                                                                   gcggcgTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGt                                                    <  5:137692/149‑1 (MQ=255)
                                                                                    cggcgTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTTTGtt                                                   <  3:223392/149‑1 (MQ=255)
                                                                                     ggcgTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTc                                                  <  2:155983/149‑1 (MQ=255)
                                                                                           tCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCctgct                                             <  5:88606/148‑1 (MQ=255)
                                                                                               tATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTg                                                                                    <  8:91676/105‑1 (MQ=255)
                                                                                               tATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTg                                                                                    >  7:91676/1‑105 (MQ=255)
                                                                                                 tCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCCTGCTGTGGaa                                       <  7:149350/148‑1 (MQ=255)
                                                                                                    gCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCAGTGGTCTGTTCCTGCTGTGGAAGGc                                    >  8:64323/1‑148 (MQ=255)
                                                                                                     cATGGCCTTGGTCATGCGACTGGCCGTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCCTGCTGTGGAAGGCGa                                  <  8:256626/149‑1 (MQ=255)
                                                                                                      aTGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCTGTGGTCTGTTCCTGCTGTGGAAGGCGAc                                 <  7:64323/149‑1 (MQ=255)
                                                                                                              ggTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCCTGCTGTGGAAGGCGACCAAGGaaa                         >  7:12950/1‑149 (MQ=255)
                                                                                                              ggTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCCTGCTGTGGAAGGCGACCAAGGa                           <  8:200990/147‑1 (MQ=255)
                                                                                                                  aTGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCCTGCTGTGGAAGGCGACCAAGGAAATTCa                     >  5:135997/1‑149 (MQ=255)
                                                                                                                  aTGCGACTGGCCTTGCTCAGCACCGTGGCCCGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCCTGCTGTGGAAGGCGACCAAGGAAATTCa                     >  8:55322/1‑149 (MQ=255)
                                                                                                                    gCGACTGGCCATGCTCAGCACCGTGGCCTGGATCGTCAAGTTGACCGACCCGGTGTTCGAAGTGTTCTGCAACGCCATCTCGTGGAAGGATCTGATCCTGAATGCAGGTGGTCTGTTCCTGCTGTCGAAGGCGCCAAAGGAAATTc                      >  2:109692/1‑146 (MQ=255)
                                                                                                                            ccTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCCTGCTGTGGAAGGCGACCAAGGAAATTCACGAAAACGTc           <  4:199189/149‑1 (MQ=255)
                                                                                                                                 cTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCCTGCTGTGGAAGGCGACCAAGGAAATTCACGAAAACGTCGAccc      <  6:135997/149‑1 (MQ=255)
                                                                                                                                     gCACCGTCGCCTGGATCGTCCAGTTGACCGACACGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGCAGGTCTGTTCCTGCTGTGGAAGGCGACCAAGGAAATACACGAAAACGTCGACCCGCAt  >  4:34427/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                             
ATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCCTGCTGTGGAAGGCGACCAAGGAAATTCACGAAAACGTCGACCCGCAT  >  NC_002947/2793005‑2793286

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: