Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,279,370 G→T A298S (GCA→TCA)  ghxP → guanine/hypoxanthine permease, high affinity; guanine/hypoxanthine:H+ symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,279,3700GT98.4% 201.1 / ‑3.4 61A298S (GCA→TCA) ghxPguanine/hypoxanthine permease, high affinity; guanine/hypoxanthine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (28/32);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

GATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGC  >  NC_000913/4279273‑4279464
                                                                                                 |                                                                                              
gatgaCGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGACGGCCAGGCGAACCTGCTGGATAAAGATGGGAAGATCCTCCACGGTGGGAAATca                                                                                              >  2:645590/1‑100 (MQ=255)
 atgaCGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATcac                                                                                             <  1:78767/100‑1 (MQ=255)
  tgaCGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATcact                                                                                            <  1:1418324/100‑1 (MQ=255)
   gaCGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATcactg                                                                                           <  1:1138908/100‑1 (MQ=255)
   gaCGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATcactg                                                                                           >  2:1069695/1‑100 (MQ=255)
     cGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATcactgac                                                                                         <  1:570356/100‑1 (MQ=255)
      ggCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATcactgacc                                                                                        <  1:1443057/100‑1 (MQ=255)
      ggCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATcactgacc                                                                                        <  1:1468552/100‑1 (MQ=255)
       gCAGTATTTGATGCCACCGGAACTAGCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATcactgacca                                                                                       <  2:1292212/100‑1 (MQ=255)
        cAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATcactgaccac                                                                                      >  1:723627/1‑100 (MQ=255)
          gTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATcactgaccactg                                                                                    >  1:652919/1‑100 (MQ=255)
          gTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGcactgaccactg                                                                                    <  1:510203/100‑1 (MQ=255)
             tttGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACt                                                                                 <  1:234256/100‑1 (MQ=255)
                gATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCa                                                                              <  2:577757/100‑1 (MQ=255)
                  tGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATg                                                                            <  2:1086029/100‑1 (MQ=255)
                   gCCACCAGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTAGATAAAGATGGGCAGATCATCGACGGCGAGAAATCACTGACCACTGACTCCATGa                                                                           >  2:316872/1‑100 (MQ=38)
                    ccACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTTGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGag                                                                          <  2:1008603/100‑1 (MQ=255)
                    ccACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATgtg                                                                          >  1:817713/1‑98 (MQ=255)
                    ccACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGag                                                                          <  1:860570/100‑1 (MQ=255)
                    ccACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGag                                                                          >  1:439998/1‑100 (MQ=255)
                     cACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGagc                                                                         <  2:999270/100‑1 (MQ=255)
                        cGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGagcagc                                                                      >  1:1364151/1‑100 (MQ=255)
                        cGGAACTATCCGTGCCGTCGCCGGCCAAGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGagcagc                                                                      >  1:109803/1‑100 (MQ=255)
                            aCTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGTAAATCACTGACCACTGACTCCATGAGCAGCGttt                                                                  <  2:460128/100‑1 (MQ=255)
                             cTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGtttt                                                                 >  2:126020/1‑100 (MQ=255)
                                 ccGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTctct                                                             <  1:291692/100‑1 (MQ=255)
                                  cGTGCCGTCGCCGGACAGGCGAACCTACTGGATAAAGAGGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTg                                                            <  1:1204242/100‑1 (MQ=255)
                                   gTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTgg                                                           <  1:363162/100‑1 (MQ=255)
                                     gCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTTACTCCATGAGCAGCGTTTTCTCTgg                                                           <  2:976423/98‑1 (MQ=255)
                                      ccGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCt                                                        >  1:1267635/1‑100 (MQ=255)
                                      ccGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCt                                                        >  1:1345364/1‑100 (MQ=255)
                                          cGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCtggtg                                                    <  2:1238389/100‑1 (MQ=255)
                                            ccGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTggtggg                                                  <  1:345397/100‑1 (MQ=255)
                                             cGGCCAGGCGAACCTGCTGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTggtgggtg                                                >  2:1299818/1‑100 (MQ=255)
                                              ggCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTggtgggtg                                                >  2:451218/1‑100 (MQ=255)
                                                ccAGGCGAACCTGCTGCGTAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCg                                              <  2:924420/100‑1 (MQ=255)
                                                 cAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCgg                                             >  2:837729/1‑100 (MQ=255)
                                                 cAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCATGGTGGGTGCgg                                             >  1:307213/1‑100 (MQ=255)
                                                   ggCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTTGGAAATCACTGACCACTGACTCCATGTGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCt                                           <  2:779719/100‑1 (MQ=255)
                                                   ggCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCt                                           >  2:93108/1‑100 (MQ=255)
                                                      gAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCg                                        <  1:146325/100‑1 (MQ=255)
                                                      gAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCg                                        <  2:531984/100‑1 (MQ=255)
                                                      gAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCg                                        <  1:1004740/100‑1 (MQ=255)
                                                       aaCCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCgg                                       <  1:1447496/100‑1 (MQ=255)
                                                           tgctgGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGc                                   >  2:499031/1‑100 (MQ=255)
                                                            gctgGATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGTGGCTCCGGCAGCg                                  <  1:1060873/100‑1 (MQ=255)
                                                               ggATAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTGTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGta                               <  2:307213/100‑1 (MQ=255)
                                                                 aTAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGtata                             <  2:458412/100‑1 (MQ=255)
                                                                 aTAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGtata                             >  2:1341382/1‑100 (MQ=255)
                                                                  tAAAGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATAc                            <  1:5407/100‑1 (MQ=255)
                                                                     aGATGGGCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATc                         <  2:1364144/100‑1 (MQ=255)
                                                                         gggCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGt                     >  2:1408447/1‑100 (MQ=255)
                                                                         gggCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGt                     >  2:888819/1‑100 (MQ=255)
                                                                          ggCAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTc                    <  2:109803/100‑1 (MQ=255)
                                                                            cAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTg                  >  1:461310/1‑100 (MQ=255)
                                                                            cAGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTg                  >  2:72459/1‑100 (MQ=255)
                                                                             aGATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTgc                 >  1:392283/1‑100 (MQ=255)
                                                                              gATCATCGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTgcg                >  2:700354/1‑100 (MQ=255)
                                                                                   tcGACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCggg           >  1:1052845/1‑100 (MQ=255)
                                                                                     gACGGTGGGAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTa         >  2:550675/1‑100 (MQ=255)
                                                                                            ggAAATCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAcggcggc  >  2:45890/1‑100 (MQ=255)
                                                                                                 |                                                                                              
GATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGC  >  NC_000913/4279273‑4279464

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