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%	129.2 / NA	40	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 (21/19);  total (21/19)

TCGTTCAGTTCATCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACA  >  NC_000913/3218229‑3218500                                                                                                                                         |                                                                                                                                        tcgttcAGTTCATCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTAt                                                                                                                                       <  2:717419/139‑1 (MQ=255)    ttcagttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGc                                                                                                                                    >  2:283945/1‑139 (MQ=255)    ttcagttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGc                                                                                                                                    >  2:392226/1‑139 (MQ=255)    ttcagttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGc                                                                                                                                    >  1:106964/1‑139 (MQ=255)    ttcagttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACC‑CACGTAATGTCAGCCCCAGCTTATCGCt                                                                                                                                   >  1:181175/1‑139 (MQ=255)      cagttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTg                                                                                                                                  >  2:523107/1‑139 (MQ=255)              cGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAg                                                                                                                          <  1:283945/139‑1 (MQ=255)               ggTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGa                                                                                                                         <  1:198375/139‑1 (MQ=255)                ttgacTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGG‑CACCTGGCTTGAGACATCGTTAATTAGCCAACGGCATATCAGGCCAAGTTGCCGTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATg                                                                                                                       <  1:74061/135‑1 (MQ=255)                 tGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGACTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATg                                                                                                                       <  2:428563/139‑1 (MQ=255)                   ccTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                                     >  1:220012/1‑139 (MQ=255)                   ccTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                                     >  1:27012/1‑139 (MQ=255)                   ccTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                                     <  1:336945/139‑1 (MQ=255)                   ccTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                                     <  1:92368/139‑1 (MQ=255)                   ccTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                                     >  2:540303/1‑139 (MQ=255)                                              cTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCa                                                                                          <  2:27012/139‑1 (MQ=255)                                                  cacTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCtt                                                                                      <  2:220012/139‑1 (MQ=255)                                                       gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTg                                                                                 >  1:497244/1‑139 (MQ=255)                                                        acacCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCCCCTTCAGTGc                                                                                >  1:671852/1‑139 (MQ=255)                                                               gCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCtt                                                                                      <  1:675997/126‑1 (MQ=255)                                                               gCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCtt                                                                                      >  2:675997/1‑126 (MQ=255)                                                                           gTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCaa                                                                                                                  >  1:633192/1‑86 (MQ=255)                                                                           gTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCaa                                                                                                                  >  1:635393/1‑86 (MQ=255)                                                                           gTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCaa                                                                                                                  <  2:633192/86‑1 (MQ=255)                                                                           gTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCaa                                                                                                                  <  2:635393/86‑1 (MQ=255)                                                                              aaTTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATc                                                          <  1:523107/139‑1 (MQ=255)                                                                                    ccAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGgcgc                                                    >  1:70607/1‑139 (MQ=255)                                                                                    ccAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGgcgc                                                    >  2:227937/1‑139 (MQ=255)                                                                                                   gCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCa                                     >  2:1117/1‑139 (MQ=255)                                                                                                                 ccGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGc                       >  2:370349/1‑139 (MQ=255)                                                                                                                    cACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCcaa                    >  1:408841/1‑139 (MQ=255)                                                                                                                        tAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAAt                <  1:1117/139‑1 (MQ=255)                                                                                                                           tGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGgcg                                                     <  2:254115/99‑1 (MQ=255)                                                                                                                           tGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGgcg                                                     >  1:254115/1‑99 (MQ=255)                                                                                                                                 ccccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGTGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAAc       >  1:589030/1‑139 (MQ=255)                                                                                                                                   ccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGc     <  2:408841/139‑1 (MQ=255)                                                                                                                                    cAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGCCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGca    <  2:497244/139‑1 (MQ=255)                                                                                                                                     aGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGcac   >  1:335490/1‑139 (MQ=255)                                                                                                                                      gCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGcaca  <  1:556708/139‑1 (MQ=255)                                                                                                                                      gCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGcaca  <  1:392226/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        TCGTTCAGTTCATCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCACA  >  NC_000913/3218229‑3218500

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