Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,274,418 A→C 9.6% Q125H (CAA→CAC)  glcD → glycolate oxidase FAD‑linked subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,274,4180AC 9.6% 147.7 / 8.1 52Q125H (CAA→CAC) glcDglycolate oxidase FAD‑linked subunit
Reads supporting (aligned to +/- strand):  ref base A (21/26);  new base C (3/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

GCGCGGCGTACCGGTGGTTGCGCGTGGCGCCGGTACCGGCCTGTCTGGTGGTGCCCTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCT  >  NC_002947/4274274‑4274564
                                                                                                                                                |                                                                                                                                                  
gcgcGGCGTACCGGTGGTTGCGCGTGGCGCCGGTACCGGCCTGTCTGGTGGTGCCCTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCgg                                                                                                                                                <  4:225632/149‑1 (MQ=255)
      cgTACCGGTGGTTGCGCGTGGCGCCGGTACCGGCCTGTCTGGTGGTGCCCTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGGTTGCACGCGGGCGTaa                                                                                                                                          <  8:181389/149‑2 (MQ=255)
       gTACCGGTGGTTGCGCGTGGCGCCGGTACCGGCCTGTCTGGTGGTGCCCTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACg                                                                                                                                         <  8:187268/149‑1 (MQ=255)
                 ttGCGCGTGGCGCCGGTACCGGCCTGTCTGGTGGTGCCCTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGcc                                                                                                                               <  1:309820/149‑1 (MQ=255)
                 ttGCGCGTGGCGCCGGTACCGGCCTGTCTGGTGGTGCCCTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGcc                                                                                                                               <  6:51739/149‑1 (MQ=255)
                       gTGGCGCCGGTACCGGCCTGTCTGGTGGTGCCCTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTcc                                                                                                                         <  7:287181/149‑1 (MQ=255)
                          gcgcCGGTACCGGCCTGTCTGGTGGTGCCCTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGGTTGCACGCGGGCGTaa                                                                                                                                          >  2:114999/1‑128 (MQ=255)
                          gcgcCGGTACCGGCCTGTCTGGTGGTGCCCTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGGTTGCACGCGGGCGTaa                                                                                                                                          <  1:114999/129‑2 (MQ=255)
                               ggTACCGGCCTGTCTGGTGGTGCCCTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGcc                                                                                                                               <  6:159184/135‑1 (MQ=255)
                               ggTACCGGCCTGTCTGGTGGTGCCCTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGcc                                                                                                                               >  5:159184/1‑135 (MQ=255)
                                          gTCTGGTGGTGCCCTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATg                                                                                                      >  7:290930/1‑149 (MQ=255)
                                           tCTGGTGGTGCCCTGCCGCTGGCCAAGGGCATCCAGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATg                                                                                                      <  5:69251/148‑1 (MQ=255)
                                               gtggtgCCCTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGt                                                                                                    <  2:318888/146‑1 (MQ=255)
                                                  gtgCCCTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGta                                                                                               >  4:14424/1‑148 (MQ=255)
                                                     cccTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGtactac                                                                                           >  6:20200/1‑149 (MQ=255)
                                                      ccTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCATCACCCGTATCCTCGAGGTCAACCCGCACGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCAATCACCCAGGCCGCCGCACCCCATGGAATGTCCTACg                                                                                          >  6:207091/1‑149 (MQ=255)
                                                             cTGGGCAAGGGCCTCCTGCTGGTGTTGGCGCGCGTGAACCGGATGCGCGAGGTCAACCCGCAGGGGCGTTAAGCCCGCGTGGAACCGGGCGGAAGCAACCGGGCCATCTCCCAGGCCGCCGCCCCCCCGGGGATGTTCTACGCCCCCgg                                                                                   >  2:212780/1‑148 (MQ=255)
                                                              tGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGAc                                                                                  >  2:196909/1‑149 (MQ=255)
                                                                gCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGAccc                                                                                <  8:290930/149‑1 (MQ=255)
                                                                gCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGAcc                                                                                 <  3:14424/148‑1 (MQ=255)
                                                                 ccAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGAATCCTCGAGGTCACCCCGCAGGGGCGTTACGCCCGCGTGCAACAGGGCGTACGCTACCTGGCCAGATCCCAGGCCGCCGCCCCCCATCGAATGTACTCAGCACCCGAccca                                                                               >  6:287544/1‑148 (MQ=255)
                                                                  cAAGGGCATGCTTCGGGTGATGGCGCGGTTGAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGCCGTACGCAACCTg                                                                                                                                  <  7:203830/97‑1 (MQ=255)
                                                                  cAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATg                                                                                                      <  4:255028/125‑1 (MQ=255)
                                                                                   tGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGc                                                             >  5:330901/1‑149 (MQ=255)
                                                                                          gcgcTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCg                                                      <  6:330901/149‑1 (MQ=255)
                                                                                          gcgcTTCAACAGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCAGATCgga                                                        >  5:200604/1‑145 (MQ=255)
                                                                                          gcgcTTCAAACGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGc                                                                 <  6:200604/138‑1 (MQ=255)
                                                                                            gcTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGt                                                    <  1:196909/149‑1 (MQ=255)
                                                                                             cTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGTg                                                   <  5:20200/149‑1 (MQ=255)
                                                                                                   ccGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGTg                                                   >  2:239040/1‑143 (MQ=255)
                                                                                                   ccGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGTg                                                   <  1:239040/143‑1 (MQ=255)
                                                                                                   ccGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGTGGCAACg                                             >  6:312111/1‑149 (MQ=255)
                                                                                                    cGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGTGGCAACGt                                            >  5:159245/1‑149 (MQ=255)
                                                                                                         ccTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATg                                                                                                      >  4:195065/1‑86 (MQ=255)
                                                                                                         ccTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATg                                                                                                      <  3:195065/86‑1 (MQ=255)
                                                                                                          cTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCAt                                                                                                       >  7:151710/1‑84 (MQ=255)
                                                                                                          cTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCAt                                                                                                       <  8:151710/84‑1 (MQ=255)
                                                                                                           tCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTg                                                              >  5:250038/1‑124 (MQ=255)
                                                                                                           tCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTg                                                              <  6:250038/124‑1 (MQ=255)
                                                                                                                      ccGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCCTCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTgg                          >  8:242125/1‑149 (MQ=255)
                                                                                                                           ggggCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCt                                                            >  7:236179/1‑110 (MQ=255)
                                                                                                                           ggggCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCt                                                            <  8:236179/110‑1 (MQ=255)
                                                                                                                              gCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACt                  >  4:313896/1‑149 (MQ=255)
                                                                                                                                gTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCt                                                                           <  5:72382/90‑1 (MQ=255)
                                                                                                                                gTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCt                                                                           >  6:72382/1‑90 (MQ=255)
                                                                                                                                 ttACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCAc                                                                                       <  2:62170/77‑1 (MQ=255)
                                                                                                                                 ttACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCAc                                                                                       >  1:62170/1‑77 (MQ=255)
                                                                                                                                 gtACGCCCGGTTGCACGCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGcc               >  4:252056/2‑149 (MQ=255)
                                                                                                                                 gtACGCCCGGTTGCACGCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGcc               >  8:318312/2‑149 (MQ=255)
                                                                                                                                          cgTGGAACCGGGCGTACGCAACCTGGCCATCTCGCAGGCCGCCGCACCCCATGGTATGTACTAGGGACCCTACCCTTCCTCGCAAATTGCCTGCTCGATGGGTGGCAACGTCGCCTAAAACGCCGGTGGCGTGCACTGCCTCAAGTACg      <  5:207091/149‑1 (MQ=255)
                                                                                                                                          cctGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACg      <  5:312111/147‑1 (MQ=255)
                                                                                                                                              cAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCt  <  2:119167/149‑1 (MQ=255)
                                                                                                                                                aCCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGc                >  2:49018/1‑133 (MQ=255)
                                                                                                                                                aCCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGc                <  1:49018/133‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GCGCGGCGTACCGGTGGTTGCGCGTGGCGCCGGTACCGGCCTGTCTGGTGGTGCCCTGCCGCTGGCCAAGGGCATCCTGCTGGTGATGGCGCGCTTCAACCGTATCCTCGAGGTCAACCCGCAGGGGCGTTACGCCCGCGTGCAACCGGGCGTACGCAACCTGGCCATCTCCCAGGCCGCCGCACCCCATGGTATGTACTACGCACCCGACCCTTCCTCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCT  >  NC_002947/4274274‑4274564

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑