Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,573,558 G→A 7.2% G105D (GGT→GAT)  PP_2257 → aerotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,573,5580GA 7.2% 156.3 / 2.1 55G105D (GGT→GAT) PP_2257aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (28/23);  new base A (2/2);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-01

CCACCCGGACATGCCGCCGGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCC  >  NC_002947/2573415‑2573706
                                                                                                                                               |                                                                                                                                                    
ccACCCGGACATGCCGCCGGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCACCTTCGAGGATGGGc                                                                                                                                                 <  5:206195/149‑1 (MQ=255)
 cACCCGGACATGCCGCCGGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCg                                                                                                                                                >  1:281988/1‑149 (MQ=255)
  aCCCGGACATGCCGCCGGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCgg                                                                                                                                               <  2:281988/149‑1 (MQ=255)
        aCATGCCGCCGGCGGTGTACGCGTTGATGGGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAACAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTg                                                                                                                                         <  6:179916/149‑1 (MQ=255)
          aTGCCGCCGGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTggg                                                                                                                                       >  4:182630/1‑149 (MQ=255)
              cgccgGCGGGGAACGAGTTGATGTGGCCCTACCTGAAACCCGCCAAGCGCTGGCTGGGCATCGTCACGAACCGCTGCCGTAACCGCGACTACTACTGGGTCACTGCCTACGTCACGCCCATCCTCGAAGGTCGGCGGGTAGTGGGTTAc                                                                                                                                   <  8:286448/149‑1 (MQ=255)
               gccgGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAAACGCTGACGTAACGGCGACTACGACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTgggggac                                                                                                                                   >  6:17918/1‑144 (MQ=255)
                ccgGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGa                                                                                                                                 <  7:214627/149‑1 (MQ=255)
                       tgtACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGt                                                                                                                                          <  1:22736/133‑1 (MQ=255)
                       tgtACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGt                                                                                                                                          >  2:22736/1‑133 (MQ=255)
                            gAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGcgcg                                                                                                                      >  4:50442/1‑148 (MQ=255)
                                  aTGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGccc                                                                                                               >  3:201698/1‑149 (MQ=255)
                                  aTGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGcc                                                                                                                <  3:182630/148‑1 (MQ=255)
                                     tgGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTAc                                                                                                            <  6:102190/149‑1 (MQ=255)
                                     tgGCGATACCTGAAAGCCGGCAAGAGCTGGATGGGTAGCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTAc                                                                                                            <  4:176860/149‑1 (MQ=255)
                                            aCCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTAcgcgatat                                                                                                     >  1:243623/1‑145 (MQ=255)
                                            aCCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTAcgc                                                                                                          <  2:243623/144‑1 (MQ=255)
                                             ccttcaaGCCGGCAAGAGCTGGAGGGGTATCGGCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACTCCCATCCTCGAAGGTGGGCGTGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGa                                                                                                      <  5:17918/142‑1 (MQ=255)
                                               tGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGc                                                                                                    >  7:206097/1‑147 (MQ=255)
                                               tGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGc                                                                                                    >  8:45468/1‑147 (MQ=255)
                                               tGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGc                                                                                                    >  2:169415/1‑147 (MQ=255)
                                                gAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCATCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAg                                                                                                 >  7:46669/1‑149 (MQ=255)
                                                 aaaGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTg                                                                                                                                                    >  3:94590/1‑97 (MQ=255)
                                                 aaaGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTg                                                                                                                                                    <  4:94590/97‑1 (MQ=255)
                                                 aaaGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAg                                                                                                 >  1:164501/1‑148 (MQ=255)
                                                 aaaGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGAGAGCAAg                                                                                                 >  8:159555/1‑148 (MQ=255)
                                                          aaGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGc                                                                                       <  3:50442/149‑1 (MQ=255)
                                                              gCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGaa                                                                                   <  3:158972/149‑1 (MQ=255)
                                                                 ggATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCAACGCGCGAGCAAGTAGCCCGTGCCGAAGcc                                                                                >  3:247545/1‑149 (MQ=255)
                                                                       gTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTAt                                                                          >  8:160311/1‑149 (MQ=255)
                                                                             tCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCgg                                                                    >  7:146726/1‑149 (MQ=255)
                                                                             tCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCgg                                                                    >  5:118081/1‑149 (MQ=255)
                                                                               aaGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCt                                                                  >  8:196309/1‑149 (MQ=255)
                                                                                aGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTg                                                                 <  7:160311/149‑1 (MQ=255)
                                                                                  aaCCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTgcg                                                               <  8:46669/149‑1 (MQ=255)
                                                                                            gTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGc                                                                                       <  8:267341/115‑1 (MQ=255)
                                                                                            gTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGc                                                                                       >  7:267341/1‑115 (MQ=255)
                                                                                               aCGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAggcg                                                  <  4:201698/149‑1 (MQ=255)
                                                                                                         ctactGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCgg                                        >  3:238319/1‑149 (MQ=255)
                                                                                                          tactGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGgc                                       >  6:126324/1‑149 (MQ=255)
                                                                                                            ctGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGgcgc                                     <  8:206097/149‑1 (MQ=255)
                                                                                                            ctGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGgcg                                      <  2:164501/148‑1 (MQ=255)
                                                                                                                          aCGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGc                       >  6:78013/1‑149 (MQ=255)
                                                                                                                            gTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCtt                     >  1:206880/1‑149 (MQ=255)
                                                                                                                                 gCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTgg                >  4:170081/1‑149 (MQ=255)
                                                                                                                                 gCCCACCTTCGAGGATGGGCGGGTACTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCATGTAGCCCGTGCAGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCAGAGCTTGCtcg                >  8:165215/1‑147 (MQ=255)
                                                                                                                                    caccttCGAGGATGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAggcggc                                                <  8:92073/109‑1 (MQ=255)
                                                                                                                                    caccttCGAGGATGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAggcggc                                                >  7:92073/6‑114 (MQ=255)
                                                                                                                                     aTCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTg                 <  7:33106/144‑1 (MQ=255)
                                                                                                                                     aTCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTg                 >  8:33106/1‑144 (MQ=255)
                                                                                                                                      tCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAg           <  7:45468/149‑1 (MQ=255)
                                                                                                                                         tCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCAGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAggcc        >  5:53134/1‑149 (MQ=255)
                                                                                                                                         tCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGACCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTCGCGCAggcc        >  6:182202/1‑149 (MQ=255)
                                                                                                                                          cGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAggcct       <  7:196309/149‑1 (MQ=255)
                                                                                                                                            aaGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAggccttg     <  5:126324/149‑1 (MQ=255)
                                                                                                                                               gTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAggccttggcc  <  1:169415/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CCACCCGGACATGCCGCCGGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCC  >  NC_002947/2573415‑2573706

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑