Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,068,629 T→C 13.5% A387A (GCA→GCG 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,6290TC13.5% 184.5 / 16.9 77A387A (GCA→GCGaerfused signal transducer for aerotaxis sensory component and methyl accepting chemotaxis component
Reads supporting (aligned to +/- strand):  ref base T (41/25);  new base C (6/4);  total (48/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTC  >  minE/2068561‑2068698
                                                                    |                                                                     
aGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTg                                                                      >  1:1268404/1‑70 (MQ=255)
aGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTg                                                                      >  1:2227967/1‑70 (MQ=255)
aGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTg                                                                      >  1:5551494/1‑70 (MQ=255)
aGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCAATTCCCCTGCCACCACTg                                                                      >  1:4773202/1‑70 (MQ=255)
   ttGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCaa                                                                   >  1:3401501/1‑70 (MQ=255)
     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:2017912/1‑68 (MQ=255)
            tCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCttt                                                            >  1:762290/1‑68 (MQ=255)
            tCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGcc                                                         >  1:1731936/1‑71 (MQ=255)
            tCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGcc                                                         <  1:5241137/71‑1 (MQ=255)
               ttGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTg                                                      <  1:2189034/71‑1 (MQ=255)
               ttGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTg                                                      <  1:509632/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)
                        ttAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCaaa                                                                  <  1:4201143/50‑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)
                          aGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACt                                                                       >  1:650250/1‑43 (MQ=255)
                          aGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCaaaa                                                                 >  1:4828554/1‑49 (MQ=255)
                          aGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCaaaa                                                                 >  1:4640921/1‑49 (MQ=255)
                          aGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCaaaa                                                                 >  1:4462562/1‑49 (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)
                                 cGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCaaaa                                                                 >  1:2461944/1‑42 (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:289426/1‑38 (MQ=255)
                                         tAAATGACGCACTTCCCCTGCCACCACTGCAAAACCtt                                                             >  1:45832/1‑38 (MQ=255)
                                         tAAATGACGCACTTCCCCTGCCACCACTGCAAAACCtt                                                             >  1:4270326/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:2145271/1‑71 (MQ=255)
                                                                ccacCGCAAAACCGCGCCCTTGCTCACCCGCACGCGCAGCCTCAACCGCCGCGTTCAAAGCCAGAATATTg     <  1:5665247/71‑1 (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)
                                                                    |                                                                     
AGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTC  >  minE/2068561‑2068698

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: