Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,573,549	T→C	11.1%	I102T (ATC→ACC)	PP_2257 →	aerotaxis receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,573,549	0	T	C	11.1%	112.7 / 5.7	45	I102T (ATC→ACC)	PP_2257	aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (20/20);  new base C (3/2);  total (23/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 = 6.38e-01

CTGGTGCGCCACCCGGACATGCCGCCGGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTG  >  NC_002947/2573407‑2573688                                                                                                                                               |                                                                                                                                            cTGGTGCGCCACCCGGACATGCCGCCGGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGa                                                                                                                                       <  2:127356/149‑1 (MQ=255)    gTGCGCCACCCGGACATGCCGCCGGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAg                                                                                                                                     <  6:102037/148‑1 (MQ=255)         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                                                                                                                             <  2:316717/104‑1 (MQ=255)                                                        tGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCgg                                                                                                                             >  1:316717/1‑104 (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       <  7:47764/141‑1 (MQ=255)                                                                                                                                         cGCCCACCTTCGAGGATGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGa       >  8:47764/1‑141 (MQ=255)                                                                                                                                              aTCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTg  >  2:242877/1‑141 (MQ=255)                                                                                                                                              aTCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTg  <  1:242877/141‑1 (MQ=255)                                                                                                                                               |                                                                                                                                            CTGGTGCGCCACCCGGACATGCCGCCGGCGGTGTACGAGTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTG  >  NC_002947/2573407‑2573688

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: ‑