Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,093,239 C→A 8.6% G207C (GGT→TGT)  gntT ← gluconate transporter, high‑affinity GNT I system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,093,2390CA 8.6% 224.1 / 4.8 69G207C (GGT→TGT) gntTgluconate transporter, high‑affinity GNT I system
Reads supporting (aligned to +/- strand):  ref base C (33/30);  new base A (2/4);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

AGAGGTCCAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTCGG  >  AP009048/4093171‑4093308
                                                                    |                                                                     
agagGTCCAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCtt                                                                     <  1:2069556/71‑1 (MQ=255)
    gTCCAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCt                                                                      <  1:2521722/66‑1 (MQ=255)
       cAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAAt                                                              >  1:161298/1‑71 (MQ=255)
       cAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAGGTTTTCGCGCTGTAGAGACATTCCGGAAtt                                                             <  1:2684691/71‑1 (MQ=255)
        aGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAAtt                                                             <  1:2391335/71‑1 (MQ=255)
        aGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCGGAAtt                                                             <  1:4204637/71‑1 (MQ=255)
         gACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTg                                                            >  1:3255493/1‑71 (MQ=255)
          aCGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTgg                                                           >  1:3439028/1‑71 (MQ=255)
            gCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGaa                                                               >  1:977623/1‑65 (MQ=255)
            gCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCt                                                         <  1:650564/71‑1 (MQ=255)
            gCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCt                                                         <  1:1845803/71‑1 (MQ=255)
              tGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAAtt                                                             >  1:3376058/1‑65 (MQ=255)
                aCGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTa                                                       <  1:1732867/69‑1 (MQ=255)
                 cGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATc                                                     >  1:4291340/1‑70 (MQ=255)
                  gCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGa                                                   >  1:4609495/1‑71 (MQ=255)
                   ccAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGa                                                   >  1:1250959/1‑70 (MQ=255)
                    cAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATa                                                 >  1:2927984/1‑71 (MQ=255)
                    cAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATa                                                 >  1:1261844/1‑71 (MQ=255)
                    cAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATa                                                 >  1:1408459/1‑71 (MQ=255)
                     aaaGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATAc                                                >  1:3414112/1‑71 (MQ=255)
                      aaGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCTCTGTAGAGACCTTACGGAATTGGCTTATCg                                                    >  1:4470787/1‑66 (MQ=255)
                      aaGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATAcc                                               >  1:976529/1‑71 (MQ=255)
                      aaGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATAcc                                               <  1:4483552/71‑1 (MQ=255)
                       aGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCt                                              >  1:2243704/1‑71 (MQ=255)
                       aGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCt                                              >  1:190639/1‑71 (MQ=255)
                        gCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTcc                                                                   <  1:2834792/49‑1 (MQ=255)
                             gCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAac                                        <  1:2794584/71‑1 (MQ=255)
                             gCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCGGAATTGGCTTATCGATACCTTTCaa                                         >  1:4159512/1‑70 (MQ=255)
                              cATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAaca                                       >  1:1200231/1‑71 (MQ=255)
                               aTCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAacac                                      >  1:1002755/1‑71 (MQ=255)
                                tctcTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAAtt                                                             <  1:1659902/47‑1 (MQ=255)
                                  tcttctTCGCTGAAGGTTTTCGCGCTGTAGAGACATTCCGGAATTGGCTTATCGATACCTTTCAACAcgc                                    >  1:3008377/1‑70 (MQ=255)
                                    ttcttcGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGaa                                                               <  1:77249/41‑1 (MQ=255)
                                    ttcttcGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATc                                                     <  1:3658003/51‑1 (MQ=255)
                                    ttcttcGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAacac                                      >  1:2183459/1‑66 (MQ=255)
                                    ttcttcGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACAcgc                                    <  1:1988532/68‑1 (MQ=255)
                                    ttcatcGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAg                                 >  1:1623254/1‑71 (MQ=255)
                                     tcttcGCTGAAGGTTTTCGCGCTGTAGAGACCTTGCGGAATTGGCTTATCGATACCTTTCAACACGCGAGc                                >  1:4157900/1‑71 (MQ=255)
                                     tcttcGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCgg                                                                 <  1:1811896/38‑1 (MQ=255)
                                      cttcGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGa                                  <  1:366000/68‑1 (MQ=255)
                                      cttcGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCg                               >  1:2564704/1‑71 (MQ=255)
                                        tcGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTc                                           >  1:3606111/1‑57 (MQ=255)
                                         cGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGt                              <  1:718226/69‑1 (MQ=255)
                                          gCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTaaa                           <  1:71254/71‑1 (MQ=255)
                                          gCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTaaa                           <  1:76735/71‑1 (MQ=255)
                                          gCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTaa                            <  1:951521/70‑1 (MQ=255)
                                            tGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTc                                           <  1:2018292/53‑1 (MQ=255)
                                             gAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAaccg                        >  1:14302/1‑71 (MQ=255)
                                              aaGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGa                                                   >  1:3524823/1‑43 (MQ=255)
                                                ggTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAaccggac                     >  1:2431388/1‑71 (MQ=255)
                                                ggTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAaccgga                      >  1:496670/1‑70 (MQ=255)
                                                 gTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAaccggacc                    <  1:2746481/71‑1 (MQ=255)
                                                  ttttCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAaccggacc                    <  1:3658524/70‑1 (MQ=255)
                                                    ttCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAaccggaccgg                  >  1:1151735/1‑70 (MQ=255)
                                                     tCGCGCTGTAGAGACATTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGc                 <  1:4443956/70‑1 (MQ=255)
                                                       gcgcTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTaa                            <  1:2676189/57‑1 (MQ=255)
                                                          cTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAaccgg                       >  1:3142748/1‑59 (MQ=255)
                                                          cTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAAt           <  1:3061771/71‑1 (MQ=255)
                                                             tAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAAt           >  1:4494890/1‑68 (MQ=255)
                                                             tAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCAc        <  1:4345893/71‑1 (MQ=255)
                                                               gagaCCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCga               <  1:4525121/62‑1 (MQ=255)
                                                               gagaCCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACgg      >  1:431669/1‑71 (MQ=255)
                                                               gagaCCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACg       <  1:3581929/70‑1 (MQ=255)
                                                                 gaCCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAacc                         >  1:2210381/1‑50 (MQ=255)
                                                                 gaCCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTc    <  1:3615053/71‑1 (MQ=255)
                                                                 gaCCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTc    >  1:196790/1‑71 (MQ=255)
                                                                 gaCATTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTc    <  1:1453850/71‑1 (MQ=255)
                                                                  aCCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTaa                            <  1:4353998/46‑1 (MQ=255)
                                                                  aCCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAaccggaccgg                  <  1:381149/56‑1 (MQ=255)
                                                                    cTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCg                >  1:1115708/1‑56 (MQ=255)
                                                                    cTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTCgg  >  1:2846543/1‑70 (MQ=255)
                                                                    cTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTCgg  >  1:2124182/1‑70 (MQ=255)
                                                                    |                                                                     
AGAGGTCCAGACGCTGACGCCAAAGCTCGGCATCTCTTCTTCGCTGAAGGTTTTCGCGCTGTAGAGACCTTCCGGAATTGGCTTATCGATACCTTTCAACACGCGAGCGTAAACCGGACCGGCGAGAATCACGGTCGG  >  AP009048/4093171‑4093308

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