Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,544,106 T→C S82P (TCC→CCC)  ppk → polyphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,544,1060TC98.6% 594.4 / ‑3.2 142S82P (TCC→CCC) ppkpolyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (74/66);  total (76/66)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 150 total aligned reads displayed.

AAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATGGCGCGTAACC  >  REL606/2544034‑2544178
                                                                        |                                                                        
aaGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGcccc                                                                       >  2:912024/1‑76 (MQ=255)
aaGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGcccc                                                                       >  2:1106176/1‑76 (MQ=255)
 aGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCg                                                                      >  2:1186547/1‑76 (MQ=255)
 aGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCg                                                                      >  2:5228621/1‑76 (MQ=255)
  gCGACGCATCATTTTTAGCGAAGAACAAGGCTCCAACTCTCATCCCCGCCATTTACTGGGCAAAATTCAGCCCCgg                                                                     <  2:5822997/76‑1 (MQ=255)
   cGACGCATCATTATTAGCGAAGAACCAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCggg                                                                    >  1:838454/1‑76 (MQ=255)
   cGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCggg                                                                    >  1:622741/1‑76 (MQ=255)
    gACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGt                                                                   <  2:2441181/76‑1 (MQ=255)
    gACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGt                                                                   >  1:6695314/1‑76 (MQ=255)
        catcatTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTg                                                               <  2:1422350/76‑1 (MQ=255)
        catcatTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTg                                                               <  2:3388355/76‑1 (MQ=255)
        catcatTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTg                                                               <  1:6976901/76‑1 (MQ=255)
         atcatTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGa                                                              >  2:249565/1‑76 (MQ=255)
          tcatTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGaa                                                             <  2:19801/76‑1 (MQ=255)
           catTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGaaa                                                            >  1:281455/1‑76 (MQ=255)
           catTATTAGCGAAGAACAAGGCTCCAACTCTCATCCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGaaa                                                            <  2:797111/76‑1 (MQ=255)
            attattAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAg                                                           >  2:4093075/1‑76 (MQ=255)
            attattAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAg                                                           <  1:3471650/76‑1 (MQ=255)
            attattAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAg                                                           <  1:5847133/76‑1 (MQ=255)
            attattAGCGAAGAACAAGGCTCCAACTCTCATCCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAg                                                           <  1:3828962/76‑1 (MQ=255)
             ttattaGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGc                                                          >  2:1420576/1‑76 (MQ=255)
             ttattaGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGc                                                          >  1:3441367/1‑76 (MQ=255)
             ttattaGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGc                                                          >  1:436128/1‑76 (MQ=255)
              tattaGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGcc                                                         >  1:6563396/1‑76 (MQ=255)
              tattaGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGcc                                                         >  2:699266/1‑76 (MQ=255)
              tattaGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGcc                                                         <  1:4793488/76‑1 (MQ=255)
               attaGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCg                                                        >  1:6060678/1‑76 (MQ=255)
                ttaGCAAAGAACAGGCCCCCAACTCTCATCCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGa                                                       <  1:4516483/76‑1 (MQ=255)
                 taGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGAt                                                      >  1:949327/1‑76 (MQ=255)
                 taGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGAt                                                      >  2:6936532/1‑76 (MQ=255)
                   gCGAAGAACAAGGCTCCACCTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCa                                                    <  2:4175698/76‑1 (MQ=255)
                     gaagaaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAgg                                                  >  1:1027867/1‑76 (MQ=255)
                     gaagaaCAAGGCGCCACCTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAgg                                                  <  2:3453630/76‑1 (MQ=255)
                     gaagaaCAAGCCTCCAACTCTCATCCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAgg                                                  <  2:4287226/76‑1 (MQ=255)
                      aagaaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGa                                                 >  2:3497871/1‑76 (MQ=255)
                      aagaaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGa                                                 >  1:2220051/1‑76 (MQ=255)
                      aagaaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGa                                                 >  2:651526/1‑76 (MQ=255)
                      aagaaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGa                                                 >  2:753335/1‑76 (MQ=255)
                      aagaaCAAGGCCCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGa                                                 <  2:2109059/76‑1 (MQ=255)
                       agaaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGaa                                                <  2:2282369/76‑1 (MQ=255)
                       agaaCAAGGCTCCAACTCTCATCCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGaa                                                <  1:3385052/76‑1 (MQ=255)
                        aaaCAAGGCTCCAACTCTCATCCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAAt                                               <  2:2324312/75‑1 (MQ=255)
                         aaCAGGGCTCCACCTCTCATCCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAAtt                                              <  2:2052257/76‑1 (MQ=255)
                         aaCAGGGCCCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAAtt                                              <  2:1115370/76‑1 (MQ=255)
                         aaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAAtt                                              >  1:3770415/1‑76 (MQ=255)
                         aaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAAtt                                              >  1:1071889/1‑76 (MQ=255)
                         aaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAAtt                                              >  1:6910526/1‑76 (MQ=255)
                         aaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAAtt                                              >  1:6922422/1‑76 (MQ=255)
                         aaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAAtt                                              <  1:2937522/76‑1 (MQ=255)
                         aaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAAtt                                              >  1:5715101/1‑76 (MQ=255)
                          aCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTc                                             >  1:5646307/1‑76 (MQ=255)
                           cAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCg                                            <  1:1871639/76‑1 (MQ=255)
                            aaGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGa                                           >  1:1948997/1‑76 (MQ=255)
                            aaGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGa                                           <  1:269454/76‑1 (MQ=255)
                              ggCTCCACCTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACg                                         <  2:5787782/76‑1 (MQ=255)
                              ggCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACg                                         >  1:1821711/1‑76 (MQ=255)
                               gCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACgg                                        >  2:1606569/1‑76 (MQ=255)
                                cTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGc                                       <  2:6075876/76‑1 (MQ=255)
                                cTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGc                                       >  2:1972701/1‑76 (MQ=255)
                                 tCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGcc                                      >  2:3739871/1‑76 (MQ=255)
                                  ccAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCct                                     >  1:1176091/1‑76 (MQ=255)
                                   cAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCctc                                    >  1:6739904/1‑76 (MQ=255)
                                   cAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCctc                                    >  1:4673566/1‑76 (MQ=255)
                                    aaCTCTCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCctct                                   >  1:5377696/1‑76 (MQ=255)
                                     aCTCTCATCCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTa                                  <  2:4894424/76‑1 (MQ=255)
                                      ctctCATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTAc                                 >  2:2507517/1‑76 (MQ=255)
                                       tctcATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACa                                >  2:2764600/1‑76 (MQ=255)
                                       tctcATCCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACa                                <  1:2705201/76‑1 (MQ=255)
                                        ctcATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACaa                               >  2:5949002/1‑76 (MQ=255)
                                         tcATTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAAc                              >  1:1615720/1‑76 (MQ=255)
                                           aTTCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGa                            >  2:1770251/1‑76 (MQ=255)
                                            ttCCCGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAg                           >  2:4707661/1‑76 (MQ=255)
                                              cccGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAgct                         >  1:4042091/1‑76 (MQ=255)
                                              cccGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAgct                         <  2:1406521/76‑1 (MQ=255)
                                              cccGCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAgct                         >  2:1740225/1‑76 (MQ=255)
                                                 gCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAgctgct                      >  1:3221673/1‑76 (MQ=255)
                                                 gCCATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAgctgct                      >  2:5045414/1‑76 (MQ=255)
                                                   cATTTACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAgctgctgc                    >  1:643243/1‑76 (MQ=255)
                                                     tttACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGctgctgctg                  <  1:6087710/76‑1 (MQ=255)
                                                      ttACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTgg                 >  1:1708096/1‑76 (MQ=255)
                                                      ttACTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTgg                 >  2:3590176/1‑76 (MQ=255)
                                                      ttACGGGGCAAAATTCAGCCCCGGGTGCTGAAACCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTgg                 <  1:4810027/76‑1 (MQ=255)
                                                        aCTGGGCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGgag               >  1:4601454/1‑76 (MQ=255)
                                                           gggCAAAATTCAGCCCCGGGTGTTGAAAGCCGATCAGAAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATg            <  2:723852/76‑1 (MQ=255)
                                                           gggCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATg            <  1:2426652/76‑1 (MQ=255)
                                                            ggaaaaaTTCGGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATgg           <  2:3353750/73‑1 (MQ=255)
                                                            gcaaaaaTTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATgg           <  2:4291991/73‑1 (MQ=255)
                                                             gaaaaaTTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATGgc          <  2:5642573/74‑1 (MQ=255)
                                                             gCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATGgc          <  1:4967108/76‑1 (MQ=255)
                                                             gCAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATGgc          <  2:2994355/76‑1 (MQ=255)
                                                              cAAAATTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATGgcg         >  2:5911451/1‑76 (MQ=255)
                                                                aaaTTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATGGcgcg       <  1:5456239/76‑1 (MQ=255)
                                                                 aaTTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATGGCGCGt      >  2:1356962/1‑76 (MQ=255)
                                                                  aTTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATGGCGCGTa     >  2:226115/1‑76 (MQ=255)
                                                                  aTTCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATGGCGCGTa     >  2:1561625/1‑76 (MQ=255)
                                                                   ttCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATGGCGCGTaa    <  1:6172930/76‑1 (MQ=255)
                                                                   ttCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATGGCGCGTaa    <  2:6073738/76‑1 (MQ=255)
                                                                   ttCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATGGCGCGTaa    <  2:477602/76‑1 (MQ=255)
                                                                   ttCAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATGGCGCGTaa    <  1:3731889/76‑1 (MQ=255)
                                                                     cAGCCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATGGCGCGTAAcc  >  2:6432073/1‑76 (MQ=255)
                                                                        |                                                                        
AAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTGCTGCTGGAGATGGCGCGTAACC  >  REL606/2544034‑2544178

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑