Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,443,053 T→A 10.3% T4S (ACC→TCC)  ghrB ← 2‑ketoaldonate reductase/hydroxypyruvate/glyoxylate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,443,0530TA10.3% 93.1 / 4.0 39T4S (ACC→TCC) ghrB2‑ketoaldonate reductase/hydroxypyruvate/glyoxylate reductase
Reads supporting (aligned to +/- strand):  ref base T (16/19);  new base A (2/2);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-01

CGCCGATCAGACCATGGGCTTCAGGCAGGGCTTCGTTGAACTGGGCGCTGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAG‑CCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAAT  >  NC_002947/1442909‑1443193
                                                                                                                                                |                                                                                                                                             
cGCCGATCAGACCATGGGCTTCAGGCAGGGCTTCGTTGAACTGGGCGCTGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCttt                                                                                                                                           <  6:128443/149‑1 (MQ=255)
       cAGACCATGGGCTTCAGGCAGGGCTTCGTTGAACTGGGCGCTGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGATCGGCCATGGCCGGGGTGATACGGCTGACGGCCAGGACGGTCTTTTTCATCg                                                                                                                                    >  8:140182/1‑149 (MQ=255)
                     caggcaggGCTTCGTTGAACTGGGCGCTGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGATCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTg                                                                                                                       <  7:140182/148‑1 (MQ=255)
                        gcaggGCTTCGTTGAACTGGGCGCTGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAg                                                                                                                   >  1:184160/1‑149 (MQ=255)
                        gcaggGCTTCGTTGAACTGGGCGCTGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAg                                                                                                                   >  2:281072/1‑149 (MQ=255)
                          aggGCTTCGTTGAACTGGGCGCTGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGc                                                                                                                 <  4:47348/149‑1 (MQ=255)
                             gCTTCGTTGAACTGGGCGCTGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCggg                                                                                                              <  2:184160/149‑1 (MQ=255)
                              cTTCGTTGAACTGGGCGCTGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCgggg                                                                                                             >  6:65134/1‑149 (MQ=255)
                                 cGTTGAACTGGGCGCTGATGTCGCCGAGCTTGGGGGTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGt                                                                                                            <  6:219898/147‑1 (MQ=255)
                                  gTTGAACTGGGCGCTGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTggg                                                                                                         >  4:46576/1‑149 (MQ=255)
                                   ttGAACTGGGCGCTGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTgggg                                                                                                        >  3:167374/1‑149 (MQ=255)
                                   ttGAACTGGGCGCTGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTgggg                                                                                                        >  3:79198/1‑149 (MQ=255)
                                      atatGGGCGCTGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCtg                                                                                                     <  1:281072/146‑1 (MQ=255)
                                       aCTGGGCGAGGATGTCGCCGAGCTTGGGTTTGGGCAGGATCACGTTGAATTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCtg                                                                                                     <  4:79198/148‑1 (MQ=255)
                                           ggCGCTGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGACTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTg                                                                                                >  2:320594/1‑149 (MQ=255)
                                                tGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTccc                                                                                           >  6:152358/1‑149 (MQ=255)
                                                  aTGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGACTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTccc                                                                                           <  1:320594/147‑1 (MQ=255)
                                                     tCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTg                                                                                                           <  1:84411/128‑1 (MQ=255)
                                                     tCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTg                                                                                                           >  2:84411/1‑128 (MQ=255)
                                                     tCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGACTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCt                                                                                                      <  5:219459/133‑1 (MQ=255)
                                                     tCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGACTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCt                                                                                                      <  5:224169/133‑1 (MQ=255)
                                                     tCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGACTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCt                                                                                                      >  6:219459/1‑133 (MQ=255)
                                                     tCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGACTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCt                                                                                                      >  6:224169/1‑133 (MQ=255)
                                                       gCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCtgt                                                                                                    >  1:172466/1‑133 (MQ=255)
                                                       gCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCtgt                                                                                                    <  2:172466/133‑1 (MQ=255)
                                                       gCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTATCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCgg                                                                                    >  6:306028/1‑149 (MQ=255)
                                                                ggggTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCgg                                                                                    <  1:27638/140‑1 (MQ=255)
                                                                ggggTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCgg                                                                                    >  2:27638/1‑140 (MQ=255)
                                                                       ggCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCggg                                                                                   <  5:154655/134‑1 (MQ=255)
                                                                       ggCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCggg                                                                                   >  6:154655/1‑134 (MQ=255)
                                                                            gATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAG‑CCCGCTCCTACggggg                                                              <  3:294565/149‑1 (MQ=255)
                                                                                aCGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAG‑CCCGCTCCTACGGGGGGgc                                                           <  3:46576/148‑1 (MQ=255)
                                                                                   ttGACGTCTTGCTCCAGGCGTGCGGCCAGGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTATCTTCGGCCAACTTCTCTGCACGGCGGGTTGGGGCTGTGCTGGTCCCTTCGCGGGCAAG‑CCCGCTCCTACggggggggcgcg                                                       <  6:53578/149‑6 (MQ=38)
                                                                                          ctggctgcAGGCGGTCGGCCATGGCCGGGGTGATCCGGCGGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCGGTGCTGGTCCCTTCGCGGGCAAG‑CCCGCTCCTACGGGGGTGCGCGGTGCCTTg                                                <  4:167374/146‑1 (MQ=255)
                                                                                              ctgcAGGCGTGCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGGTGTGCTGGTCCCTTCGCGGGCAAGCCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTa                                              <  1:202421/148‑1 (MQ=255)
                                                                                                 cAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAG‑CCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGc                                         >  5:80633/1‑149 (MQ=255)
                                                                                                   ggCGTTCGGCCATGGCCGGGGTGATCCGCCTGAAGGCCATGACGGTCTTTTTCATCGGGCACCAGCTCTGCAAGGCGGGGGGGGGCTGTGCTGGTCCCTTCGCGGGCAAG‑CCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCgg                                       <  5:306028/149‑1 (MQ=39)
                                                                                                                   cGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTGCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAG‑CCCGCTCCTACAAGGCACCGCGCCCCCCCGTAGGAGCGGGCTTGCCCGCGAAggg                       <  5:152358/149‑2 (MQ=21)
                                                                                                                                        cAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAG‑CCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAAt  >  4:245653/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                             
CGCCGATCAGACCATGGGCTTCAGGCAGGGCTTCGTTGAACTGGGCGCTGATGTCGCCGAGCTTGGGGTTGGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAG‑CCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAAT  >  NC_002947/1442909‑1443193

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: ‑