Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,068,638 T→A 12.4% K384N (AAA→AAT) ‡ aer ← fused signal transducer for aerotaxis sensory component and methyl accepting chemotaxis component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,068,6380TA12.4% 170.7 / 14.9 69K384N (AAA→AAT) ‡aerfused signal transducer for aerotaxis sensory component and methyl accepting chemotaxis component
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base A (5/3);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCA  >  minE/2068566‑2068707
                                                                        |                                                                     
gCGAATATCGTTCGCAGCATTAGCGCTGCGGCTGGCTAAATGA‑‑CACTTCCCCTGCCACCACTGCAAAACCt                                                                       >  1:5477889/1‑71 (MQ=255)
   aaTATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCtt                                                                      <  1:5128433/71‑1 (MQ=255)
    atatCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCttt                                                                     >  1:642705/1‑71 (MQ=255)
    atatCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCttt                                                                     <  1:3440332/71‑1 (MQ=255)
     tatCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTg                                                                    <  1:3794506/71‑1 (MQ=255)
     tagcGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCttt                                                                     <  1:1950387/67‑1 (MQ=255)
       tCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCttt                                                                     >  1:1496479/1‑68 (MQ=255)
       tCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCttt                                                                     >  1:762290/1‑68 (MQ=255)
       tCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCttt                                                                     >  1:2017912/1‑68 (MQ=255)
       tCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGcc                                                                  <  1:5241137/71‑1 (MQ=255)
       tCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGcc                                                                  >  1:1731936/1‑71 (MQ=255)
          ttGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTg                                                               <  1:509632/71‑1 (MQ=255)
          ttGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTg                                                               <  1:2189034/71‑1 (MQ=255)
            ggCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGtt                                                             >  1:1260742/1‑71 (MQ=255)
            ggCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGtt                                                             >  1:2271021/1‑71 (MQ=255)
                  aTTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCgg                                                       <  1:1934680/71‑1 (MQ=255)
                   ttAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGcc                                                         <  1:2562609/68‑1 (MQ=255)
                   ttAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCgg                                                       <  1:4988295/70‑1 (MQ=255)
                    tAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCa                                                     >  1:4627060/1‑71 (MQ=255)
                     aGCGCTGCGGCTGGCTAAATGCCGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCg                                                        >  1:3394627/1‑67 (MQ=255)
                     aGCGCTGCGGCTGGCTAAATGCCGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCAc                                                    >  1:1497492/1‑71 (MQ=255)
                     aGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCg                                                        >  1:288336/1‑67 (MQ=255)
                       cgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCAcgc                                                  <  1:2795825/71‑1 (MQ=255)
                        gcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCAcgcg                                                 <  1:3653829/71‑1 (MQ=255)
                         cTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACgcgc                                                <  1:31476/71‑1 (MQ=255)
                          tacGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCgagc                                                >  1:949004/3‑70 (MQ=255)
                            cGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGc                                                                   >  1:4130272/1‑49 (MQ=255)
                            cGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGc                                                          >  1:2561294/1‑58 (MQ=255)
                               ctggctAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCGTc                                          >  1:5514198/1‑71 (MQ=255)
                                   ctAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc                                                            <  1:5691084/49‑1 (MQ=255)
                                    tAAATGACGCACTTCCCCTGCCACCACTGCAAAACCtt                                                                      >  1:45832/1‑38 (MQ=255)
                                    tAAATGACGCACTTCCCCTGCCACCACTGCAAAACCtt                                                                      >  1:4270326/1‑38 (MQ=255)
                                    tAAATGACGCACTTCCCCTGCCACCACTGCAAAACCtt                                                                      >  1:289426/1‑38 (MQ=255)
                                     aaaTGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc                                                            <  1:96151/47‑1 (MQ=255)
                                       aTGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGt                                                              >  1:591358/1‑43 (MQ=255)
                                       aTGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGc                                                          >  1:3498637/1‑47 (MQ=255)
                                        tGACGCACTTCCCCTGCCACCACTGCAAAACCTTTg                                                                    <  1:4226627/36‑1 (MQ=255)
                                        tGACGCACTTCCCCTGCCACCACTGCAAAACCTTTg                                                                    <  1:1648839/36‑1 (MQ=255)
                                        tGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCAccgccg                                  >  1:2579055/1‑70 (MQ=255)
                                        tGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCAccgcc                                   >  1:4293917/1‑69 (MQ=255)
                                          aCGCACTTCCCCTGCCACCACTGCAAAACCTTTGc                                                                   >  1:2773435/1‑35 (MQ=255)
                                          aCGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGcgcc                                               <  1:5137532/55‑1 (MQ=255)
                                          aCGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGcgcc                                               <  1:3152915/55‑1 (MQ=255)
                                          aCGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGcgcc                                               <  1:3223092/55‑1 (MQ=255)
                                           cGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                              <  1:4839485/71‑1 (MQ=255)
                                                ttCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAggg                         >  1:4406751/1‑71 (MQ=255)
                                                 tcaccTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCa                            >  1:3808880/4‑67 (MQ=255)
                                                 tcaccTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGc                        >  1:1607488/4‑71 (MQ=255)
                                                 tCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCa                                        >  1:3851941/1‑55 (MQ=255)
                                                 tCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCa                                        >  1:5501552/1‑55 (MQ=255)
                                                  ccccTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGc                        >  1:3133454/1‑70 (MQ=255)
                                                    ccTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAg                     <  1:3115023/71‑1 (MQ=255)
                                                      tGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                              <  1:1804907/60‑1 (MQ=255)
                                                      tGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGaa                   >  1:4876798/1‑71 (MQ=255)
                                                       gCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGaa                   >  1:551201/1‑70 (MQ=255)
                                                       gCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGaa                   >  1:3128799/1‑70 (MQ=255)
                                                        ccactaCCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGaa                   <  1:1596229/69‑1 (MQ=255)
                                                        ccaccaCTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCAc                                       >  1:1004843/1‑49 (MQ=255)
                                                         caccacTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAatat                <  1:3381062/71‑1 (MQ=255)
                                                         caccacAGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAatat                >  1:1366067/1‑71 (MQ=255)
                                                           ccacTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTg              <  1:3391711/71‑1 (MQ=255)
                                                           ccacCGCAAAACCGCGCCCTTGCTCACCCGCACGCGCAGCCTCAACCGCCGCGTTCAAAGCCAGAATATTg              <  1:5665247/71‑1 (MQ=255)
                                                           ccacCGCAAAACCGCGCCCTTGCTCACCCGCACGCGCAGCCTCAACCGCCGCGTTCAAAGCCAGAATATTg              >  1:2145271/1‑71 (MQ=255)
                                                             acTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCa                            >  1:872275/1‑55 (MQ=255)
                                                             acTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGt            >  1:1418717/1‑71 (MQ=255)
                                                              cTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTc                                          <  1:983823/40‑1 (MQ=255)
                                                              cTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTc           >  1:4854925/1‑71 (MQ=255)
                                                                 cAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAg                           >  1:560931/1‑52 (MQ=255)
                                                                 cAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTg         <  1:3944400/70‑1 (MQ=255)
                                                                  aaaaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAata                 >  1:5113625/1‑61 (MQ=255)
                                                                  aaaaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAata                 >  1:915605/1‑61 (MQ=255)
                                                                    aaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTc                             <  1:4460276/47‑1 (MQ=255)
                                                                    aaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAAc     <  1:3913188/71‑1 (MQ=255)
                                                                    aaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAAc     <  1:310011/71‑1 (MQ=255)
                                                                    aaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAAc     >  1:2298716/1‑71 (MQ=255)
                                                                    aaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAAc     <  1:210974/71‑1 (MQ=255)
                                                                      ccTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACg    >  1:2514079/1‑70 (MQ=255)
                                                                      ccTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGc   <  1:61465/71‑1 (MQ=255)
                                                                      ccTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGc   <  1:1503705/71‑1 (MQ=255)
                                                                       cTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCa  <  1:3548251/71‑1 (MQ=255)
                                                                        |                                                                     
GCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCA  >  minE/2068566‑2068707

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: