Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,218,365 C→T E238K (GAG→AAG)  aer ← fused signal transducer for aerotaxis sensory component/methyl accepting chemotaxis component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,218,3650CT98.4% 195.1 / NA 63E238K (GAG→AAG) aerfused signal transducer for aerotaxis sensory component/methyl accepting chemotaxis component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (29/33);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TTCAGTTCATCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAA  >  NC_000913/3218232‑3218494
                                                                                                                                     |                                                                                                                                 
ttcagttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGc                                                                                                                              <  1:526754/139‑1 (MQ=255)
     ttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGa                                                                                                                         <  2:1117548/139‑1 (MQ=255)
             tGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATg                                                                                                                 >  2:378949/1‑139 (MQ=255)
               ccTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                               >  2:123421/1‑139 (MQ=255)
                  ttCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCaa                                                                                                            <  1:572086/139‑1 (MQ=255)
                          cGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTa                                                                                                    >  2:97038/1‑139 (MQ=255)
                                   gCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCg                                                                                           <  2:41266/139‑1 (MQ=255)
                                      aTTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTc                                                                                        <  1:1478344/139‑1 (MQ=255)
                                      aTTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTc                                                                                        <  1:216678/139‑1 (MQ=255)
                                      aTTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTc                                                                                        <  1:378949/139‑1 (MQ=255)
                                          cTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCa                                                                                    <  1:123421/139‑1 (MQ=255)
                                               acTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTc                                                                               <  1:97038/139‑1 (MQ=255)
                                                cTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCa                                                                              <  1:1252302/139‑1 (MQ=255)
                                                 tGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTtct                                                                                              >  1:595971/1‑122 (MQ=255)
                                                 tGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTtct                                                                                              <  2:595971/122‑1 (MQ=255)
                                                 tGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCt                                                                                 <  1:918146/135‑1 (MQ=255)
                                                 tGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCt                                                                                 >  2:918146/1‑135 (MQ=255)
                                                   gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCaa                                                                                                            <  2:1177904/106‑1 (MQ=255)
                                                   gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACTTTCTCCGGTCGCCACCTTCAGt                                                                            <  1:1224290/138‑1 (MQ=255)
                                                   gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGt                                                                            >  2:1224290/1‑138 (MQ=255)
                                                   gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTCATCGCTGCGATTCAGATGCTCaa                                                                                                            >  1:1177904/1‑106 (MQ=255)
                                                       ccTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTg                                                                       <  2:1351551/139‑1 (MQ=255)
                                                          ggCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGc                                                                                                                           <  1:437694/84‑1 (MQ=255)
                                                          ggCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGc                                                                                                                           >  2:437694/1‑84 (MQ=255)
                                                           gCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATgg                                                                   >  1:1254622/1‑139 (MQ=255)
                                                            cTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATggg                                                                  <  2:268385/139‑1 (MQ=255)
                                                              tGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCGCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATgggta                                                                <  2:1223450/139‑3 (MQ=255)
                                                                      cGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTc                                                        >  2:735694/1‑139 (MQ=255)
                                                                      cGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTc                                                        >  2:632833/1‑139 (MQ=255)
                                                                          aaTTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAg                                                                                                                    <  1:814254/75‑1 (MQ=255)
                                                                          aaTTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAg                                                                                                                    >  2:814254/1‑75 (MQ=255)
                                                                              aGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAaca                                                                                                          <  2:1169951/81‑1 (MQ=255)
                                                                              aGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAaca                                                                                                          >  1:1169951/1‑81 (MQ=255)
                                                                               gCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGgcg                                               <  1:735694/139‑1 (MQ=255)
                                                                                ccAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGgcg                                               >  2:628596/1‑138 (MQ=255)
                                                                                ccAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGgcg                                               <  1:628596/138‑1 (MQ=255)
                                                                                    cGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGTCTg                                                                       >  2:750815/1‑110 (MQ=255)
                                                                                    cGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGTCTg                                                                       <  1:750815/110‑1 (MQ=255)
                                                                                        acaTCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTg                                      <  1:1197298/139‑1 (MQ=255)
                                                                                          aTCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGcc                                    >  2:622352/1‑139 (MQ=255)
                                                                                          aTCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGcc                                    >  2:52056/1‑139 (MQ=255)
                                                                                             aGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCAtt                                 <  1:632833/139‑1 (MQ=255)
                                                                                                    gTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATg                                                                                                                 >  2:1017328/1‑52 (MQ=255)
                                                                                                    gTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATg                                                                                                                 <  1:1017328/52‑1 (MQ=255)
                                                                                                    gTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGc                                     <  1:754601/128‑1 (MQ=255)
                                                                                                    gTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGc                                     >  2:754601/1‑128 (MQ=255)
                                                                                                    gTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACa                          >  2:1314421/1‑139 (MQ=255)
                                                                                                              cGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTc                                                                                        >  2:957274/1‑67 (MQ=255)
                                                                                                              cGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTc                                                                                        <  1:957274/67‑1 (MQ=255)
                                                                                                                cACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCcaa              <  2:1254622/139‑1 (MQ=255)
                                                                                                                cACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCcaa              >  1:1459496/1‑139 (MQ=255)
                                                                                                                  cGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCcaaca            >  1:916472/1‑139 (MQ=255)
                                                                                                                  cGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCcaaca            >  2:115025/1‑139 (MQ=255)
                                                                                                                       tGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATAcc       >  2:1369830/1‑139 (MQ=255)
                                                                                                                           aGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGgcg                                               <  2:1178870/95‑1 (MQ=255)
                                                                                                                           aGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGgcg                                               >  1:1178870/1‑95 (MQ=255)
                                                                                                                            gCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTaa  >  1:418182/1‑139 (MQ=255)
                                                                                                                             ccccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTAt                                                     >  2:46611/1‑87 (MQ=255)
                                                                                                                             ccccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTAt                                                     <  1:46611/87‑1 (MQ=255)
                                                                                                                              cccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTc                                                                               <  1:791688/60‑1 (MQ=255)
                                                                                                                              cccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTc                                                                               >  2:791688/1‑60 (MQ=255)
                                                                                                                                 aGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTAt                                                     >  2:216493/1‑83 (MQ=255)
                                                                                                                                 aGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTAt                                                     <  1:216493/83‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                 
TTCAGTTCATCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAA  >  NC_000913/3218232‑3218494

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: