Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,068,626 C→T 10.4% V388V (GTG→GTA 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,6260CT10.4% 180.4 / 9.1 67V388V (GTG→GTAaerfused signal transducer for aerotaxis sensory component and methyl accepting chemotaxis component
Reads supporting (aligned to +/- strand):  ref base C (38/22);  new base T (5/2);  total (43/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

AATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTG  >  minE/2068557‑2068695
                                                                     |                                                                     
aatcaGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGccacca                                                                      >  1:1268924/1‑71 (MQ=255)
aatcaGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGccacca                                                                      >  1:3698807/1‑71 (MQ=255)
aatcaGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGccacca                                                                      >  1:4334973/1‑71 (MQ=255)
aatcaGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGccacca                                                                      >  1:493288/1‑71 (MQ=255)
aatcaGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGccacca                                                                      >  1:5227666/1‑71 (MQ=255)
    aGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTg                                                                   >  1:1268404/1‑70 (MQ=255)
    aGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTg                                                                   >  1:5551494/1‑70 (MQ=255)
    aGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTg                                                                   >  1:2227967/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)
                tCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGccacca                                                                      >  1:604336/1‑55 (MQ=255)
                tCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCttt                                                         >  1:2017912/1‑68 (MQ=255)
                tCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCttt                                                         >  1:1496479/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: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)
                            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:4462562/1‑49 (MQ=255)
                              aGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCaaaa                                                              >  1:4640921/1‑49 (MQ=255)
                              aGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCaaaa                                                              >  1:4828554/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:45832/1‑38 (MQ=255)
                                             tAAATGACGCACTTCCCCTGCCACCACTGCAAAACCtt                                                          >  1:289426/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:3152915/55‑1 (MQ=255)
                                                   aCGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGcgcc                                   <  1:3223092/55‑1 (MQ=255)
                                                   aCGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGcgcc                                   <  1:5137532/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)
                                                                     |                                                                     
AATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTG  >  minE/2068557‑2068695

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: