Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,360,306 T→A D656E (GAT→GAA pta → phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,360,3060TA99.5% 862.4 / ‑4.2 205D656E (GAT→GAAptaphosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (125/79);  total (126/79)
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.71e-01

Only 100 of 208 total aligned reads displayed.

TGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCGGGCCGATGCTG  >  REL606/2360233‑2360378
                                                                         |                                                                        
tGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACt                                                                        <  1:5542912/76‑1 (MQ=255)
tGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACt                                                                        >  2:2681505/1‑76 (MQ=255)
 ggCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTg                                                                       <  1:6977990/76‑1 (MQ=255)
 ggCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTg                                                                       >  2:4633194/1‑76 (MQ=255)
 ggCTGACGTGGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTg                                                                       <  1:3769632/76‑1 (MQ=255)
   cTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTgca                                                                     >  1:2374360/1‑74 (MQ=255)
     gACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAaca                                                                   >  1:5304051/1‑76 (MQ=255)
     gACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAaca                                                                   <  2:5863092/76‑1 (MQ=255)
        gTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCg                                                                >  1:6481847/1‑76 (MQ=255)
         ttGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCgg                                                               >  1:498740/1‑76 (MQ=255)
         ttGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCgg                                                               >  1:146608/1‑76 (MQ=255)
         ttGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCgg                                                               <  2:3445045/76‑1 (MQ=255)
         ttGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCgg                                                               <  1:2656706/76‑1 (MQ=255)
             gAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAac                                                           >  1:1071124/1‑76 (MQ=255)
             gAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAac                                                           <  1:2245098/76‑1 (MQ=255)
              aaaTCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAaca                                                          >  1:4322705/1‑76 (MQ=255)
              aaaTCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAaca                                                          >  1:6689015/1‑76 (MQ=255)
               aaTCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAacac                                                         >  1:126089/1‑76 (MQ=255)
                aTCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAAcacc                                                        <  1:2018409/76‑1 (MQ=255)
                aTCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAAcacc                                                        <  2:5049363/76‑1 (MQ=255)
                aTCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAAcacc                                                        >  1:4631632/1‑76 (MQ=255)
                  ccAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAAcaccac                                                      <  2:3353486/76‑1 (MQ=255)
                   cAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACaccacc                                                     >  2:1148237/1‑76 (MQ=255)
                   cAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACaccacc                                                     <  1:1715752/76‑1 (MQ=255)
                    aaaGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCt                                                    >  2:3061685/1‑76 (MQ=255)
                    aaaGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCt                                                    <  2:3737912/76‑1 (MQ=255)
                     aaGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTa                                                   >  2:2739845/1‑76 (MQ=255)
                     aaGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTa                                                   <  2:2240537/76‑1 (MQ=255)
                      aGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTAc                                                  >  1:11610/1‑76 (MQ=255)
                         gcCGAACTCTCCGGTTGCAGGTCGCTCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACaaa                                               <  1:7099120/76‑1 (MQ=255)
                         gcCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACaaa                                               >  1:5975871/1‑76 (MQ=255)
                         gcCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACaaa                                               >  1:2274080/1‑76 (MQ=255)
                          ccGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAg                                              <  2:6758457/76‑1 (MQ=255)
                          ccGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAg                                              >  1:1130470/1‑76 (MQ=255)
                          ccGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAg                                              >  1:4713420/1‑76 (MQ=255)
                             aaCTCTCCGGTTGCAGGTCGCGCTACCGTGTTCGTCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCgg                                           >  2:5139390/1‑76 (MQ=255)
                             aaCTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCgg                                           <  2:4969363/76‑1 (MQ=255)
                              aCTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGt                                          >  2:3764309/1‑76 (MQ=255)
                              aCTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGt                                          >  2:5338729/1‑76 (MQ=255)
                               ctctCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTa                                         <  2:691240/76‑1 (MQ=255)
                               ctctCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTa                                         >  1:3686381/1‑76 (MQ=255)
                                tctcCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTAc                                        >  1:5936578/1‑76 (MQ=255)
                                tctcCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTAc                                        >  1:1914393/1‑76 (MQ=255)
                                tctcCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTAc                                        >  1:507654/1‑76 (MQ=255)
                                tctcCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTAc                                        <  2:3259773/76‑1 (MQ=255)
                                tctcCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTAc                                        >  2:6271434/1‑76 (MQ=255)
                                 ctcCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACa                                       >  1:6718510/1‑76 (MQ=255)
                                 ctcCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACa                                       <  2:7031663/76‑1 (MQ=255)
                                 ctcCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACAGGTAACACCACCTACAAAGCGGTACa                                       >  1:5510067/1‑76 (MQ=255)
                                   ccGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGc                                     <  1:184610/76‑1 (MQ=255)
                                   ccGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGc                                     >  2:2443722/1‑76 (MQ=255)
                                   ccGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGc                                     >  2:2121654/1‑76 (MQ=255)
                                    cGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCg                                    <  2:5646271/76‑1 (MQ=255)
                                    cGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCg                                    >  1:683726/1‑76 (MQ=255)
                                     ggTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGt                                   >  1:402664/1‑76 (MQ=255)
                                     ggTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGt                                   >  2:777742/1‑76 (MQ=255)
                                       ttGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTc                                 <  2:1953915/76‑1 (MQ=255)
                                       ttGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTc                                 >  2:3480239/1‑76 (MQ=255)
                                       ttGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTc                                 <  2:2843526/76‑1 (MQ=255)
                                       ttGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGCACAGCGTTc                                 >  1:3626038/1‑76 (MQ=255)
                                        tGCAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCt                                >  1:1766065/1‑76 (MQ=255)
                                          cAGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGc                              >  2:3647675/1‑76 (MQ=255)
                                           aGGTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGcc                             >  2:852351/1‑76 (MQ=255)
                                            ggTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCg                            <  2:5605169/76‑1 (MQ=255)
                                            ggTCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCg                            >  2:1571270/1‑76 (MQ=255)
                                              tCGCGCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGAc                          <  1:4474929/76‑1 (MQ=255)
                                               cgcgCTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGAcc                         >  2:2500109/1‑76 (MQ=255)
                                                gcgcTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCt                        <  2:289110/76‑1 (MQ=255)
                                                gcgcTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCt                        <  1:4572814/76‑1 (MQ=255)
                                                  gcTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGa                      >  1:1376922/1‑76 (MQ=255)
                                                  gcTACCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGa                      <  2:495110/76‑1 (MQ=255)
                                                    tACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGAtc                    >  2:657679/1‑76 (MQ=255)
                                                     aCCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGAtct                   >  2:5171819/1‑76 (MQ=255)
                                                     aCCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGAtct                   >  2:6400600/1‑76 (MQ=255)
                                                     aCCGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGAtct                   >  1:2806910/1‑76 (MQ=255)
                                                      ccGTGTTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGAtctc                  >  2:738298/1‑76 (MQ=255)
                                                        gtgtTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTTCCGACCTGATCTCCa                <  2:321525/76‑1 (MQ=255)
                                                        gtgtTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCa                >  1:1681528/1‑76 (MQ=255)
                                                        gtgtTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCa                >  1:1547925/1‑76 (MQ=255)
                                                        gtgtTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCa                >  1:3268657/1‑76 (MQ=255)
                                                         tgtTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCAt               >  1:497796/1‑76 (MQ=255)
                                                         tgtTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCAt               >  2:1894767/1‑76 (MQ=255)
                                                          gtTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATc              <  2:4589541/76‑1 (MQ=255)
                                                          gtTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATc              >  1:6900018/1‑76 (MQ=255)
                                                          gtTCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATc              >  1:437463/1‑76 (MQ=255)
                                                            tCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGTTCTCCATCgg            <  2:4876635/76‑1 (MQ=255)
                                                            tCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCgg            >  1:3144938/1‑76 (MQ=255)
                                                            tCATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCgg            >  1:4528465/1‑76 (MQ=255)
                                                             cATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCggg           >  2:4106064/1‑76 (MQ=255)
                                                             cATCTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCggg           >  2:4496388/1‑76 (MQ=255)
                                                                cTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCGGGCCg        <  2:6845260/76‑1 (MQ=255)
                                                                cTTCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCGGGCCg        <  1:5772735/76‑1 (MQ=255)
                                                                  tCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCGGGCCGAt      >  1:3387129/1‑76 (MQ=255)
                                                                  tCCCGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCGGGCCGAt      >  1:5232793/1‑76 (MQ=255)
                                                                   cccGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCGGGCCGAtg     >  2:3901152/1‑76 (MQ=255)
                                                                    ccGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCGGGCCGAtgc    >  2:1440220/1‑76 (MQ=255)
                                                                     cGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCGGGCCGAtgct   >  2:3107274/1‑76 (MQ=255)
                                                                     cGGAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCGGGCCGAtgct   <  2:708870/76‑1 (MQ=255)
                                                                      ggAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCGGGCCGAtgctg  <  2:6941797/76‑1 (MQ=255)
                                                                      ggAACTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCGGGCCGAtgctg  >  1:5195044/1‑76 (MQ=255)
                                                                         |                                                                        
TGGCTGACGTTGCGAAATCCAAAGCGCCGAACTCTCCGGTTGCAGGTCGCGCTACCGTGTTCATCTTCCCGGATCTGAACACCGGTAACACCACCTACAAAGCGGTACAGCGTTCTGCCGACCTGATCTCCATCGGGCCGATGCTG  >  REL606/2360233‑2360378

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