Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,068,671 C→T 16.1% A373A (GCG→GCA 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,6710CT16.1% 143.5 / 17.0 62A373A (GCG→GCAaerfused signal transducer for aerotaxis sensory component and methyl accepting chemotaxis component
Reads supporting (aligned to +/- strand):  ref base C (28/24);  new base T (5/5);  total (33/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.67e-01

TGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCAATGC  >  Exported/2068606‑2068741
                                                                 |                                                                      
tGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCAccgccg                                                                    >  1:2579055/1‑70 (MQ=255)
tGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCAccgcc                                                                     >  1:4293917/1‑69 (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)
          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)
                 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)
                      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: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:3532087/71‑1 (MQ=255)
                                         ttCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTaa                          >  1:3250491/1‑71 (MQ=255)
                                           cGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCATTGATc                        >  1:5670102/1‑71 (MQ=255)
                                           cGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCATTGATc                        >  1:2974824/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:3486042/71‑1 (MQ=255)
                                               ggCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCg                    <  1:4441120/71‑1 (MQ=255)
                                               ggCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCg                    <  1:1418521/71‑1 (MQ=255)
                                                gCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAatat                                >  1:4307586/1‑58 (MQ=255)
                                                  aCGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCg                    >  1:4356910/1‑68 (MQ=255)
                                                      gccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAAt             <  1:5300292/71‑1 (MQ=255)
                                                       ccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGtt                            <  1:5531113/55‑1 (MQ=255)
                                                       ccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGtt                            <  1:257843/55‑1 (MQ=255)
                                                       ccgcTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATg            >  1:3617225/1‑71 (MQ=255)
                                                            tCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCaa                                   >  1:2266063/1‑43 (MQ=255)
                                                            tCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATc                        >  1:40211/1‑54 (MQ=255)
                                                            tCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATgg           >  1:4345646/1‑67 (MQ=255)
                                                              cACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCaa     >  1:994310/1‑71 (MQ=255)
                                                              cACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCaa     <  1:810295/71‑1 (MQ=255)
                                                              cACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCaa     <  1:1210731/71‑1 (MQ=255)
                                                               aCCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGt                <  1:2488465/59‑1 (MQ=255)
                                                                 cgccgcATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCAATgc  <  1:5189686/71‑1 (MQ=255)
                                                                 |                                                                      
TGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCAATGC  >  Exported/2068606‑2068741

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: