Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,068,647 T→C 20.6% E381E (GAA→GAG 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,6470TC20.6% 121.1 / 26.9 63E381E (GAA→GAGaerfused signal transducer for aerotaxis sensory component and methyl accepting chemotaxis component
Reads supporting (aligned to +/- strand):  ref base T (27/23);  new base C (7/6);  total (34/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 = 8.76e-01

GGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAA  >  Exported/2068578‑2068717
                                                                     |                                                                      
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)
                cGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGc                                                                    >  1:2561294/1‑58 (MQ=255)
                   ctggctAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCGTc                                                    >  1:5514198/1‑71 (MQ=255)
                       ctAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc                                                                      <  1:5691084/49‑1 (MQ=255)
                         aaaTGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc                                                                      <  1:96151/47‑1 (MQ=255)
                           aTGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGt                                                                        >  1:591358/1‑43 (MQ=255)
                           aTGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGc                                                                    >  1:3498637/1‑47 (MQ=255)
                            tGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCAccgccg                                            >  1:2579055/1‑70 (MQ=255)
                            tGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCAccgcc                                             >  1:4293917/1‑69 (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)
                                                 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:210974/71‑1 (MQ=255)
                                                        aaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAAc               <  1:310011/71‑1 (MQ=255)
                                                        aaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAAc               >  1:2298716/1‑71 (MQ=255)
                                                          ccTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACg              >  1:2514079/1‑70 (MQ=255)
                                                          ccTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGc             <  1:1503705/71‑1 (MQ=255)
                                                          ccTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGc             <  1:61465/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)
                                                                     |                                                                      
GGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAA  >  Exported/2068578‑2068717

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: