Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,134,986	T→A	12.5%	intergenic (‑47/+18)	PP_4520 ← / ← PP_4521	hypothetical protein/aerotaxis receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,134,986	0	T	A	12.5%	99.3 / 6.4	40	intergenic (‑47/+18)	PP_4520/PP_4521	hypothetical protein/aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (20/15);  new base A (3/2);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

GCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGCAATTCTCACGATCTG  >  NC_002947/5134843‑5135132                                                                                                                                                |                                                                                                                                                   gcagGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCg                                                                                                                                               >  5:50015/1‑149 (MQ=255)       tCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCAAAGCGGGcccc                                                                                                                                         <  4:51996/149‑1 (MQ=255)          cacGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCcgc                                                                                                                                       >  1:59971/1‑148 (MQ=255)              ggTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGAttt                                                                                                                                  <  4:268361/149‑1 (MQ=255)                   ggCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGAGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCg                                                                                                                             >  7:87588/1‑149 (MQ=255)                                      cAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAgg                                                                                                          >  1:143854/1‑149 (MQ=255)                                        gTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGc                                                                                                         >  8:194536/1‑148 (MQ=255)                                           ggctaaGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGGCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTCCCACATACCCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATg                                                                                                     <  2:297509/144‑1 (MQ=255)                                             cATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCggg                                                                                                                                             <  3:227518/107‑1 (MQ=255)                                             cATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCggg                                                                                                                                             >  4:227518/1‑107 (MQ=255)                                               tAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAgg                                                                                                 >  3:124010/1‑149 (MQ=255)                                                aGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCt                                                                                                        >  2:14367/1‑141 (MQ=255)                                                aGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCt                                                                                                        <  1:14367/141‑1 (MQ=255)                                                aGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAgg                                                                                                 >  3:162392/1‑148 (MQ=255)                                                     aTGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTcc                                                                                           >  1:31963/1‑149 (MQ=255)                                                         tCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCa                                                                                       >  5:94222/1‑149 (MQ=255)                                                         cgctgCGAGCGCTGCCAGCGCCAGCAGCGGGGTCTGGGCATGGCGGTCTCCTTCGCCGGGACAGTTTGCCACATAACCCGCGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGATGAAc                                                                                                                      >  7:193970/5‑118 (MQ=255)                                                                       gCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGAGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTc                                                                          <  8:87588/148‑1 (MQ=255)                                                                         aGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTg                                                                       <  7:194536/149‑1 (MQ=255)                                                                          gcgcTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGACAGGCTATGc                                                                                                    >  7:290421/1‑119 (MQ=255)                                                                           cgcTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGc                                                                     >  7:123392/1‑149 (MQ=255)                                                                              tCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGa                                                                  >  4:75238/1‑149 (MQ=255)                                                                                gCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGAt                                                                 <  2:59971/148‑1 (MQ=255)                                                                                      gggTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCcgcg                                                          <  2:31963/149‑1 (MQ=255)                                                                                            gggCATGGCGGTCTCCTTGGCGGGGACAGTTTGCAAAATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGc                                                    <  4:26365/149‑1 (MQ=255)                                                                                             ggCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGc                                                    <  3:75238/148‑1 (MQ=255)                                                                                               cATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGAt                                                                 >  8:82955/1‑133 (MQ=255)                                                                                               cATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGAt                                                                 <  7:82955/133‑1 (MQ=255)                                                                                                 tGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTgg                                                >  6:86489/1‑148 (MQ=255)                                                                                                 tGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTgg                                                <  5:86489/148‑1 (MQ=255)                                                                                                           cTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCa                                     >  6:4866/1‑149 (MQ=255)                                                                                                           cTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCa                                     <  8:123392/149‑1 (MQ=255)                                                                                                                   ggACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCCCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTc                             <  6:94222/149‑1 (MQ=255)                                                                                                                         tttGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGCTTTCCGCGCCTTGCCTGGCTTGCCCGGCGCTCTGGTCGcaccg                       >  8:16622/1‑147 (MQ=255)                                                                                                                         tttGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAg                       >  4:75766/1‑149 (MQ=255)                                                                                                                                      cGGGGCCCGCAAAGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGcc                                          <  1:299236/117‑1 (MQ=255)                                                                                                                                      cGGGGCCCGCAAAGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGcc                                          >  2:299236/1‑117 (MQ=255)                                                                                                                                      cGGGGCCCGCAAAGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGCAATtct           >  8:266494/1‑148 (MQ=255)                                                                                                                                      cGGGGCCCGCAAAGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGCAATtct           >  7:29072/1‑148 (MQ=255)                                                                                                                                              gCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGCAATTCTCACGATctg  <  4:162392/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   GCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGCAATTCTCACGATCTG  >  NC_002947/5134843‑5135132

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