Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,499,798 T→G T336P (ACC→CCC)  glgC ← glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,499,7980TG98.0% 197.3 / ‑4.9 50T336P (ACC→CCC) glgCglucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (20/29);  total (20/30)
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.93e-01

CGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTG  >  REL606/3499725‑3499869
                                                                         |                                                                       
cgcgAGAACAGACCGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCa                                                                       <  2:1020770/76‑1 (MQ=255)
cgcgAAAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCa                                                                       <  1:1339077/76‑1 (MQ=255)
 gcgAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCAt                                                                      >  1:376373/1‑76 (MQ=255)
 gcgAGAACAGAACGGACTGCACCACCACCGAACCAGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCAt                                                                      >  2:283743/1‑76 (MQ=255)
  cgAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATc                                                                     >  1:875334/1‑76 (MQ=255)
    agaacagaacGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATccc                                                                   >  1:777970/1‑76 (MQ=255)
       acagaacGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTg                                                                <  1:244865/76‑1 (MQ=255)
        cagaacGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTgg                                                               >  2:1209154/1‑76 (MQ=255)
           aacGGACTGCCCCACCCCCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTa                                                            <  1:200752/76‑1 (MQ=255)
           aacGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTa                                                            <  2:459250/76‑1 (MQ=255)
           aacGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTa                                                            <  1:1671177/76‑1 (MQ=255)
              ggACTGCCCCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCg                                                         <  2:803014/76‑1 (MQ=255)
                 cTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAg                                                      <  1:154513/76‑1 (MQ=255)
                 cTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAg                                                      >  1:895468/1‑76 (MQ=255)
                 cTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAg                                                      >  2:620236/1‑76 (MQ=255)
                  tGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGc                                                     >  1:745534/1‑76 (MQ=255)
                   gCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCg                                                    >  1:1894111/1‑76 (MQ=255)
                    caccaccacCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGa                                                   >  1:970105/1‑76 (MQ=255)
                     accaccaccGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGAt                                                  <  2:1249354/76‑1 (MQ=255)
                        cccaccGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCt                                               <  2:1504237/75‑1 (MQ=255)
                        accaccGAATCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGCGGCTACCGGAGCGATccc                                               >  1:1677857/1‑75 (MQ=255)
                        accaccGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCt                                               <  2:1663146/76‑1 (MQ=255)
                           cccGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCa                                            <  1:1415767/75‑1 (MQ=255)
                             cGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACg                                          >  1:1751925/1‑76 (MQ=255)
                                  cGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGGTCCTGCACGAAttt                                     >  2:108341/1‑76 (MQ=255)
                                  cGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAAttt                                     <  2:1296295/76‑1 (MQ=255)
                                   ggAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTc                                    <  2:1340186/76‑1 (MQ=255)
                                    gagaTCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCg                                   <  1:1539907/76‑1 (MQ=255)
                                       aTCACACAACCTCCCGCAACCAGTGAGTGAAAGGGCATCCCGCGGCTATCATTGAGATTCTGCACGAATTTCggcg                                >  2:5668/1‑73 (MQ=255)
                                         cacacaACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGc                              <  2:513073/76‑1 (MQ=255)
                                            acaACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGt                           >  1:1209684/1‑76 (MQ=255)
                                            acaACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGt                           <  1:155122/76‑1 (MQ=255)
                                             caACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTa                          >  1:1523666/1‑76 (MQ=255)
                                             caACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTa                          <  1:1546688/76‑1 (MQ=255)
                                             accccgccgGAAACCAGTGAGTGAAGGGGCACCTCGTCGCTACCGGAGCGATCCTGCGCGAAGTTCCCTGGCGGTa                          <  2:227277/73‑1 (MQ=255)
                                              aaCCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTaa                         >  1:260297/1‑76 (MQ=255)
                                                  gccgGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGAt                     <  2:1417386/76‑1 (MQ=255)
                                                    cgGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGGAATGATTc                   >  2:443019/1‑76 (MQ=255)
                                                      gAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCAt                 >  1:1448402/1‑76 (MQ=255)
                                                       aaaCCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCAtt                >  2:915066/1‑76 (MQ=255)
                                                         aCCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGt              <  1:699442/76‑1 (MQ=255)
                                                         aCCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGt              <  1:699436/76‑1 (MQ=255)
                                                           cAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAg            <  1:1028668/76‑1 (MQ=255)
                                                           cAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAg            <  1:859590/76‑1 (MQ=255)
                                                             gTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGt          <  2:929936/76‑1 (MQ=255)
                                                              tGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTg         <  1:130281/76‑1 (MQ=255)
                                                                 gTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGcca      >  1:1021458/1‑74 (MQ=255)
                                                                  ttAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGaa     <  1:1209154/76‑1 (MQ=255)
                                                                   tAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAAt    <  2:421177/76‑1 (MQ=255)
                                                                   tAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAAt    <  2:971931/76‑1 (MQ=255)
                                                                     aGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTg  <  1:1508944/76‑1 (MQ=255)
                                                                     aGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTg  <  2:590673/76‑1 (MQ=255)
                                                                         |                                                                       
CGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTG  >  REL606/3499725‑3499869

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