Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,951,650 G→A 8.8% L455L (CTG→TTG)  ppk ← polyphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,951,6500GA 8.8% 156.3 / 5.4 57L455L (CTG→TTG) ppkpolyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (26/26);  new base A (3/2);  total (29/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 = 7.61e-01

CGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGGCGCACG  >  NC_002947/5951505‑5951798
                                                                                                                                                 |                                                                                                                                                    
cgggcgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggcgg                                                                                                                                                  <  2:963762/150‑1 (MQ=255)
cgggcgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggcgg                                                                                                                                                  >  2:565299/1‑150 (MQ=255)
     gTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACAGCACCAGCCGCCTGCAAGCGGCTggt                                                                                                                                             >  2:616864/1‑149 (MQ=255)
       tGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCATCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCa                                                                                                                                           <  1:861601/150‑1 (MQ=255)
           ggCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCCTGCAAGCGGCTggtgcggta                                                                                                                                       <  1:616864/150‑8 (MQ=255)
               tGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAgg                                                                                                                                   <  1:277764/150‑1 (MQ=255)
                 gTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGtt                                                                                                                                 <  2:897911/150‑1 (MQ=255)
                  tAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTg                                                                                                                                <  1:964731/150‑1 (MQ=255)
                    gTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGa                                                                                                                              >  1:1250572/1‑150 (MQ=255)
                     ttGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGAc                                                                                                                             <  1:473949/150‑1 (MQ=255)
                                 aGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCaaaa                                                                                                                 >  2:1018703/1‑150 (MQ=255)
                                    tGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTTGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAAcgc                                                                                                              <  2:757651/150‑1 (MQ=255)
                                                  gATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGAt                                                                                                <  1:559962/150‑1 (MQ=255)
                                                  gATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGAt                                                                                                <  2:481067/150‑1 (MQ=255)
                                                    tctcGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTCGAGCCCGCTACACCGTAGATAACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCa                                                                                              >  1:517276/1‑150 (MQ=255)
                                                      tcGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCAcc                                                                                            <  2:525321/150‑1 (MQ=255)
                                                        gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTCGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGc                                                                                          >  1:515174/1‑150 (MQ=255)
                                                        gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGc                                                                                                                                      >  2:587770/1‑106 (MQ=255)
                                                        gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGc                                                                                                                                      <  1:587770/106‑1 (MQ=255)
                                                        gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGc                                                                                          >  2:808454/1‑150 (MQ=255)
                                                         cccTGCTCGCCGCGCAGAATCAGCATCATCTTTGCGTGGGTCTTGAAGCCAACCACACCCTAAATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTTGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCa                                                                                         <  1:16582/150‑1 (MQ=255)
                                                          ccTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAg                                                                                        >  1:1157742/1‑150 (MQ=255)
                                                          ccTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTAGAATATTCGTGAAGACTCGTGCGGAACTCGAGCACCGCAg                                                                                        >  2:976761/1‑150 (MQ=255)
                                                           cTGCTCGCGCCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGt                                                                                       <  1:1018703/150‑1 (MQ=255)
                                                               tCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACACAACACTATGTCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAcc                                                                                   >  2:314631/1‑150 (MQ=255)
                                                                 gcggcgCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTc                                                                                                                        >  1:209105/1‑111 (MQ=255)
                                                                 gcggcgCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTc                                                                                                                        <  2:209105/111‑1 (MQ=255)
                                                                                    atcTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAg                                                                                        >  1:1058859/1‑124 (MQ=255)
                                                                                    atcTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAg                                                                                        <  2:1058859/124‑1 (MQ=255)
                                                                                     tcTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAAcgc                                                                                                              <  1:524651/101‑1 (MQ=255)
                                                                                     tcTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAAcgc                                                                                                              >  2:524651/1‑101 (MQ=255)
                                                                                     tcTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAt                                                             <  1:128228/150‑1 (MQ=255)
                                                                                            gTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTTTCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAg                                                      <  1:656078/150‑1 (MQ=255)
                                                                                                 tCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGt                                                 >  1:713133/1‑150 (MQ=255)
                                                                                                     gAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCAcc                                                                                            >  1:95227/1‑103 (MQ=255)
                                                                                                     gAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCAcc                                                                                            <  2:95227/103‑1 (MQ=255)
                                                                                                       aGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTcgacga                                           >  1:853121/1‑150 (MQ=255)
                                                                                                           aaCCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATccacca                                                       <  1:934761/134‑1 (MQ=255)
                                                                                                           aaCCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATccacca                                                       >  2:934761/1‑134 (MQ=255)
                                                                                                           aaCCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTc                                       >  2:312097/1‑150 (MQ=255)
                                                                                                             ccACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCg                                     >  2:641687/1‑150 (MQ=255)
                                                                                                             ccACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTAGATAACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCg                                     >  2:174480/1‑150 (MQ=255)
                                                                                                              cacaCCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGa                                    <  2:1250572/150‑1 (MQ=255)
                                                                                                              cacaCCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGa                                    >  2:987328/1‑150 (MQ=255)
                                                                                                              accacCGTGGATCTCTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGa                                    >  2:690749/3‑150 (MQ=255)
                                                                                                                 acCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGtt                                 >  2:459145/1‑150 (MQ=255)
                                                                                                                     tAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCAGTTACGCGCCGCATcc                                                           >  2:330481/1‑120 (MQ=255)
                                                                                                                     tAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCAGTTACGCGCCGCATcc                                                           <  1:330481/120‑1 (MQ=255)
                                                                                                                      aGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGccc                            <  2:464219/150‑1 (MQ=255)
                                                                                                                                cGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTa                  >  2:418304/1‑150 (MQ=255)
                                                                                                                                  cACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACa                >  2:1145815/1‑150 (MQ=255)
                                                                                                                                    ccAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAgg              >  1:795991/1‑150 (MQ=255)
                                                                                                                                     cAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTTAACACGCGCAAGCAACTGGATCACCGCAGTGACATCCTTGCCGTTACGCGCCGCATCCATCCGGGCTTCGACGATTTTTGATTTAACACACGAGCGGTACactg             >  2:225456/1‑147 (MQ=255)
                                                                                                                                     cAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCa                                                                                         >  1:1032542/1‑74 (MQ=39)
                                                                                                                                     cAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCa                                                                                         <  2:1032542/74‑1 (MQ=39)
                                                                                                                                       gCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTc           <  1:418304/150‑1 (MQ=255)
                                                                                                                                         cGCTTGCAGGCGGCTGGCGATCTGCAGGTTGGACTCTTCGTCACAACTCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCACCGCATCCACCAGGGCGTCGACCATTTCCGCGTTAGCCCCCGAGCGGTACAGGGTCTg         <  2:517276/150‑1 (MQ=255)
                                                                                                                                          gCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTgg        <  2:1124439/150‑1 (MQ=255)
                                                                                                                                                aGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGGCGCACg  <  2:1157742/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGGCGCACG  >  NC_002947/5951505‑5951798

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