Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 882,022 G→T 8.5% A512S (GCC→TCC)  PP_0765 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947882,0220GT 8.5% 173.9 / 3.0 55A512S (GCC→TCC) PP_0765hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/25);  new base T (3/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGTGAAGTGGGTGTTACCCATGTCGGCGGCCTGGAGAATGCCCACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGCCGTGTCATCTGGGACTACAACGACGTATTCGCGGGTATCAACCTGAAGCCGAGCGTG  >  NC_002947/881879‑882171
                                                                                                                                               |                                                                                                                                                     
tCGGTGAAGTGGGTGTTACCCATGTCGGCGGCCTGGAGAATGCCCACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATcgccggcg                                                                                                                                                 <  1:406404/150‑1 (MQ=255)
         tGGGTGTTACCCATGTCGGCGGCCTGGAGAATGCCCACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGCTACCTGCGAGGAGCTCTACGGCAGCACCAGCTCCGGCGCCGGTGCgg                                                                                                                                        >  2:270807/1‑150 (MQ=255)
          gggTGTTACCCATGTCGGCGGCCTGGAGAATGCCCACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTTAATGGCAGCACCATCGCCGGCGCCGGTGCggg                                                                                                                                       >  2:369604/1‑150 (MQ=255)
          gggTGTTACCCATGTCGGCGGCCTGGAGAATGCCCACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCggc                                                                                                                                       <  2:515561/150‑2 (MQ=255)
           ggTGTTACCCATGTCGGCGGCCTGGAGAATGCCCACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCgggg                                                                                                                                      >  2:677873/1‑150 (MQ=255)
                tACCCATGTCGGCGGCCTGGAGAATGACCAGGTTACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTc                                                                                                                                 <  2:94358/150‑1 (MQ=255)
                   ccATGTCGGCGGCCTGGAGAATGCCCACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCAc                                                                                                                              >  1:247664/1‑150 (MQ=255)
                   ccATGTCGGCGGCCTGGAGAATGCCCACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACATGCCAGGCGCGCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCAc                                                                                                                              >  2:2421/1‑150 (MQ=255)
                     aTGTCGGCGGCCTGGAGAATGCCCACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCg                                                                                                                            >  2:120277/1‑150 (MQ=255)
                     aTGTCGGCGGCCTGGAGAATGCCCACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCg                                                                                                                            <  1:677873/150‑1 (MQ=255)
                                  ggAGAATGCCCACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTGCa                                                                                                                               <  1:751859/134‑1 (MQ=255)
                                  ggAGAATGCCCACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCa                                                                                                                               >  2:751859/1‑134 (MQ=255)
                                   gagaATGCCCACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGc                                                                                                              >  1:528810/1‑150 (MQ=255)
                                      aaTGCCCACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCaaa                                                                                                           <  1:270807/150‑1 (MQ=255)
                                          cccACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCg                                                                                                       >  2:68705/1‑150 (MQ=255)
                                            cACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCgggg                                                                                                                                      >  2:572814/1‑117 (MQ=255)
                                            cACGATACCCGCTACGGTCGCGAACCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCgggg                                                                                                                                      <  1:572814/117‑1 (MQ=255)
                                                aTACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGCCCACCGCCAACCTGAACCGCAAATGCGAGAacg                                                                                                 >  1:391504/1‑150 (MQ=255)
                                                     cGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGc                                                                                            <  1:1936/150‑1 (MQ=255)
                                                     cGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGc                                                                                            <  1:120277/150‑1 (MQ=255)
                                                     cGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCATCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGc                                                                                            <  1:448030/150‑1 (MQ=255)
                                                             ccgcgACCGGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTAcaccac                                                                                    <  2:104649/149‑1 (MQ=255)
                                                               gcgACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACAccacca                                                                                  <  2:525713/150‑1 (MQ=255)
                                                                    ccGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACt                                                                             >  2:187922/1‑150 (MQ=255)
                                                                     cGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACAccacca                                                                                  >  1:525652/1‑144 (MQ=255)
                                                                     cGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACAccacca                                                                                  <  2:525652/144‑1 (MQ=255)
                                                                         gtTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTacac                                                                                       <  1:258326/135‑1 (MQ=255)
                                                                         gtTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTacac                                                                                       >  2:258326/1‑135 (MQ=255)
                                                                         gtTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGCCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTg                                                                        <  2:391504/150‑1 (MQ=255)
                                                                           tCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTggg                                                                      <  1:460085/150‑1 (MQ=255)
                                                                           tCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTggg                                                                      <  1:68705/150‑1 (MQ=255)
                                                                                   ggTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTGAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACcgc                                                              >  1:335813/1‑150 (MQ=255)
                                                                                    gTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCATCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACcgcg                                                             <  1:546394/150‑1 (MQ=255)
                                                                                        gTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCTCGGCCg                                                         >  2:448539/1‑150 (MQ=255)
                                                                                        gTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGCCg                                                         >  1:388277/1‑150 (MQ=255)
                                                                                             cATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCTCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACaccacc                                                                                   >  1:567591/1‑119 (MQ=255)
                                                                                             cATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCTCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACaccacc                                                                                   <  2:567610/119‑1 (MQ=255)
                                                                                             cATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCTCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACaccacc                                                                                   <  2:567591/119‑1 (MQ=255)
                                                                                             cATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCTCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACaccacc                                                                                   >  1:567610/1‑119 (MQ=255)
                                                                                              aTCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAAc                                                                                                                       <  2:598599/82‑1 (MQ=255)
                                                                                              aTCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAAc                                                                                                                       >  1:598599/1‑82 (MQ=255)
                                                                                                ccACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGCCgtg                                                       <  2:678712/144‑1 (MQ=255)
                                                                                                ccACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGCCgtg                                                       >  1:678712/1‑144 (MQ=255)
                                                                                                 cACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGc                                                           >  2:342256/1‑139 (MQ=255)
                                                                                                 cACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGc                                                           <  1:342256/139‑1 (MQ=255)
                                                                                                  aCCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCTCGGCCGTGTCATCTgg                                               <  1:448539/150‑1 (MQ=255)
                                                                                                          cgcCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGCCGTGTCATCTGGGACTACaa                                       <  1:18685/150‑1 (MQ=255)
                                                                                                            ccGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGCCGTGTCATCTGGGACTACAacg                                     >  2:208757/1‑150 (MQ=255)
                                                                                                              gATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGCCGTGTCATCTGGGACTACAacgac                                   <  1:208757/150‑1 (MQ=255)
                                                                                                                      ccAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGCCGTGTCATCTGGGACTACAACGACGTATTcgc                           >  2:382084/1‑150 (MQ=255)
                                                                                                                        aGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGCCGTGTCATCTGGGACTACAACGACGTATTCGCgg                         <  2:247664/150‑1 (MQ=255)
                                                                                                                         ggCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGCCGTGTCATCTGGGACTACAACGACGTATTCGCggg                        >  2:418932/1‑150 (MQ=255)
                                                                                                                                        caCCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGCCGTGTCATCTgg                                               <  2:138856/112‑1 (MQ=255)
                                                                                                                                        caCCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGCCGTGTCATCTgg                                               >  1:138856/1‑112 (MQ=255)
                                                                                                                                           cATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGGACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGGCGTGTCATCTGGGACTACAACGACGTATTCGCGGGTATCAACCTGAAGCCGAg      >  1:716768/1‑150 (MQ=255)
                                                                                                                                             tCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGATACCGCGGCCGTATTAGTGTGGACTACAACGACGTATTCGCGGGTATCAACCTGAAGCCGAGCg    >  1:356019/1‑150 (MQ=255)
                                                                                                                                               gccggcgccggTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGCCGTGTCATCTGGGACTACAACGACGTATTCGCGGGTATCAACCTGAAGCCGAGCGTg  <  2:528810/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
TCGGTGAAGTGGGTGTTACCCATGTCGGCGGCCTGGAGAATGCCCACGATACCCGCTACGGTCGCGACCCGGTGTTCGGCCCAGGTCCGTTGCCATCCACCGGTGGCGCCGATACCTGCCAGGCGCTCAATGGCAGCACCATCGCCGGCGCCGGTGCGGGGGCGTCCACCGCCAACCTGAACCGCAAATGCGAGAACGACGGCTACACCACCAGCACTTCCTGGGGTTACCGCGGCCGTGTCATCTGGGACTACAACGACGTATTCGCGGGTATCAACCTGAAGCCGAGCGTG  >  NC_002947/881879‑882171

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