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%	168.8 / NA	53	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 (28/25);  total (28/25)

TTCAGTTCATCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCAC  >  NC_000913/3218232‑3218499                                                                                                                                      |                                                                                                                                       ttcagttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGc                                                                                                                                   >  2:381745/1‑139 (MQ=255)  tcagttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATGAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCt                                                                                                                                  >  1:979628/1‑139 (MQ=255)       tcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGAt                                                                                                                             >  1:1326223/1‑139 (MQ=255)         aTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTc                                                                                                                           <  2:979628/139‑1 (MQ=255)              tGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATg                                                                                                                      <  1:381745/139‑1 (MQ=255)                ccTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                                    >  1:1236838/1‑139 (MQ=255)                     cGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAaca                                                                                                               >  2:1090390/1‑139 (MQ=255)                               tcACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGtt                                                                                                     >  1:1269999/1‑139 (MQ=255)                                  cTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTtctc                                                                                                  >  1:151768/1‑139 (MQ=255)                                       aTTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTc                                                                                             <  1:1090390/139‑1 (MQ=255)                                         ttCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGc                                                                                           <  1:133201/139‑1 (MQ=255)                                          tCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGcc                                                                                          <  2:1269999/139‑1 (MQ=255)                                           cTGACACTGGACACCTGGCTTGAGACATCGTTAATT‑‑CCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCAcc                                                                                       <  1:1040750/139‑1 (MQ=255)                                                    gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAaca                                                                                                               >  2:608327/1‑108 (MQ=255)                                                    gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAaca                                                                                                               <  1:608327/108‑1 (MQ=255)                                                    gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTc                                                                                             <  2:806238/126‑1 (MQ=255)                                                    gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTc                                                                                             >  1:806238/1‑126 (MQ=255)                                                    gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTg                                                                                >  1:39587/1‑139 (MQ=255)                                                    gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTg                                                                                >  1:167231/1‑139 (MQ=255)                                                              ttGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGtt                                                                                                     <  2:697831/108‑1 (MQ=255)                                                              ttGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGtt                                                                                                     >  1:697831/1‑108 (MQ=255)                                                               tGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCa                                                                     <  1:1062616/139‑1 (MQ=255)                                                                               aGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTg                                                                            >  2:927536/1‑116 (MQ=255)                                                                               aGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTg                                                                            <  1:927536/116‑1 (MQ=255)                                                                                 ccAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCa                                                                     <  2:1321314/121‑1 (MQ=255)                                                                                 ccAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCa                                                                     >  1:1321314/1‑121 (MQ=255)                                                                                  cAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGca                                                  <  2:568517/139‑1 (MQ=255)                                                                                   aaCGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATggg                                                                       >  2:534583/1‑117 (MQ=255)                                                                                   aaCGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATggg                                                                       <  1:534583/117‑1 (MQ=255)                                                                                     cGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGc                                                                                           >  1:1001721/1‑95 (MQ=255)                                                                                     cGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGc                                                                                           <  2:1001721/95‑1 (MQ=255)                                                                                           aTCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAaca                                                                                                               >  2:1021068/1‑69 (MQ=255)                                                                                           aTCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAaca                                                                                                               <  1:1021068/69‑1 (MQ=255)                                                                                                gCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTAt                                                                                                                                      >  2:760647/1‑41 (MQ=255)                                                                                                     gTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                                    <  2:89322/54‑1 (MQ=255)                                                                                                     gTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACa                               >  1:1073374/1‑139 (MQ=255)                                                                                                     gTTGCCCTACCGCACGTAATGCCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                                    >  1:89322/1‑54 (MQ=255)                                                                                                       tGCCCTGCCGCACGTAATGTCAGCCCCGGCTTATCGCTGCGATTCAGATGCTCAACACTAGTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAgg                             <  1:348223/139‑1 (MQ=255)                                                                                                                gCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCca                    >  2:652300/1‑139 (MQ=255)                                                                                                                      aaTGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATa              >  2:271530/1‑139 (MQ=255)                                                                                                                        tGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATAc             >  2:1151789/1‑138 (MQ=255)                                                                                                                        tGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATAc             <  1:1151789/138‑1 (MQ=255)                                                                                                                            aGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCcaa                   <  1:804824/128‑1 (MQ=255)                                                                                                                            aGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCcaa                   >  2:804824/1‑128 (MQ=255)                                                                                                                              ccccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGgcg                                                    >  2:509221/1‑93 (MQ=255)                                                                                                                              ccccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGgcg                                                    <  1:509221/93‑1 (MQ=255)                                                                                                                              ccccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATAc             <  2:250562/132‑1 (MQ=255)                                                                                                                              ccccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATAc             >  1:250562/1‑132 (MQ=255)                                                                                                                              ccccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAAc      >  1:297002/1‑139 (MQ=255)                                                                                                                                 cAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGca   <  2:151768/139‑1 (MQ=255)                                                                                                                                  aGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGcac  <  1:435522/139‑1 (MQ=255)                                                                                                                                   gCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCg                                                        <  2:1319788/84‑1 (MQ=255)                                                                                                                                   gCTTATCGCTGCGATTCAGATGCTAAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCg                                                        >  1:1319788/1‑84 (MQ=255)                                                                                                                                      |                                                                                                                                       TTCAGTTCATCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTAACGCAC  >  NC_000913/3218232‑3218499

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