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_000913	3,218,365	0	C	T	100.0%	192.7 / NA	59	E238K (GAG→AAG)	aer	fused signal transducer for aerotaxis sensory component/methyl accepting chemotaxis component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/30);  total (29/30)

TCGTTCAGTTCATCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACA  >  NC_000913/3218229‑3218500                                                                                                                                         |                                                                                                                                        tcgttcAGTTCATCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCG‑‑CGTAATGTCAGCCCCAGCTTATCg                                                                                                                                     <  1:788610/139‑1 (MQ=255)    ttcagttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGc                                                                                                                                    <  1:1005149/139‑1 (MQ=255)    ttcagttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGc                                                                                                                                    >  1:687394/1‑139 (MQ=255)            aTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTCGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTc                                                                                                                            >  2:185644/1‑139 (MQ=255)                   ccTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGTCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                                     >  2:183884/1‑139 (MQ=255)                   ccTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                                     <  1:1015156/139‑1 (MQ=255)                   ccTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                                     >  2:413354/1‑139 (MQ=255)                    cTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTc                                                                                                                    <  2:346878/139‑1 (MQ=255)                      ttCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCaa                                                                                                                  >  2:171374/1‑139 (MQ=255)                      ttCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATCTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCaa                                                                                                                  <  1:386410/139‑1 (MQ=255)                             gCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTGTCGCTGCGATTCAGATGCTCAACACTAtt                                                                                                           >  2:681333/1‑139 (MQ=255)                             gCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTAtt                                                                                                           <  1:1017811/139‑1 (MQ=255)                                           tttCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCg                                                                                             <  2:1011517/139‑1 (MQ=255)                                             tCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGcc                                                                                           >  2:22043/1‑139 (MQ=255)                                              cTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCa                                                                                          <  1:222988/139‑1 (MQ=255)                                                  cacTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCtt                                                                                      <  2:906047/139‑1 (MQ=255)                                                    cTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGTTTCAg                                                                                                                          <  1:347837/101‑1 (MQ=255)                                                    cTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGTTTCAg                                                                                                                          >  2:347837/1‑101 (MQ=255)                                                      ggACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGt                                                                                  >  1:331638/1‑139 (MQ=255)                                                       gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTAt                                                                                                                                       <  2:135335/85‑1 (MQ=255)                                                       gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTAt                                                                                                                                       >  1:135335/1‑85 (MQ=255)                                                       gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGc                                                                                            >  1:1038146/1‑128 (MQ=255)                                                       gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGc                                                                                            <  2:1038146/128‑1 (MQ=255)                                                       gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGt                                                                                  >  2:907225/1‑138 (MQ=255)                                                       gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGt                                                                                  <  1:907225/138‑1 (MQ=255)                                                       gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTg                                                                                 >  1:68106/1‑139 (MQ=255)                                                       gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTg                                                                                 >  2:467618/1‑139 (MQ=255)                                                       gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTg                                                                                 >  2:142304/1‑139 (MQ=255)                                                                 ttGAGACATCGTTAATTAGCCAACGGCACATCAGGTCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGc                                                                       <  1:183884/139‑1 (MQ=255)                                                                 ttGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGc                                                                       <  1:413354/139‑1 (MQ=255)                                                                  tGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCa                                                                      <  1:467618/139‑1 (MQ=255)                                                                  tGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCa                                                                      <  1:104610/139‑1 (MQ=255)                                                                   gagaCATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCaa                                                                     >  2:266080/1‑139 (MQ=255)                                                                         tCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCTCTACCGCACGTAATGTCAGCCCCAGCATATCGCTGCGATTCAGATGCTCAGCACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACAtttt                                                               >  1:672298/1‑139 (MQ=255)                                                                         tCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACAtttt                                                               >  1:228035/1‑139 (MQ=255)                                                                               aTTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCg                                                         <  1:545617/139‑1 (MQ=255)                                                                                     cAACGGCACATCAGGCCAAGTCGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGca                                                   <  1:185644/139‑1 (MQ=255)                                                                                           cacaTCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCt                                             >  1:777385/1‑139 (MQ=255)                                                                                            acaTCAGGCCAAGTTGCCCTACCGCATGTAATGTCAGTCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTg                                            <  2:672298/139‑1 (MQ=255)                                                                                            acaTCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTg                                            <  2:228035/139‑1 (MQ=255)                                                                                             caTCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACt                                                                                                              >  1:535115/1‑72 (MQ=255)                                                                                             caTCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACt                                                                                                              <  2:535115/72‑1 (MQ=255)                                                                                              aTCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCa                                                                      <  1:78381/111‑1 (MQ=255)                                                                                              aTCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCa                                                                      >  2:78381/1‑111 (MQ=255)                                                                                                   gCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCa                                     >  2:493651/1‑139 (MQ=255)                                                                                                      aaGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCaa                                                                                                                  >  1:623191/1‑59 (MQ=255)                                                                                                      aaGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCaa                                                                                                                  <  2:623191/59‑1 (MQ=255)                                                                                                          tGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAgg                              <  1:142304/139‑1 (MQ=255)                                                                                                             ccTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGgcgc                           >  1:940558/1‑139 (MQ=255)                                                                                                             ccTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGgcgc                           <  1:681333/139‑1 (MQ=255)                                                                                                                 ccGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGc                       <  1:22043/139‑1 (MQ=255)                                                                                                                        tAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAAt                <  2:68106/139‑1 (MQ=255)                                                                                                                           tGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATAcc             >  1:389711/1‑139 (MQ=255)                                                                                                                                 ccccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAAc       >  2:341223/1‑139 (MQ=255)                                                                                                                                 ccccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAAc       >  2:125436/1‑139 (MQ=255)                                                                                                                                     aGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGgcg                            >  1:867805/1‑114 (MQ=255)                                                                                                                                     aGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGgcg                            <  2:867805/114‑1 (MQ=255)                                                                                                                                     aGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGcac   >  1:190958/1‑139 (MQ=255)                                                                                                                                      gCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGcaca  <  1:493651/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        TCGTTCAGTTCATCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACA  >  NC_000913/3218229‑3218500

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