Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,573,558	G→A	10.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_002947	2,573,558	0	G	A	10.2%	133.3 / 4.9	49	G105D (GGT→GAT)	PP_2257	aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (24/20);  new base A (3/2);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

CCACCCGGACATGCCGCCGGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCC  >  NC_002947/2573415‑2573700                                                                                                                                                |                                                                                                                                               ccACCCGGACATGCCGCCGGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGc                                                                                                                                           >  2:85770/1‑149 (MQ=255)   aCCCGGACATGCCGCCGGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCgg                                                                                                                                         <  1:85770/149‑1 (MQ=255)                        tgtACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTg                                                                                                                    >  4:202531/1‑149 (MQ=255)                            cGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTgcgc                                                                                                                 <  4:152909/148‑1 (MQ=255)                            cGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTgcgc                                                                                                                 <  6:55256/148‑1 (MQ=255)                             gAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGcgcg                                                                                                                >  6:235114/1‑148 (MQ=255)                               gTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTAt                                                                                                             >  5:284692/1‑149 (MQ=255)                                     gtgGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTa                                                                                                       >  8:159202/1‑149 (MQ=255)                                      tgGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTAc                                                                                                      >  1:347923/1‑149 (MQ=255)                                              ccTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGa                                                                                                <  5:2277/147‑1 (MQ=255)                                              acTGAAAGCCGGCAAGAGCTGGATGGGTATCGACAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCAACGTCACGCCCATCCTCGCAGGTGGGCGGGGAGTGGGTTACAAGTCGGTGCGCGTATGCACTACGCGCGAg                                                                                               >  6:127600/2‑148 (MQ=255)                                                tGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCgg                                                                                                                                         >  1:316717/1‑104 (MQ=255)                                                tGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCgg                                                                                                                                         <  2:316717/104‑1 (MQ=255)                                                 gAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAg                                                                                           >  7:47700/1‑149 (MQ=255)                                                  aacgCCGGCAAGAGCTGGTTGGGTATCGTCAAGAACGGCTGCCGTAACTGCGACTACTACTTGGTCAATGCCTTCGTCACGCCCATCCTCGAATGGGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGATCaa                                                                                            <  5:127600/144‑1 (MQ=255)                                                  aaaGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCaa                                                                                            <  3:202531/147‑1 (MQ=255)                                                  aaaGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAg                                                                                           <  6:284692/148‑1 (MQ=255)                                                   aaGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTa                                                                                         <  8:47700/149‑1 (MQ=255)                                                      ccGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGcc                                                                                      >  1:254961/1‑149 (MQ=255)                                                        ggCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCg                                                                                    >  4:324674/1‑149 (MQ=255)                                                                  ggATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGt                                                                                          <  3:74987/133‑1 (MQ=255)                                                                  ggATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGt                                                                                          >  4:74987/1‑133 (MQ=255)                                                                   gATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGct                                                                          <  2:254961/148‑2 (MQ=255)                                                                    aTGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCt                                                                        <  5:235114/149‑1 (MQ=255)                                                                             gTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGcgcg                                                               >  4:161080/1‑149 (MQ=255)                                                                              tCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTActa                                                                                                                           <  4:198585/88‑3 (MQ=255)                                                                              tCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGa                                                                                                                           >  3:198585/1‑88 (MQ=255)                                                                                       gCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAAc                                                     <  5:341253/149‑1 (MQ=255)                                                                                        cTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACg                                                    <  7:159202/149‑1 (MQ=255)                                                                                        cTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACg                                                    >  7:24547/1‑149 (MQ=255)                                                                                        cTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCAGCCGCGCata                                                     >  1:147315/1‑146 (MQ=255)                                                                                           ccGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCAGCCGCGCataca                                                   <  2:147315/147‑5 (MQ=255)                                                                                                        tactactGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTc                                    >  8:48679/1‑149 (MQ=255)                                                                                                          ctactGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTc                                    <  3:161080/147‑1 (MQ=255)                                                                                                                     aTGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAACAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCgg                       >  2:251299/1‑149 (MQ=255)                                                                                                                         cTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCt                     >  5:165071/1‑147 (MQ=255)                                                                                                                              acccGCCCACCTTCGAGGATGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTg              >  7:328646/4‑149 (MQ=255)                                                                                                                                ccGCCCACCTTCGAGGATGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGgcgccgc                            >  2:58628/2‑133 (MQ=255)                                                                                                                                ccGCCCACCTTCGAGGATGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGgcgccgc                            <  1:58628/132‑1 (MQ=255)                                                                                                                                 cGCCCACCTTCGAGGATGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGa                   >  8:47764/1‑141 (MQ=255)                                                                                                                                 cGCCCACCTTCGAGGATGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGa                   <  7:47764/141‑1 (MQ=255)                                                                                                                                      aTCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTg              >  2:242877/1‑141 (MQ=255)                                                                                                                                      aTCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTg              <  1:242877/141‑1 (MQ=255)                                                                                                                                          tCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGa                   <  2:215157/132‑1 (MQ=255)                                                                                                                                          tCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGgc         >  2:110664/1‑142 (MQ=255)                                                                                                                                          tCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGgc         <  1:110664/142‑1 (MQ=255)                                                                                                                                          tCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAggcc  >  3:34234/1‑149 (MQ=255)                                                                                                                                          tCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAggcc  >  5:115121/1‑149 (MQ=255)                                                                                                                                          tCGAAGGTGGGCGGGTAGTGGGTTACCAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGa                   >  1:215157/1‑132 (MQ=255)                                                                                                                                                |                                                                                                                                               CCACCCGGACATGCCGCCGGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCC  >  NC_002947/2573415‑2573700

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