Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,590,087 A→T 14.6% intergenic (‑763/+36) gadE ← / ← hdeD DNA‑binding transcriptional activator/acid‑resistance membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,590,0870AT14.6% 137.2 / 16.2 55intergenic (‑763/+36)gadE/hdeDDNA‑binding transcriptional activator/acid‑resistance membrane protein
Reads supporting (aligned to +/- strand):  ref base A (27/20);  new base T (4/4);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

TTTCAATGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTG  >  Exported/2590017‑2590156
                                                                      |                                                                     
tttCAATGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTcaca                                                                       <  1:2409530/71‑1 (MQ=255)
   cAATGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATc                                                                     >  1:3009801/1‑70 (MQ=255)
   cAATGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATc                                                                     >  1:2900653/1‑70 (MQ=255)
   cAATGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATc                                                                     >  1:1356221/1‑70 (MQ=255)
   cAATGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATc                                                                     >  1:1475256/1‑70 (MQ=255)
    aaTGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCaa                                                                   <  1:1774916/71‑1 (MQ=255)
      tGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGt                                                                 <  1:3697609/71‑1 (MQ=255)
      tGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGt                                                                 >  1:1508633/1‑71 (MQ=255)
         tAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACTTGATGTGAc                                                              <  1:551641/71‑1 (MQ=38)
          aGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGAc                                                              >  1:958681/1‑70 (MQ=255)
            tGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTAt                                                           >  1:1770637/1‑71 (MQ=255)
              cgcAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATg                                                          <  1:1733957/70‑1 (MQ=255)
              cgcAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGa                                                         >  1:2865140/1‑71 (MQ=255)
                taaTTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATc                                                       <  1:3554410/70‑1 (MQ=255)
                cAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATc                                                       <  1:1289165/71‑1 (MQ=255)
                 aaTTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATcc                                                      >  1:2736671/1‑71 (MQ=255)
                 aaTTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATcc                                                      >  1:2387378/1‑71 (MQ=255)
                 aaTTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATcc                                                      <  1:1887424/71‑1 (MQ=255)
                    tGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACTTGATGTGACTAt                                                           <  1:413759/63‑1 (MQ=25)
                    tGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACTTGATGTGACTAt                                                           <  1:1348074/63‑1 (MQ=25)
                     gATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCggg                                                   <  1:1135893/70‑1 (MQ=255)
                      aTCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACTTGATGTGACTATGAt                                                        >  1:1482704/1‑64 (MQ=25)
                      aTCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACTTGATGTGACTATGAt                                                        >  1:174023/1‑64 (MQ=25)
                      aTCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGt                                                  <  1:2937992/70‑1 (MQ=255)
                           caacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACAt                                                                      >  1:3240357/1‑45 (MQ=255)
                           caacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACAt                                                                      >  1:2430534/1‑45 (MQ=255)
                           caacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACAt                                                                      >  1:3550635/1‑45 (MQ=255)
                           caacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACAt                                                                      >  1:2635906/1‑45 (MQ=255)
                           caacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACa                                              <  1:2441114/69‑1 (MQ=255)
                            aacaCTGCGTAGCGGAGAGAGTATTAATCTGATCATAGTCACATc                                                                     >  1:1500100/1‑45 (MQ=255)
                            aacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAAcc                                           >  1:3008340/1‑71 (MQ=255)
                            aacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAAcc                                           >  1:635911/1‑71 (MQ=255)
                              cacTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTcaca                                                                       <  1:1228488/41‑1 (MQ=255)
                              cacTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACTTGATGTGACt                                                             <  1:1598542/51‑1 (MQ=21)
                              cacTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCg                                          >  1:1481494/1‑70 (MQ=255)
                               acTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGt                                                                 <  1:2813525/46‑1 (MQ=255)
                                cTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCgggg                                       >  1:1820606/1‑71 (MQ=255)
                                  gCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCgg                                                    <  1:360923/56‑1 (MQ=255)
                                   cGTAGCGGAGAGAGTATTAATCGGATCATAGTCACTTGATGTGACTATGa                                                         >  1:2501433/1‑50 (MQ=18)
                                   cGTAGCGGAGAGAGTATTAATCGGATCATAGTCACTTGATGTGACTATGa                                                         >  1:1929548/1‑50 (MQ=18)
                                    gTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAat                                   <  1:352939/71‑1 (MQ=255)
                                     tAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAatt                                  <  1:80342/71‑1 (MQ=255)
                                       gCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAattat                                <  1:2055885/71‑1 (MQ=255)
                                       gCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGAGACTATGATCCGGGTGACAACCGGGGTAattat                                >  1:1306696/1‑71 (MQ=255)
                                          gagagaGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATtgc                             <  1:3732275/71‑1 (MQ=255)
                                            gagagTATTAATCGGATCATAGTCACATCAAGTGACTa                                                            >  1:3349736/1‑38 (MQ=255)
                                             agagTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATtgctgc                          >  1:1446948/1‑71 (MQ=255)
                                                gTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTa                       <  1:2124720/71‑1 (MQ=255)
                                                 tATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTaa                      >  1:2770297/1‑71 (MQ=255)
                                                     aaTCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGa                   >  1:1868235/1‑70 (MQ=255)
                                                      aTCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGaa                  >  1:3273105/1‑70 (MQ=255)
                                                         ggATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACaaa              <  1:1455562/71‑1 (MQ=255)
                                                          gATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAAc             <  1:2632683/71‑1 (MQ=255)
                                                          gATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAAc             <  1:1171747/71‑1 (MQ=255)
                                                             cATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCtt                        <  1:3148782/57‑1 (MQ=255)
                                                                aGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCg        >  1:930581/1‑70 (MQ=255)
                                                                 gTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCg        >  1:2874615/1‑69 (MQ=255)
                                                                    acaTCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGc    <  1:3252652/70‑1 (MQ=255)
                                                                    acaTCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGc    <  1:597434/70‑1 (MQ=255)
                                                                    acaTCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCt   <  1:2859995/71‑1 (MQ=255)
                                                                    acaTCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCt   <  1:2065624/71‑1 (MQ=255)
                                                                     caTCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTg  <  1:3632356/71‑1 (MQ=255)
                                                                      |                                                                     
TTTCAATGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTG  >  Exported/2590017‑2590156

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: