Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,064,080 C→G 6.5% T107T (ACC→ACG 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,0800CG 6.5% 304.3 / 8.6 92T107T (ACC→ACGrpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base C (59/27);  new base G (2/4);  total (61/31)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 109 total aligned reads displayed.

TTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTC  >  Exported/2064010‑2064150
                                                                      |                                                                      
tttCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGGAcc                                                                        >  1:5506179/1‑71 (MQ=255)
tttCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCAcc                                                                        >  1:3976680/1‑71 (MQ=255)
tttCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCAcc                                                                        >  1:129816/1‑71 (MQ=255)
tttCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCAcc                                                                        >  1:1645377/1‑71 (MQ=255)
tttCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCAcc                                                                        >  1:3762845/1‑71 (MQ=255)
tttCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCAcc                                                                        >  1:3651550/1‑71 (MQ=255)
   ccAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGtt                                                                     >  1:678672/1‑71 (MQ=255)
   ccAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGtt                                                                     >  1:2834270/1‑71 (MQ=255)
    cAGCGTGGAATCTTAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTg                                                                    <  1:6506342/71‑1 (MQ=255)
      gCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATATACATGCGTGAAATGGGCACCGTTGaa                                                                  >  1:1744990/1‑71 (MQ=255)
        gTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACt                                                                >  1:4532968/1‑71 (MQ=255)
        gTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACt                                                                <  1:4393210/71‑1 (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)
             aTCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACt                                                                >  1:6254662/1‑66 (MQ=255)
             aTCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACt                                                                >  1:5988735/1‑66 (MQ=255)
             aTCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACt                                                                >  1:5903568/1‑66 (MQ=255)
             aTCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACt                                                                >  1:3232491/1‑66 (MQ=255)
              tCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTg                                                                    <  1:6277889/61‑1 (MQ=255)
              tCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTg                                                                    <  1:3238560/61‑1 (MQ=255)
                tGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGAccc                                                        >  1:1193875/1‑71 (MQ=255)
                  aaaTCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGaa                                                                  >  1:6274242/1‑59 (MQ=255)
                  aaaTCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGaa                                                                  >  1:2118820/1‑59 (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)
                             aCGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGaa                                                                  >  1:1699176/1‑48 (MQ=255)
                             aCGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAAt                                           >  1:4572189/1‑71 (MQ=255)
                              cGACTGACCCGGTACTCATGTACATGCGTGAAATGGGCACCGTTGAACt                                                                <  1:5591224/49‑1 (MQ=255)
                               gACTGACCCGGTACGCATGTACATGCGTGAAATGGGTACCGTTGAACTGTTGACCCGCGAAGGCGAAATTg                                         <  1:3514912/71‑1 (MQ=255)
                               gACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTg                                         <  1:1976413/71‑1 (MQ=255)
                                 cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                                        >  1:1759616/1‑70 (MQ=255)
                                 cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                                        >  1:5067565/1‑70 (MQ=255)
                                 cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                                        >  1:1989583/1‑70 (MQ=255)
                                 cTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                                        >  1:5769790/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:2749929/71‑3 (MQ=255)
                                        ggTCCGCATGTATATGCGCGAGATGGGCTCGGTCGAGCTGCTCACGCGCGAAGGCGAGATCGAGATCGcga                                <  1:4710248/71‑3 (MQ=255)
                                        ggTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTa                                <  1:1136591/71‑1 (MQ=255)
                                            cGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCg                            <  1:2212401/71‑1 (MQ=255)
                                             gCATGTACATGCGCGAGATGGGCGCGGTGTCGCTGCTGACGCGCGAAGGCGAGATCGAGATCGCCAAGCg                            >  1:4794699/1‑70 (MQ=255)
                                                tGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTAtt                        >  1:648202/1‑71 (MQ=255)
                                                 gTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTAtt                        >  1:4802033/1‑70 (MQ=255)
                                                  tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGaaa                                            >  1:5745828/1‑49 (MQ=255)
                                                  tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGaaa                                            >  1:4381273/1‑49 (MQ=255)
                                                  tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGaaa                                            >  1:112708/1‑49 (MQ=255)
                                                  tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGaaa                                            >  1:1163336/1‑49 (MQ=255)
                                                  tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGaaa                                            >  1:6827173/1‑49 (MQ=255)
                                                  tACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGa                      <  1:4693496/71‑1 (MQ=255)
                                                   aCATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTg                       >  1:3931574/1‑69 (MQ=255)
                                                   aCATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTg                       >  1:550127/1‑69 (MQ=255)
                                                   aCATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTg                       >  1:2275490/1‑69 (MQ=255)
                                                   aCATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTg                       >  1:1748459/1‑69 (MQ=255)
                                                    cATGCGTGAAATGGGCACCGTTGAACTGTTGACCcgcg                                                     <  1:2722520/38‑1 (MQ=255)
                                                    cATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTa                          >  1:3862035/1‑65 (MQ=255)
                                                       gCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACg                 >  1:5501822/1‑71 (MQ=255)
                                                       gCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACg                 >  1:3237163/1‑71 (MQ=255)
                                                       gCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACg                 >  1:4387304/1‑71 (MQ=255)
                                                        cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg                                                  >  1:5343440/1‑37 (MQ=255)
                                                        cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg                                                  >  1:1877126/1‑37 (MQ=255)
                                                        cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg                                                  >  1:6608286/1‑37 (MQ=255)
                                                        cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg                                                  >  1:2489841/1‑37 (MQ=255)
                                                         gTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGa                                        <  1:2718488/46‑1 (MQ=255)
                                                         gTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACggg               <  1:4926132/71‑1 (MQ=255)
                                                         gTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACgg                >  1:5131990/1‑70 (MQ=255)
                                                         gTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACgg                >  1:2610697/1‑70 (MQ=255)
                                                           gAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACa                                      >  1:4439339/1‑46 (MQ=255)
                                                           gAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACa                                      >  1:5451769/1‑46 (MQ=255)
                                                            aaaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGAt             >  1:2570542/1‑70 (MQ=255)
                                                            aaaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGAt             >  1:4077446/1‑70 (MQ=255)
                                                             aaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATc            >  1:3640458/1‑70 (MQ=255)
                                                             aaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATc            >  1:5389874/1‑70 (MQ=255)
                                                             aaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATc            >  1:6543466/1‑70 (MQ=255)
                                                             aaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATc            >  1:503257/1‑70 (MQ=255)
                                                               tGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCa           <  1:3682112/69‑1 (MQ=255)
                                                               tGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCa           <  1:6633275/69‑1 (MQ=255)
                                                                  gCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAg      <  1:586824/71‑1 (MQ=255)
                                                                  gCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAg      <  1:4979151/71‑1 (MQ=255)
                                                                     ccGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTg                       <  1:980058/51‑1 (MQ=255)
                                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAg                              <  1:941492/43‑1 (MQ=255)
                                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAg                              <  1:94000/43‑1 (MQ=255)
                                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAg                              <  1:2907331/43‑1 (MQ=255)
                                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAg                              <  1:680362/43‑1 (MQ=255)
                                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAg                              <  1:6797563/43‑1 (MQ=255)
                                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAg                              <  1:2463551/43‑1 (MQ=255)
                                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAg                              <  1:2488474/43‑1 (MQ=255)
                                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTc  <  1:2521589/71‑1 (MQ=255)
                                                                      cGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCATATTGAAGACGGGATCAACCAGGTTc  <  1:4017084/71‑1 (MQ=255)
                                                                      cGTTGAACTGTTGACACGCGAAGGCGAAATTGACATCTCTAAGCGTATTGAAGACGGGATCAACCAGGTTc  <  1:6208232/71‑1 (MQ=255)
                                                                      |                                                                      
TTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTC  >  Exported/2064010‑2064150

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: