Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,068,656 A→G 18.2% R378R (CGT→CGC 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
*Exported2,068,6560AG18.2% 139.5 / 23.8 66R378R (CGT→CGCaerfused signal transducer for aerotaxis sensory component and methyl accepting chemotaxis component
Reads supporting (aligned to +/- strand):  ref base A (29/25);  new base G (7/5);  total (36/30)
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

TAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCG  >  Exported/2068586‑2068723
                                                                      |                                                                   
tAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCa                                                                     >  1:4627060/1‑71 (MQ=255)
 aGCGCTGCGGCTGGCTAAATGCCGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCAc                                                                    >  1:1497492/1‑71 (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)
           ctggctAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCGTc                                                          >  1:5514198/1‑71 (MQ=255)
                    tGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCAccgccg                                                  >  1:2579055/1‑70 (MQ=255)
                    tGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCAccgcc                                                   >  1:4293917/1‑69 (MQ=255)
                      aCGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGcgcc                                                               <  1:5137532/55‑1 (MQ=255)
                      aCGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGcgcc                                                               <  1:3223092/55‑1 (MQ=255)
                      aCGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGcgcc                                                               <  1:3152915/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:5501552/1‑55 (MQ=255)
                             tCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCa                                                        >  1:3851941/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)
                                             cAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAg                                           >  1:560931/1‑52 (MQ=255)
                                             cAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTg                         <  1:3944400/70‑1 (MQ=255)
                                              aaaaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAata                                 >  1:915605/1‑61 (MQ=255)
                                              aaaaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAata                                 >  1:5113625/1‑61 (MQ=255)
                                                aaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTc                                             <  1:4460276/47‑1 (MQ=255)
                                                aaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAAc                     <  1:310011/71‑1 (MQ=255)
                                                aaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAAc                     <  1:210974/71‑1 (MQ=255)
                                                aaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAAc                     >  1:2298716/1‑71 (MQ=255)
                                                aaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAAc                     <  1:3913188/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)
                                                      ggCCCTGCTCGCCGGCGCGCGCCGCTTCGACcgccgc                                                 >  1:1004657/2‑37 (MQ=25)
                                                        cccTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAattcc             <  1:2684443/71‑4 (MQ=255)
                                                         ccTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCg            >  1:5478507/1‑71 (MQ=255)
                                                          cTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                                              <  1:74979/36‑1 (MQ=16)
                                                           tGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGt           >  1:5684660/1‑70 (MQ=255)
                                                            gTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGtt          >  1:4153744/1‑70 (MQ=255)
                                                             ttCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATAtt                               >  1:4971986/1‑48 (MQ=255)
                                                             ttCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATAtt                               >  1:1426970/1‑48 (MQ=255)
                                                             ttCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTaa        >  1:3250491/1‑71 (MQ=255)
                                                             ttCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTaa        <  1:3532087/71‑1 (MQ=255)
                                                               cGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCATTGATc      >  1:2974824/1‑71 (MQ=255)
                                                               cGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCATTGATc      >  1:5670102/1‑71 (MQ=255)
                                                                gCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCa     <  1:4874683/71‑1 (MQ=255)
                                                                 ccGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAg    <  1:265479/71‑1 (MQ=255)
                                                                  cGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCt                          >  1:3075071/1‑48 (MQ=38)
                                                                   ggCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCg  <  1:1418521/71‑1 (MQ=255)
                                                                   ggCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCg  <  1:3486042/71‑1 (MQ=255)
                                                                   ggCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCg  <  1:4441120/71‑1 (MQ=255)
                                                                    gCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAatat              >  1:4307586/1‑58 (MQ=255)
                                                                      aCGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCg  >  1:4356910/1‑68 (MQ=255)
                                                                      |                                                                   
TAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCG  >  Exported/2068586‑2068723

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: