Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,123,585 +G coding (953/2286 nt) glcF ← glycolate oxidase iron‑sulfur subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,123,5821.G93.8% 196.3 / 9.1 80coding (956/2286 nt)glcFglycolate oxidase iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base G (46/29);  total (49/31)
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.98e-01

ACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGA  >  NC_000913/3123507‑3123651
                                                                               |                                                                      
agcGGGGTTAAACAGCCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGc                                                                         <  1:2275171/74‑1 (MQ=255)
aCCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGc                                                                         <  1:1551985/76‑1 (MQ=255)
  cGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCgg                                                                      <  1:1778973/76‑1 (MQ=255)
   ggggTTAAAGCGGCCGCAGGGGTCGGGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCggg                                                                     <  2:497099/76‑1 (MQ=255)
   ggggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCggg                                                                     <  2:4150715/76‑1 (MQ=255)
   ggggTAAACCACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCggg                                                                     <  1:2073012/76‑1 (MQ=255)
    gggTTAGCCACGCCGCAAGGGTCGGGCGGCTGTTTAAGCTGGTGGTGATAGCGGAATAAAGGAGCCGATAGCgggg                                                                    <  2:2347474/76‑2 (MQ=255)
    gggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCgggg                                                                    >  2:2508719/1‑75 (MQ=255)
    gggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCgg                                                                      >  1:5122824/1‑74 (MQ=255)
    gggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCgg                                                                      <  2:5122823/74‑1 (MQ=255)
    gggTGAAACTCGCCGCAAGGGGCGGGCTGCGGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCgggg                                                                    <  2:5390662/76‑2 (MQ=255)
     ggTTAAACACGCGGCAGGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCggggc                                                                   <  2:2825305/76‑3 (MQ=255)
     ggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCggggc                                                                   >  1:2736822/1‑74 (MQ=255)
          aaCCCGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAGGGAGCCGATAGCGGGGCAAAGc                                                              <  2:1108827/76‑1 (MQ=255)
           acacGCCGCAGGGGTCGGGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGcc                                                             <  2:3806072/76‑1 (MQ=255)
           acacGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGcc                                                             >  1:4129203/1‑76 (MQ=255)
           acacGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGcc                                                             >  2:4495112/1‑76 (MQ=255)
            cacGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGccc                                                            >  2:1549088/1‑75 (MQ=255)
             acGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGccac                                                           >  1:86243/1‑76 (MQ=255)
               gccgcAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGccacca                                                         >  1:5399183/1‑76 (MQ=255)
               gccgcAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGACAGCGGGGCCAAGccaccc                                                         >  2:5512904/1‑75 (MQ=255)
                 cgcAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCCtc                                                       >  2:2888849/1‑76 (MQ=255)
                 cgcAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCAtc                                                       >  1:3018425/1‑76 (MQ=255)
                 cgcAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCAtc                                                       <  1:5741049/76‑1 (MQ=255)
                 cgcAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCAtc                                                       >  1:5616235/1‑76 (MQ=255)
                  gcAGGGGTCGAGCGGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGGCACCAtct                                                      <  2:1178292/76‑1 (MQ=255)
                   caggggTGGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCAtctc                                                     <  1:5061755/73‑1 (MQ=255)
                   cAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCCtctc                                                     >  2:1523967/1‑76 (MQ=255)
                    aaGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTcc                                                    >  2:1862968/1‑76 (MQ=255)
                     aGGGTGGAGTTGCTGTTTAGGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCg                                                   <  1:4572845/76‑1 (MQ=255)
                     aGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCAATAGCGGGGCAAAGCCACCATCTCCg                                                   >  1:1933455/1‑76 (MQ=255)
                      gggTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCgg                                                  <  1:5095661/76‑1 (MQ=255)
                         tCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCACAGCCACCATCTCCCGccc                                               >  2:705052/1‑74 (MQ=255)
                         tCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGccc                                               >  2:5807498/1‑74 (MQ=255)
                          cGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCCAAGCCACCATCTCCGGCACt                                              >  1:2961876/1‑76 (MQ=255)
                            aGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTaa                                            >  1:5190075/1‑76 (MQ=255)
                             gctgctGTTTAGGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTaaa                                           <  1:629937/76‑1 (MQ=255)
                              gtgctgTTTAAGGGGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAg                                          <  2:5274966/75‑1 (MQ=255)
                              ctgctgTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAg                                          >  2:4535848/1‑76 (MQ=255)
                                 ctgTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCgg                                       >  1:1766034/1‑76 (MQ=255)
                                 ctgTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGAGGGGCAAAGCCAACATCACCGGCACTAAAGCgg                                       >  2:5359145/1‑76 (MQ=255)
                                  tgTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCggg                                      <  1:2189306/76‑1 (MQ=255)
                                  tgTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCggg                                      >  2:1104788/1‑76 (MQ=255)
                                    tttAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTc                                    >  2:4160969/1‑76 (MQ=255)
                                    tttAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTc                                    >  1:2289356/1‑76 (MQ=255)
                                     ttAAGCTGGGGGTGATAGCGGAATAAGGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCg                                   <  2:3922379/76‑1 (MQ=255)
                                     ttAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCg                                   >  2:3550675/1‑76 (MQ=255)
                                     ttAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCg                                   <  1:3129612/76‑1 (MQ=255)
                                     ttAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCAGGTCg                                   >  1:3291355/1‑76 (MQ=255)
                                     ttAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACCAAAGCGGGTCg                                   >  2:3557185/1‑76 (MQ=255)
                                      tAAGGTGCGGGTGAGGGGGGAATAAAGGCGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGc                                  <  2:2842789/76‑1 (MQ=255)
                                      tAAGCTGCTGGTGATGGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGc                                  <  1:300553/76‑1 (MQ=255)
                                      tAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGc                                  <  1:649881/76‑1 (MQ=255)
                                        aGCTGGTGGTGAGAGCGGAATAAGGGGGCCGATAGCGGGGCAAAGCCACCATCTCCGGCCCTAAAGCGGGTCGCAt                                <  2:1037715/76‑1 (MQ=255)
                                        aGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCAt                                <  1:1644509/76‑1 (MQ=255)
                                        aGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCAt                                <  2:4325332/76‑1 (MQ=255)
                                        aGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCAt                                >  2:1796429/1‑76 (MQ=255)
                                         gcggctgGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCCCCGGCACTAAAGCGGGTCGCATg                               <  2:3362059/73‑1 (MQ=255)
                                          ctgctgGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGa                              >  1:5443335/1‑76 (MQ=255)
                                           tgctgGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGAc                             >  1:4044151/1‑76 (MQ=255)
                                            gctgGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGAcc                            <  1:4480112/76‑1 (MQ=255)
                                            gctgGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGAcc                            <  1:3741345/76‑1 (MQ=255)
                                             gtgGTGATAGCGAAATAAAGGATCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCg                           <  1:409992/75‑1 (MQ=255)
                                             ctgGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCg                           >  1:2898372/1‑76 (MQ=255)
                                             ctgGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCg                           >  2:5218209/1‑76 (MQ=255)
                                              tgGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGc                          <  1:4026817/76‑1 (MQ=255)
                                              tgGTGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGc                          >  2:3745178/1‑76 (MQ=255)
                                               ggTAATAGCGGAAGAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATCACCGcc                         <  1:375082/76‑1 (MQ=255)
                                                gTGATAGGGAAATAAAGGAGCCGATGGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCa                        <  1:5669496/76‑1 (MQ=255)
                                                 tGATAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAg                       >  1:2314780/1‑76 (MQ=255)
                                                 tGATAGCGGAATAAAGGAGCCAAAACCGGGCAAAACCCACCATCCCCGGCACTAAACCGGGCCCCATGCCCGcaac                       >  2:1400067/1‑73 (MQ=255)
                                                 tGATAGCGAAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAg                       <  1:4283838/76‑1 (MQ=255)
                                                    tAGCGGAATAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACCAAAGCGGGCCGCATGACCGCCACCGt                    >  2:258430/1‑76 (MQ=255)
                                                      ggggAAAAAGGGGGCGGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGGGGGTCGCATGACCGCCAGCGTTg                  <  2:3094627/74‑1 (MQ=255)
                                                      ggggAAAAAAGGGAGCAGATGGCGGGGAAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGtt                   <  1:438927/74‑1 (MQ=255)
                                                      agggAAAAAAGGGAGCCGATAGGGGGGAAAAGCCCCCATCCCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGtt                   <  1:2993424/74‑1 (MQ=255)
                                                       cgcaaTAAAGGAGCCGATAGCCGGGAAAAGCCACCATCTCCGGCGCTAAAGGGCGTGGCATCACAGCCAGCGCTGc                 >  2:4390881/4‑76 (MQ=255)
                                                          aaTAAAGGAGCCGATAGAGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGAGTTGTggg              >  1:4556240/1‑76 (MQ=255)
                                                          aaTAAAGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCggg              >  1:1687711/1‑76 (MQ=255)
                                                             aaaGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcga           >  1:2042237/1‑76 (MQ=255)
                                                              aaGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgat          >  1:2586683/1‑76 (MQ=255)
                                                              aaGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgag          >  2:3672543/1‑75 (MQ=255)
                                                              aaGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgag          >  1:986464/1‑75 (MQ=255)
                                                               aGGAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgatg         >  2:3048783/1‑76 (MQ=255)
                                                                  gAGCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgatgcg       >  1:1156027/1‑76 (MQ=255)
                                                                  gAGCCGATAGCGGGGCAAAGCCACCATCGCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGCTGCGGgcgatgcg       >  2:4328/1‑76 (MQ=255)
                                                                    gCCGATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgatgcgat     >  1:5007983/1‑76 (MQ=255)
                                                                      cGATAGCGGGGCAAAGCCACCATCTCCGGCACTaaa                                           <  2:1230618/36‑1 (MQ=37)
                                                                      cGATAGCGGGGCAAAGCCACCATCTCCGGCACTaaa                                           >  1:1230618/1‑36 (MQ=37)
                                                                      cGAGAGCGGGAAAAACCCACCATCCCCGGCACAAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATgg   <  1:5612288/76‑1 (MQ=255)
                                                                       gATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGa  >  1:2753602/1‑76 (MQ=255)
                                                                       gATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGCTGCGGGCGCGGCGACggg  >  2:2076298/1‑75 (MQ=255)
                                                                       gATAGCGGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGCCCCCTGACCGCCCGCGCTCCGGGCGATGCGCTGGa  >  2:5704895/1‑76 (MQ=255)
                                                                               |                                                                      
ACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGGCAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGA  >  NC_000913/3123507‑3123651

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 4 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: