Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,345,860 T→G N397H (AAC→CAC)  aspA ← aspartate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,345,8600TG99.1% 1406.7 / ‑0.3 321N397H (AAC→CAC) aspAaspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (132/186);  total (133/188)
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.95e-01

Only 100 of 328 total aligned reads displayed.

TGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAA  >  REL606/4345787‑4345932
                                                                         |                                                                        
tgACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAg                                                                        >  2:6307914/1‑76 (MQ=255)
tgACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAg                                                                        >  1:6137840/1‑76 (MQ=255)
   ccGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAg                                                                     <  1:1733883/76‑1 (MQ=255)
    cGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGt                                                                    <  2:395327/76‑1 (MQ=255)
    cGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGt                                                                    >  1:2525690/1‑76 (MQ=255)
    cGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAACAGt                                                                    <  1:4724658/76‑1 (MQ=255)
     gATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTg                                                                   <  1:6952567/76‑1 (MQ=255)
      aTGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGCAACCTTCGCACACTTCTTTGTGAGCAGTGa                                                                  >  1:1100036/1‑76 (MQ=255)
       tGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGAt                                                                 >  1:882368/1‑76 (MQ=255)
       tGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGAt                                                                 >  1:1461969/1‑76 (MQ=255)
        gAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATg                                                                <  2:7520794/76‑1 (MQ=255)
        gAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATg                                                                >  2:578062/1‑76 (MQ=255)
        gAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATg                                                                >  1:3100026/1‑76 (MQ=255)
         aaCGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGc                                                               <  2:5669446/76‑1 (MQ=255)
           cGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCg                                                             <  1:7240108/76‑1 (MQ=255)
            gggTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGt                                                            >  1:6456699/1‑76 (MQ=255)
            gggTGCAGGTAAGTAACGATACCGTTAGAGTTGTAAACGTAACGTCCGCACACTTCTTTGTGCGCAGTGATGCCGt                                                            >  2:7461595/1‑76 (MQ=255)
             ggTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGtt                                                           <  2:135738/76‑1 (MQ=255)
             ggTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGtt                                                           <  1:6389097/76‑1 (MQ=255)
             ggTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGtt                                                           <  1:3691178/76‑1 (MQ=255)
               ttCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTaa                                                         <  1:7375575/76‑1 (MQ=255)
                tCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAAt                                                        <  1:1973414/76‑1 (MQ=255)
                 cAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATg                                                       <  1:2427046/76‑1 (MQ=255)
                 cAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATg                                                       <  1:183241/76‑1 (MQ=255)
                 cAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATg                                                       <  1:6839105/76‑1 (MQ=255)
                  aGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGc                                                      <  1:6240200/76‑1 (MQ=255)
                   ggTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCa                                                     <  1:5710591/76‑1 (MQ=255)
                   ggTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCa                                                     <  1:5444882/76‑1 (MQ=255)
                     taagtaaCGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCAtt                                                   <  2:2852735/76‑1 (MQ=255)
                     taagtaaCGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCAtt                                                   <  1:5370825/76‑1 (MQ=255)
                     taagtaaCGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCAtt                                                   <  1:4229723/76‑1 (MQ=255)
                      aagtaaCGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCAttt                                                  >  1:481212/1‑76 (MQ=255)
                        gtaaCGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCAttttt                                                <  1:4438405/76‑1 (MQ=255)
                          aaCGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTcc                                              <  1:7047666/76‑1 (MQ=255)
                          aaCGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTcc                                              <  1:1280102/76‑1 (MQ=255)
                          aaCGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTcc                                              <  2:3568988/76‑1 (MQ=255)
                          aaCGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTcc                                              >  1:1596406/1‑76 (MQ=255)
                          aaCGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTcc                                              >  1:2904217/1‑76 (MQ=255)
                          aaCGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTcc                                              <  1:6013269/76‑1 (MQ=255)
                          aaCGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTcc                                              >  1:6299057/1‑76 (MQ=255)
                          aaCGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTcc                                              >  1:4128/1‑76 (MQ=255)
                            cGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCcag                                            >  1:7066156/1‑76 (MQ=255)
                             gATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCcagc                                           >  1:525467/1‑76 (MQ=255)
                             gATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCcagc                                           >  1:4976080/1‑76 (MQ=255)
                               tACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCcagcag                                         <  1:2218872/76‑1 (MQ=255)
                                aCCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAgg                                        <  2:6121911/76‑1 (MQ=255)
                                aCCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAgg                                        >  2:2636491/1‑76 (MQ=255)
                                 ccGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGt                                       >  2:5094459/1‑76 (MQ=255)
                                 ccGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGt                                       >  2:7393138/1‑76 (MQ=255)
                                 ccGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGt                                       <  1:6799238/76‑1 (MQ=255)
                                   gATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTg                                     <  1:2573724/76‑1 (MQ=255)
                                   gATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTg                                     <  2:5348482/76‑1 (MQ=255)
                                   gATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTg                                     <  1:992919/76‑1 (MQ=255)
                                   gATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTg                                     <  2:3137003/76‑1 (MQ=255)
                                   gATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTg                                     <  1:4039721/76‑1 (MQ=255)
                                    aTAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGt                                    >  2:6516429/1‑76 (MQ=255)
                                    aTAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGt                                    <  2:5461602/76‑1 (MQ=255)
                                     tAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTa                                   <  1:712932/76‑1 (MQ=255)
                                     tAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTa                                   <  1:2990203/76‑1 (MQ=255)
                                     tAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTa                                   >  2:7718391/1‑76 (MQ=255)
                                      agagTTGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAg                                  >  2:2789594/1‑76 (MQ=255)
                                          ttGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAg                              <  1:3508028/76‑1 (MQ=255)
                                          ttGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAg                              >  2:2271863/1‑76 (MQ=255)
                                           tGTAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGc                             >  1:5221038/1‑76 (MQ=255)
                                             tAAACGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGt                           >  1:7550731/1‑76 (MQ=255)
                                              aaaCGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGtt                          <  2:4159740/76‑1 (MQ=255)
                                                aCGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTgg                        <  2:7358954/76‑1 (MQ=255)
                                                 cGTAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGt                       >  1:3972366/1‑76 (MQ=255)
                                                   tAACCTTCGCACACTTCTTTGTGAGCAGTGATGCTGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCa                     <  2:6392425/76‑1 (MQ=255)
                                                   tAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCa                     <  1:2348448/76‑1 (MQ=255)
                                                   tAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCa                     <  2:2332137/76‑1 (MQ=255)
                                                   tAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCa                     >  1:6719941/1‑76 (MQ=255)
                                                   tAACCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCa                     >  1:7449046/1‑76 (MQ=255)
                                                    aaCCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAg                    >  2:669400/1‑76 (MQ=255)
                                                    aaCCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAg                    >  1:2171593/1‑76 (MQ=255)
                                                    aaCCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAg                    >  1:6368183/1‑76 (MQ=255)
                                                    aaCCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAg                    <  1:5145043/76‑1 (MQ=255)
                                                     aCCTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGa                   >  2:1891764/1‑76 (MQ=255)
                                                      ccTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGaa                  <  2:4898481/76‑1 (MQ=255)
                                                      ccTTCGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGaa                  >  1:1683722/1‑76 (MQ=255)
                                                       cTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAAt                 <  2:3352741/76‑1 (MQ=255)
                                                          cGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAAtgtg              >  2:1740802/1‑76 (MQ=255)
                                                          cGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAAtgtg              <  2:1622576/76‑1 (MQ=255)
                                                          cGCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAAtgtg              <  2:7626451/76‑1 (MQ=255)
                                                           gCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGa             <  1:5518997/76‑1 (MQ=255)
                                                           gCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGa             >  2:5996601/1‑76 (MQ=255)
                                                           gCACACTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGa             >  2:7045652/1‑76 (MQ=255)
                                                            cacaCTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGaa            >  2:7169949/1‑76 (MQ=255)
                                                            cacaCTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGaa            >  1:3303487/1‑76 (MQ=255)
                                                            cacaCTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGaa            >  2:459495/1‑76 (MQ=255)
                                                            cacaCTACTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGaa            <  1:6553103/76‑1 (MQ=255)
                                                             acacTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAAc           <  1:1178096/76‑1 (MQ=255)
                                                             acacTTCTTTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAAc           >  2:3479917/1‑76 (MQ=255)
                                                                cttcttTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGa        >  1:2721317/1‑76 (MQ=255)
                                                                  tcttTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGAtt      <  2:4331943/76‑1 (MQ=255)
                                                                  tcttTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGAtt      >  1:1997337/1‑76 (MQ=255)
                                                                   cttTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTc     >  2:6934799/1‑76 (MQ=255)
                                                                   cttTGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTc     <  1:7645067/76‑1 (MQ=255)
                                                                    tttGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCg    >  1:1015586/1‑76 (MQ=255)
                                                                     ttGTGAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGa   >  2:2111412/1‑76 (MQ=255)
                                                                         |                                                                        
TGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAA  >  REL606/4345787‑4345932

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