Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,395 T→C S354P (TCC→CCC)  thrA → bifunctional aspartokinase I/homeserine dehydrogenase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,3950TC99.5% 927.2 / ‑4.4 219S354P (TCC→CCC) thrAbifunctional aspartokinase I/homeserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (109/109);  total (109/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 230 total aligned reads displayed.

GGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGT  >  REL606/1323‑1467
                                                                        |                                                                        
ggCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCg                                                                       <  1:5463975/76‑1 (MQ=255)
ggCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCg                                                                       <  1:2546027/76‑1 (MQ=255)
  cATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGaa                                                                     <  2:3195836/76‑1 (MQ=255)
   aTGGCGGCGCGCGTGTTTGCAGCGATGTCACGCTCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAAt                                                                    <  2:332870/76‑1 (MQ=255)
     ggcggcGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATAc                                                                  >  1:3699849/1‑76 (MQ=255)
         gcgcgcGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTa                                                              >  2:6740969/1‑76 (MQ=255)
          cgcgcgTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTAt                                                             >  2:3934676/1‑76 (MQ=255)
             gcgTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAg                                                          <  2:6484757/76‑1 (MQ=255)
              cgTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGt                                                         <  2:1714359/76‑1 (MQ=255)
              cgTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGt                                                         <  1:3412603/76‑1 (MQ=255)
               gtgtTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGtt                                                        >  2:189301/1‑76 (MQ=255)
                  tttGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGGTTACGCAATCATCTCCCGAATACAGTATCAGTTTCt                                                     >  2:5971540/1‑76 (MQ=255)
                  tttGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCt                                                     <  2:5243581/76‑1 (MQ=255)
                         cGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTcc                                              >  1:65835/1‑76 (MQ=255)
                          gATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCg                                             >  2:7622628/1‑76 (MQ=255)
                          gATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCg                                             >  2:1603673/1‑76 (MQ=255)
                            tGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCa                                           >  2:2062678/1‑76 (MQ=255)
                            tGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCa                                           >  1:5636189/1‑76 (MQ=255)
                            tGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCa                                           >  2:2892084/1‑76 (MQ=255)
                              tCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGCATACAGTATCAGTTTCTGCGTTCCGCaac                                         >  1:3455773/1‑75 (MQ=255)
                              tCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCaaa                                         <  2:4185407/76‑1 (MQ=255)
                              tCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCaaa                                         <  2:4432657/76‑1 (MQ=255)
                              tCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCaaa                                         <  2:5065696/76‑1 (MQ=255)
                               cACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAg                                        >  2:722591/1‑76 (MQ=255)
                                aCGCTCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGc                                       <  1:6505705/76‑1 (MQ=255)
                                aCGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGc                                       <  1:7509298/76‑1 (MQ=255)
                                    gcCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACt                                   >  2:152916/1‑76 (MQ=255)
                                    gcCCGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAACCGACt                                   >  2:1633021/1‑76 (MQ=255)
                                     cccGTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACtg                                  <  1:1248771/76‑1 (MQ=255)
                                        gTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACtgtgt                               >  2:6596661/1‑76 (MQ=255)
                                        gTATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACtgtgt                               >  1:6043701/1‑76 (MQ=255)
                                         tATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACtgtgtg                              >  2:530212/1‑76 (MQ=255)
                                         tATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACtgtgtg                              >  1:3504044/1‑76 (MQ=255)
                                         tATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACtgtgtg                              >  2:4789081/1‑76 (MQ=255)
                                         tATTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACtgtgtg                              >  2:3819533/1‑76 (MQ=255)
                                          aTTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGc                             >  1:4528922/1‑76 (MQ=255)
                                          aTTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGc                             >  1:7708857/1‑76 (MQ=255)
                                          aTTTCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGc                             <  1:4159380/76‑1 (MQ=255)
                                           tttCCGTGGTGCTGATTACGCAATCGTCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCg                            <  1:7820043/76‑1 (MQ=255)
                                           tttCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCg                            <  2:397807/76‑1 (MQ=255)
                                           tttCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCg                            <  1:6794013/76‑1 (MQ=255)
                                             tCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAg                          >  2:6655974/1‑76 (MQ=255)
                                             tCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAg                          <  2:5398628/76‑1 (MQ=255)
                                             tCCGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAg                          >  1:4761526/1‑76 (MQ=255)
                                              ccGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGc                         >  2:5724022/1‑76 (MQ=255)
                                              ccGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGc                         >  1:783592/1‑76 (MQ=255)
                                              ccGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGc                         >  1:1971531/1‑76 (MQ=255)
                                               cGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCt                        <  2:576754/76‑1 (MQ=255)
                                               cGTGGTGCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCt                        >  1:7264291/1‑76 (MQ=255)
                                                gtggtgCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTg                       <  1:1049502/76‑1 (MQ=255)
                                                gtggtgCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTg                       >  2:4697083/1‑76 (MQ=255)
                                                gtggtgCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTg                       <  1:1466464/76‑1 (MQ=255)
                                                 tggtgCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGa                      <  1:5726175/76‑1 (MQ=255)
                                                  ggtgCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGaa                     <  1:7121668/76‑1 (MQ=255)
                                                  ggtgCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGaa                     <  1:3335741/76‑1 (MQ=255)
                                                   gtgCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAAc                    >  1:7684854/1‑76 (MQ=255)
                                                   gtgCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAAc                    <  1:6963333/76‑1 (MQ=255)
                                                   gtgCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAAc                    <  1:2968610/76‑1 (MQ=255)
                                                    tgCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACg                   >  1:4404013/1‑76 (MQ=255)
                                                     gCTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACgg                  >  1:6601583/1‑76 (MQ=255)
                                                      cTGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACggg                 <  1:1700521/76‑1 (MQ=255)
                                                       tGATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGc                >  1:5299679/1‑76 (MQ=255)
                                                        gATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCa               >  1:4887743/1‑76 (MQ=255)
                                                        gATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCa               <  1:4320790/76‑1 (MQ=255)
                                                        gATTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCa               >  1:1987873/1‑76 (MQ=255)
                                                        gATTACGCAATCATCTCCAGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCa               >  2:7672274/1‑76 (MQ=255)
                                                         aTTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCaa              >  1:220740/1‑76 (MQ=255)
                                                         aTTACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGCGTGCGAGCTGAACGGGCaa              >  2:6461183/1‑76 (MQ=255)
                                                          ttACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAAt             >  1:2637554/1‑76 (MQ=255)
                                                           tACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATg            >  1:2815201/1‑76 (MQ=255)
                                                           tACGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATg            >  1:242945/1‑76 (MQ=255)
                                                            aCGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGc           <  2:352155/76‑1 (MQ=255)
                                                            aCGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGc           <  2:273784/76‑1 (MQ=255)
                                                            aCGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGc           <  1:6740969/76‑1 (MQ=255)
                                                            aCGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGc           >  2:6922464/1‑76 (MQ=255)
                                                            aCGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGc           <  1:5360647/76‑1 (MQ=255)
                                                            aCGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGc           <  1:2084106/76‑1 (MQ=255)
                                                             cGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCa          >  2:6112296/1‑76 (MQ=255)
                                                             cGCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCa          >  1:1330024/1‑76 (MQ=255)
                                                              gCAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGCCTGTGTGCGAGCTGAACGGGCAATGCAg         <  2:6039863/76‑1 (MQ=255)
                                                               cAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAgg        >  2:5187902/1‑76 (MQ=255)
                                                               cAATCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAgg        <  1:3107993/76‑1 (MQ=255)
                                                                aaTCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGa       <  2:4065218/76‑1 (MQ=255)
                                                                aaTCATCTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGAAATGCAGGa       <  2:71508/76‑1 (MQ=255)
                                                                 atcatcTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGaa      <  1:6262110/76‑1 (MQ=255)
                                                                 atcatcTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGaa      >  1:1180090/1‑76 (MQ=255)
                                                                 atcatcTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGaa      <  2:3070734/76‑1 (MQ=255)
                                                                  tcatcTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAag     <  2:2951933/76‑1 (MQ=255)
                                                                  tcatcTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAag     <  1:352014/76‑1 (MQ=255)
                                                                  tcatcTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAag     <  1:4092383/76‑1 (MQ=255)
                                                                   catcTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAaga    >  1:3274063/1‑76 (MQ=255)
                                                                   catcTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAaga    <  1:5876711/76‑1 (MQ=255)
                                                                   catcTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAaga    >  2:4541623/1‑76 (MQ=255)
                                                                   catcTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAaga    >  1:1589333/1‑76 (MQ=255)
                                                                   catcTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAaga    >  2:1212234/1‑76 (MQ=255)
                                                                   catcTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAaga    >  2:1424154/1‑76 (MQ=255)
                                                                    atcTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAagag   <  1:4941993/76‑1 (MQ=255)
                                                                    atcTCCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAagag   >  2:499088/1‑76 (MQ=255)
                                                                     tctcCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGt  <  1:2790244/76‑1 (MQ=255)
                                                                     tctcCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGt  >  2:5691167/1‑76 (MQ=255)
                                                                        |                                                                        
GGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGT  >  REL606/1323‑1467

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