Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,064,071 A→G 6.1% E104E (GAA→GAG rpoD → RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,064,0710AG 6.1% 290.6 / 8.6 99E104E (GAA→GAGrpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base A (68/25);  new base G (2/4);  total (70/29)
Fisher's exact test for biased strand distribution p-value = 5.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 116 total aligned reads displayed.

CAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATC  >  Exported/2064003‑2064140
                                                                    |                                                                     
cAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAAt                                                                     <  1:5415987/71‑1 (MQ=255)
     gCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                                                      >  1:3359160/1‑65 (MQ=255)
     gCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                                                      >  1:1237863/1‑65 (MQ=255)
     gCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                                                      >  1:1405950/1‑65 (MQ=255)
     gCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                                                      >  1:4220715/1‑65 (MQ=255)
     gCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCa                                                                >  1:2896334/1‑71 (MQ=255)
     gCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCa                                                                >  1:6871811/1‑71 (MQ=255)
     gCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCa                                                                >  1:4131072/1‑71 (MQ=255)
       tttCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGGAcc                                                              >  1:5506179/1‑71 (MQ=255)
       tttCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCAcc                                                              >  1:3762845/1‑71 (MQ=255)
       tttCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCAcc                                                              >  1:1645377/1‑71 (MQ=255)
       tttCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCAcc                                                              >  1:3976680/1‑71 (MQ=255)
         tCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGt                                                            >  1:4299813/1‑71 (MQ=255)
          ccAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGtt                                                           >  1:2834270/1‑71 (MQ=255)
           cAGCGTGGAATCTTAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTg                                                          <  1:6506342/71‑1 (MQ=255)
            aGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTg                                                          >  1:3250635/1‑70 (MQ=255)
             gCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATATACATGCGTGAAATGGGCACCGTTGaa                                                        >  1:1744990/1‑71 (MQ=255)
               gTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACt                                                      >  1:4532968/1‑71 (MQ=255)
               gTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACt                                                      <  1:5583495/71‑1 (MQ=255)
               gTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACt                                                      <  1:3010090/71‑1 (MQ=255)
               gTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGATGAACt                                                      <  1:886874/71‑1 (MQ=255)
                tGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTg                                                     >  1:823330/1‑71 (MQ=255)
                   aaTCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTg                                                  <  1:1877256/71‑1 (MQ=255)
                    aTCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGct                                                                >  1:4345447/1‑55 (MQ=255)
                    aTCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACt                                                      >  1:5903568/1‑66 (MQ=255)
                    aTCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACt                                                      >  1:6254662/1‑66 (MQ=255)
                     tCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTg                                                          <  1:6277889/61‑1 (MQ=255)
                     tCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTg                                                          <  1:3238560/61‑1 (MQ=255)
                       tGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                                                      >  1:4902331/1‑47 (MQ=255)
                       tGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                                                      >  1:6530002/1‑47 (MQ=255)
                       tGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                                                      >  1:5740928/1‑47 (MQ=255)
                       tGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                                                      >  1:5881589/1‑47 (MQ=255)
                       tGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                                                      >  1:1212482/1‑47 (MQ=255)
                       tGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                                                      >  1:37330/1‑47 (MQ=255)
                       tGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAAt                                                                     >  1:4634645/1‑48 (MQ=255)
                       tGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAAt                                                                     >  1:2442416/1‑48 (MQ=255)
                       tGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAAt                                                                     >  1:6790396/1‑48 (MQ=255)
                       tGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGAccc                                              >  1:1193875/1‑71 (MQ=255)
                       tGAAATCGGGCGCACGACTGACACGGTACGCATGTACATGCGTGaaa                                                                      >  1:718573/1‑47 (MQ=255)
                         aaaTCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGaa                                                        >  1:2118820/1‑59 (MQ=255)
                         aaaTCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGaa                                                        >  1:6274242/1‑59 (MQ=255)
                         aaaTCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCcgc                                            >  1:656230/1‑71 (MQ=255)
                          aaTCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCcgc                                            >  1:5743850/1‑70 (MQ=255)
                          aaTCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCcgc                                            >  1:3110790/1‑70 (MQ=255)
                             cGGCCGCACCACCGACCCGGTCCGCATGTACATGCGCGAGATGGGCGCGGTGTCGCTGCTGACGCGCGAAg                                        >  1:4729303/1‑71 (MQ=255)
                              gggCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAgg                                       >  1:5016737/1‑71 (MQ=255)
                                gcgcACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGc                                      <  1:6257689/70‑1 (MQ=255)
                                 cgcACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGa                                    <  1:6080671/71‑1 (MQ=255)
                                    aCGACTGACCCGGTACGCATGTACATGCGTGAAATg                                                                    >  1:1137915/1‑36 (MQ=255)
                                    aCGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGaa                                                        >  1:1699176/1‑48 (MQ=255)
                                    aCGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAAt                                 >  1:4572189/1‑71 (MQ=255)
                                      gACTGACCCGGTACGCATGTACATGCGTGAAATGGGTACCGTTGAACTGTTGACCCGCGAAGGCGAAATTg                               <  1:3514912/71‑1 (MQ=255)
                                      gACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTg                               <  1:1976413/71‑1 (MQ=255)
                                        cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                              >  1:1950799/1‑70 (MQ=255)
                                        cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                              >  1:1759616/1‑70 (MQ=255)
                                        cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                              >  1:2579748/1‑70 (MQ=255)
                                        cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                              >  1:5769790/1‑70 (MQ=255)
                                        cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                              >  1:1989583/1‑70 (MQ=255)
                                        cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGAc                             <  1:6385041/71‑1 (MQ=255)
                                           aCCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATc                          >  1:519326/1‑71 (MQ=255)
                                            cccGGTCCGCATGTACATGCGCGAGATGGGCGCGGTGTCGCTGCTGACGCGCGAAGGCGagatcgagatcg                         <  1:3969888/71‑1 (MQ=255)
                                            cccGGTCCGCATGTACATGCGCGAGATGGGCGCGGTGTCGCTGCTGACGCGCGAAGGCGagatcgagatcg                         <  1:2558098/71‑1 (MQ=255)
                                            cccGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCg                         <  1:5218307/71‑1 (MQ=255)
                                             ccGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGc                        >  1:5151413/1‑71 (MQ=255)
                                               ggTCCGCATGTATATGCGCGAGATGGGCTCGGTCGAGCTGCTCACGCGCGAAGGCGAGATCGAGATCGcga                      <  1:4710248/71‑3 (MQ=255)
                                               ggTCCGCATGTATATGCGCGAGATGGGCTCGGTCGAGCTGCTCACGCGCGAAGGCGAGATCGAGATCGcga                      <  1:2749929/71‑3 (MQ=255)
                                               ggTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTa                      <  1:1136591/71‑1 (MQ=255)
                                                   cGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCg                  <  1:2212401/71‑1 (MQ=255)
                                                       tGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTAtt              >  1:648202/1‑71 (MQ=255)
                                                        gTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTAtt              >  1:4802033/1‑70 (MQ=255)
                                                         tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGaaa                                  >  1:112708/1‑49 (MQ=255)
                                                         tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGaaa                                  >  1:1163336/1‑49 (MQ=255)
                                                         tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGaaa                                  >  1:6827173/1‑49 (MQ=255)
                                                         tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGaaa                                  >  1:4381273/1‑49 (MQ=255)
                                                         tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGaaa                                  >  1:5745828/1‑49 (MQ=255)
                                                         tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGa            <  1:4693496/71‑1 (MQ=255)
                                                          aCATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTg             >  1:3931574/1‑69 (MQ=255)
                                                          aCATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTg             >  1:1748459/1‑69 (MQ=255)
                                                          aCATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTg             >  1:550127/1‑69 (MQ=255)
                                                           cATGCGTGAAATGGGCACCGTTGAACTGTTGACCcgcg                                           <  1:2722520/38‑1 (MQ=255)
                                                           cATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTa                >  1:3862035/1‑65 (MQ=255)
                                                              gCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACg       >  1:4387304/1‑71 (MQ=255)
                                                              gCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACg       >  1:5501822/1‑71 (MQ=255)
                                                              gCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACg       >  1:3237163/1‑71 (MQ=255)
                                                               cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg                                        >  1:2489841/1‑37 (MQ=255)
                                                               cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg                                        >  1:6608286/1‑37 (MQ=255)
                                                               cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg                                        >  1:5343440/1‑37 (MQ=255)
                                                               cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg                                        >  1:594988/1‑37 (MQ=255)
                                                                gTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                              <  1:2718488/46‑1 (MQ=255)
                                                                gTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACggg     <  1:4926132/71‑1 (MQ=255)
                                                                gTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACggg     <  1:4918071/71‑1 (MQ=255)
                                                                gTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACgg      >  1:5131990/1‑70 (MQ=255)
                                                                gTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACgg      >  1:2610697/1‑70 (MQ=255)
                                                                  gAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACa                            >  1:5451769/1‑46 (MQ=255)
                                                                  gAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACa                            >  1:4439339/1‑46 (MQ=255)
                                                                   aaaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGAt   >  1:2570542/1‑70 (MQ=255)
                                                                    aaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATc  >  1:6543466/1‑70 (MQ=255)
                                                                    aaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATc  >  1:503257/1‑70 (MQ=255)
                                                                    aaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATc  >  1:3640458/1‑70 (MQ=255)
                                                                    aaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATc  >  1:5389874/1‑70 (MQ=255)
                                                                    |                                                                     
CAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATC  >  Exported/2064003‑2064140

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