Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,499,829 G→T F325L (TTC→TTA glgC ← glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,499,8290GT97.8% 164.5 / ‑2.4 45F325L (TTC→TTAglgCglucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (25/19);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.90e-01

AACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACATATCCAGCTCCGG  >  REL606/3499756‑3499900
                                                                         |                                                                       
acccGAAGATCACACACCCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTaa                                                                       <  2:605700/74‑1 (MQ=255)
aaCCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGACCCTGCACTaa                                                                       >  2:1636156/1‑76 (MQ=255)
  ccGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAAtt                                                                     >  1:1256437/1‑76 (MQ=255)
     gagaTCACACAACCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCg                                                                  >  2:1585763/1‑76 (MQ=255)
       gaTCACACAACCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCt                                                                >  1:1605438/1‑76 (MQ=255)
              caACCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTa                                                         >  2:1633477/1‑76 (MQ=255)
              caACCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTa                                                         >  2:1017965/1‑76 (MQ=255)
                    ccgGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGAtt                                                   >  2:209598/1‑76 (MQ=255)
                     cgGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTc                                                  >  2:1641993/1‑76 (MQ=255)
                             aGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAgg                                          <  1:1492953/76‑1 (MQ=255)
                              gTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTTATTGTAggg                                         >  1:1015038/1‑75 (MQ=255)
                               tGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTg                                        >  1:43781/1‑76 (MQ=255)
                               tGAGTTAAGGCTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTg                                        <  2:1390945/76‑1 (MQ=255)
                                    tAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAAt                                   <  2:1788850/76‑1 (MQ=255)
                                    tAAGGGTCATCCCGCGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAAt                                   <  2:1883453/76‑1 (MQ=255)
                                     aaGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTTATTGTAGGTGCGAAtt                                  >  1:1707113/1‑76 (MQ=255)
                                     aaGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAAtt                                  >  2:430459/1‑76 (MQ=255)
                                      aGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTg                                 >  2:957680/1‑76 (MQ=255)
                                      aGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTg                                 <  1:522913/76‑1 (MQ=255)
                                      aGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTg                                 >  1:513145/1‑76 (MQ=255)
                                        ggTCATCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGCAATGATTCATTGTAGGTGGGAATTGGc                               >  2:859888/1‑76 (MQ=255)
                                             tCCCGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAAtt                          <  2:1490370/76‑1 (MQ=255)
                                              cccGTGGCTACCGGAGCGATCCTGCCCTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTg                         <  2:255027/76‑1 (MQ=255)
                                               ccGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGc                        >  1:1643224/1‑76 (MQ=255)
                                                cGTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTTATTGTAGGGGCGCATTGGCCAATTGCg                       >  2:465189/1‑76 (MQ=255)
                                                 gTGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGGAATGATTCATTGTAGGTGCGAAATGGCCAATTGCGa                      >  2:274170/1‑76 (MQ=255)
                                                  tGGCTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGAt                     >  2:1737912/1‑76 (MQ=255)
                                                  tGGCTACCGGAGCGAACCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGAt                     <  1:948892/76‑1 (MQ=255)
                                                     cTACCGGAGCTATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGt                  <  2:1375587/76‑1 (MQ=255)
                                                     cTACCGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGt                  <  1:1107214/76‑1 (MQ=255)
                                                      tgccGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTa                 >  2:826767/3‑76 (MQ=255)
                                                        ccGGAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACa               >  2:1279723/1‑76 (MQ=255)
                                                         cGGAGCGATCGTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACat              <  2:1393514/76‑1 (MQ=255)
                                                          ggAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACata             <  1:537225/76‑1 (MQ=255)
                                                           gAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACatat            <  1:798665/76‑1 (MQ=255)
                                                           gAGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACatat            <  1:1655183/76‑1 (MQ=255)
                                                            aGCGATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACATATc           <  1:1966503/76‑1 (MQ=255)
                                                               gATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACATATCCAg        <  2:1000138/76‑1 (MQ=255)
                                                               gATCCTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACATATCCAg        >  2:231808/1‑76 (MQ=255)
                                                                  ccTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACATATCCAGCTc     <  1:31955/76‑1 (MQ=255)
                                                                  ccTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACATATCCAGCTc     >  1:1351120/1‑76 (MQ=255)
                                                                   cTGCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACATATCCAGCTcc    >  2:246661/1‑76 (MQ=255)
                                                                    tGCACTAATTCCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACATATCCAGCTCCg   <  1:845122/76‑1 (MQ=255)
                                                                     gCACTAATTTCGCTGGCGGTAATGGTTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACATATCCAGCTCCgg  <  2:1711802/76‑1 (MQ=255)
                                                                     gCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACATATCCAGCTCCgg  >  2:1170350/1‑76 (MQ=255)
                                                                     gCACTAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACATATCCAGCTCCgg  >  2:338235/1‑76 (MQ=255)
                                                                     gCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACATATCCAGCTCCgg  >  2:1218569/1‑76 (MQ=255)
                                                                         |                                                                       
AACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTGGCCAATTGCGATCGTACATATCCAGCTCCGG  >  REL606/3499756‑3499900

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