| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | NC_000913 | 543,091 | G→C | 18.5% | intergenic (+58/+170) | glxK → / ← allE | glycerate kinase II/S‑ureidoglycine aminohydrolase | 
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 543,091 | 0 | G | C | 18.5% | 168.6 / 34.1 | 65 | intergenic (+58/+170) | glxK/allE | glycerate kinase II/S‑ureidoglycine aminohydrolase | 
| Reads supporting (aligned to +/- strand): ref base G (26/27); new base C (5/7); total (31/34) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.54e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01 | |||||||||||
GCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAG  >  NC_000913/542949‑543232                                                                                                                                              |                                                                                                                                             gCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCgg                                                                                                                                         <  4:1864038/149‑1 (MQ=255) cTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGAGCAGTACTCATGCCGCATCCGGc                                                                                                                                        >  4:1537795/1‑149 (MQ=255)     ccAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGa                                                                                                                                    <  2:1500098/149‑1 (MQ=255)     ccAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGa                                                                                                                                    <  3:1537795/149‑1 (MQ=255)     ccAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGa                                                                                                                                    <  1:38045/149‑1 (MQ=255)     ccAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGa                                                                                                                                    <  1:307318/149‑1 (MQ=255)        gCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGAg                                                                                                                                   >  2:239886/1‑147 (MQ=255)        gCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGAg                                                                                                                                   <  1:239886/147‑1 (MQ=255)                  cAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCa                                                                                                                                       >  1:1047353/1‑133 (MQ=255)                  cAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCa                                                                                                                                       <  2:1047353/133‑1 (MQ=255)                    aTCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGa                                                                                                                                    <  4:1518337/134‑1 (MQ=255)                    aTCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGa                                                                                                                                    >  3:1518337/1‑134 (MQ=255)                     tctcTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCAtc                                                                                                                                            >  4:1603523/1‑124 (MQ=255)                     tctcTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCAtc                                                                                                                                            <  3:1603523/125‑2 (MQ=255)                               aGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATcgcg                                                                                                          >  1:1710031/1‑149 (MQ=255)                               aGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATcgcg                                                                                                          >  3:1869499/1‑149 (MQ=255)                                        aaTATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCta                                                                                                 >  2:1177220/1‑149 (MQ=255)                                                   cgcCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCaga                                                                                      <  1:1177220/149‑1 (MQ=255)                                                      cATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCAt                                                                                                                                      >  1:1720308/1‑98 (MQ=255)                                                      cATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCAt                                                                                                                                      <  2:1720308/98‑1 (MQ=255)                                                      cATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGAt                                                                                   >  2:1552239/1‑149 (MQ=255)                                                           aaaTAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtc                                                                                                                          <  1:570605/105‑1 (MQ=255)                                                           aaaTAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtc                                                                                                                          >  2:570605/1‑105 (MQ=255)                                                               aGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGaaaa                                                                                                               <  1:1472166/112‑1 (MQ=255)                                                               aGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGaaaa                                                                                                               >  2:1472166/1‑112 (MQ=255)                                                                    aGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGAt                                                                     <  1:1552239/149‑1 (MQ=255)                                                                     gggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGc                                                                                                                                        <  4:1232816/81‑1 (MQ=255)                                                                     gggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATc                                                                                                             >  1:1038426/1‑108 (MQ=255)                                                                     gggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATc                                                                                                             <  2:1038426/108‑1 (MQ=255)                                                                     gggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATATATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGc                                                                                                                                        >  3:1232816/1‑81 (MQ=255)                                                                        aaTTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagc                                                                 >  1:1081147/1‑149 (MQ=255)                                                                             aaaaCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgaagc                                                            <  2:1081147/149‑1 (MQ=255)                                                                                  tAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTg                                                                                                                   >  3:235635/1‑89 (MQ=255)                                                                                  tAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTg                                                                                                                   <  4:235635/89‑1 (MQ=255)                                                                                      ccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCAt                                                                                                                                      <  2:441649/66‑1 (MQ=255)                                                                                      ccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCAt                                                                                                                                      >  1:441649/1‑66 (MQ=255)                                                                                      ccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTg                                                                                          >  3:640762/1‑110 (MQ=255)                                                                                      ccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTg                                                                                          <  4:640762/110‑1 (MQ=255)                                                                                      ccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgaa                                                              <  4:1352639/138‑1 (MQ=255)                                                                                      ccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgaa                                                              >  3:1352639/1‑138 (MQ=255)                                                                                        ttACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCaa                                                 >  1:310491/1‑149 (MQ=255)                                                                                           ccTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCa                                                                                        >  1:1338873/1‑107 (MQ=255)                                                                                           ccTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCa                                                                                        <  2:1338873/107‑1 (MQ=255)                                                                                             tttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATcgc                                                                                                           >  4:1369974/1‑86 (MQ=255)                                                                                             tttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATcgc                                                                                                           <  3:1369974/86‑1 (MQ=255)                                                                                             tttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTa                                            <  2:1710031/149‑1 (MQ=255)                                                                                             tttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTa                                            <  4:1869499/149‑1 (MQ=255)                                                                                                   gCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGAt                                                                     >  2:1225846/1‑118 (MQ=255)                                                                                                   gCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGAt                                                                     <  1:1225846/118‑1 (MQ=255)                                                                                                     ggATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaaa                                    <  2:310491/149‑1 (MQ=255)                                                                                                           gATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCa                                                                                 <  1:559609/98‑1 (MQ=255)                                                                                                           gATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCa                                                                                 >  2:559609/1‑98 (MQ=255)                                                                                                           gATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGa                                                                                                                                    <  4:1829890/47‑1 (MQ=21)                                                                                                           gATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGa                                                                                                                                    >  3:1829890/1‑47 (MQ=21)                                                                                                                       aTAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCta                                                                                                 >  1:1753924/1‑70 (MQ=255)                                                                                                                       aTAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCta                                                                                                 <  2:1753924/70‑1 (MQ=255)                                                                                                                                      tCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATCAGCGAAGCGCATCAg                                                     <  4:1400792/99‑1 (MQ=255)                                                                                                                                      tCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAg                                                     >  3:1400792/1‑99 (MQ=255)                                                                                                                                      tCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCaa                                                 <  2:1502492/103‑1 (MQ=255)                                                                                                                                      tCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCaa                                                 >  1:1502492/1‑103 (MQ=255)                                                                                                                                      tCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGtat                          <  2:723804/126‑1 (MQ=255)                                                                                                                                      tCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGtat                          >  1:723804/1‑126 (MQ=255)                                                                                                                                       cATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAg  >  4:1388344/1‑149 (MQ=255)                                                                                                                                        aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTatatat                                                                                            <  3:192492/58‑1 (MQ=255)                                                                                                                                        aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTatatat                                                                                            >  4:192492/1‑58 (MQ=255)                                                                                                                                             ggATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCaaaa            >  4:1533372/1‑133 (MQ=255)                                                                                                                                             ggATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCaaaa            <  3:1533372/133‑1 (MQ=255)                                                                                                                                              |                                                                                                                                             GCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAG  >  NC_000913/542949‑543232 | 
| Alignment Legend | 
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG | 
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑ |