Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,500,634 T→G K57T (AAG→ACG)  glgC ← glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,500,6340TG97.0% 250.9 / ‑2.8 66K57T (AAG→ACG) glgCglucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (27/37);  total (28/38)
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.52e-01

GGGTTATCACGCCCATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCACGTCCTCCC  >  REL606/3500563‑3500705
                                                                       |                                                                       
gggTTATCACGCCCATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTAcc                                                                     >  2:979841/1‑76 (MQ=255)
gggTTATCACGCCCATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTAcc                                                                     >  2:492551/1‑76 (MQ=255)
  gTTATCACGCCCATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACTTAccgc                                                                   >  1:1250111/1‑76 (MQ=255)
   ttATCACGCCCATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTAccgcc                                                                  <  1:992004/76‑1 (MQ=255)
     aTCACGCCCATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGa                                                                <  2:1504043/76‑1 (MQ=255)
        aCGCCCATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAgt                                                             <  1:1424061/76‑1 (MQ=255)
         cGCCCATACGACGGATCCCGGAGATGCNGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAgtg                                                            >  2:7828/1‑76 (MQ=255)
          gCCCATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAgtgt                                                           <  2:895653/76‑1 (MQ=255)
          gCCCATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAgtgt                                                           >  1:1602636/1‑76 (MQ=255)
          gCCCATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAgtgt                                                           <  2:1535424/76‑1 (MQ=255)
           cccATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTa                                                          >  1:628081/1‑76 (MQ=255)
            ccATACGACGGATCCCGGAGTTGATGCAGTTAGTCAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTAc                                                         <  1:570838/76‑1 (MQ=255)
            ccATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTAc                                                         <  1:806536/76‑1 (MQ=255)
            ccATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTAc                                                         <  1:1016737/76‑1 (MQ=255)
             cATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACg                                                        <  1:653786/76‑1 (MQ=255)
             cATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGACGAAGTGTACg                                                        <  1:116372/76‑1 (MQ=255)
             cANACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGTAGTGGNCg                                                        <  2:1660364/76‑1 (MQ=255)
              aTACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACgg                                                       <  1:440693/76‑1 (MQ=255)
              aTACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACgg                                                       >  1:684860/1‑76 (MQ=255)
                 cgacgGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCg                                                    <  2:1229362/76‑1 (MQ=255)
                 cgacgGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCg                                                    >  2:1033067/1‑76 (MQ=255)
                 cgacgGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCg                                                    <  2:357205/76‑1 (MQ=255)
                      gATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGtttt                                               >  2:1291161/1‑76 (MQ=255)
                      gATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGtttt                                               <  2:1498410/76‑1 (MQ=255)
                      gATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGtttt                                               <  1:1402552/76‑1 (MQ=255)
                         cccGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATTATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCt                                            <  2:1758094/76‑1 (MQ=255)
                         cccGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCt                                            >  1:1293747/1‑76 (MQ=255)
                         cccGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCt                                            <  2:684859/76‑1 (MQ=255)
                           cGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCg                                          <  2:929370/76‑1 (MQ=255)
                           cGGAGTTGATGCAGCTAGATAGCGCAAAGTCGATAATGCGTAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCg                                          <  2:44107/76‑1 (MQ=255)
                            ggAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGc                                         <  1:886309/76‑1 (MQ=255)
                             gAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCt                                        <  2:607200/76‑1 (MQ=255)
                              aGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCtt                                       <  2:1268540/76‑1 (MQ=255)
                              aGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACTGCCGAAGTGTACGACCGGTTTTTTTCGCtt                                       >  1:36020/1‑76 (MQ=255)
                               gTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTa                                      >  2:1724314/1‑76 (MQ=255)
                                    tGCAGTGAGACAGCGCAAAGTGGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGNTTATTTGt                                 <  1:1960580/76‑1 (MQ=255)
                                     gCAGTTAGACGGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGtt                                <  1:1570051/76‑1 (MQ=255)
                                      cAGTTAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTa                               <  1:484303/76‑1 (MQ=255)
                                          tAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATc                           <  2:1177159/76‑1 (MQ=255)
                                          tAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATc                           <  1:1633837/76‑1 (MQ=255)
                                          tAGACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATc                           >  2:452633/1‑76 (MQ=255)
                                            gACAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCt                         <  1:190158/76‑1 (MQ=255)
                                             aCAGCGCAAAGTCGNTAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATGGGTTAAATCCtt                        <  1:1746000/76‑1 (MQ=255)
                                              cAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTc                       <  1:1492128/76‑1 (MQ=255)
                                              cAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTc                       <  2:812412/76‑1 (MQ=255)
                                              cAGCGCAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTc                       >  2:1282796/1‑76 (MQ=255)
                                                   cAAAGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCg                  >  2:1959003/1‑76 (MQ=255)
                                                    aaaGTCGATCATGCGGAGCGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTNAANTCTTCAGGCgg                 >  2:328146/1‑76 (MQ=255)
                                                    aaaGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCgg                 <  1:595008/76‑1 (MQ=255)
                                                    aaaGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCgg                 >  2:68376/1‑76 (MQ=255)
                                                    aaaGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCgg                 >  2:1868294/1‑76 (MQ=255)
                                                    aaaGTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGCCTTTGCTCGATTATTGGTTAAGTCCNTCAGGCgg                 >  2:348924/1‑76 (MQ=255)
                                                       gTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCgggca              >  2:461711/1‑74 (MQ=255)
                                                       gTCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTa              <  2:657233/76‑1 (MQ=255)
                                                       gTCGATAATGAGGAACGGACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTa              >  1:159635/1‑76 (MQ=255)
                                                        tCGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTAc             >  2:1825275/1‑76 (MQ=255)
                                                         cGATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTAcc            >  1:834503/1‑76 (MQ=255)
                                                          gATAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCa           <  1:973535/76‑1 (MQ=255)
                                                          gATAATGCCGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCa           >  1:1889501/1‑76 (MQ=255)
                                                           aTAATGCGGAATGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCAc          >  1:51065/1‑76 (MQ=255)
                                                           aTAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTACATCCTTCAGGCGGGTACCAc          >  1:2029983/1‑76 (MQ=255)
                                                           aTAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCAc          <  1:488991/76‑1 (MQ=255)
                                                           aTAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCAc          >  2:1204620/1‑76 (MQ=255)
                                                            tAATGCGGAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCACg         <  1:636803/76‑1 (MQ=255)
                                                                   gAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCACGTCCTccc  >  1:1179724/1‑76 (MQ=255)
                                                                   gAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCACGTCCTccc  <  1:671124/76‑1 (MQ=255)
                                                                   gAACGTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCACGTCCTccc  <  2:1178087/76‑1 (MQ=255)
                                                                       |                                                                       
GGGTTATCACGCCCATACGACGGATCCCGGAGTTGATGCAGTTAGACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCACGTCCTCCC  >  REL606/3500563‑3500705

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: