Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,218,665	T→C	T138A (ACC→GCC)	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,665	0	T	C	100.0%	200.3 / NA	54	T138A (ACC→GCC)	aer	fused signal transducer for aerotaxis sensory component/methyl accepting chemotaxis component
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/27);  total (27/27)

CCAAAGCATGGCCGCCAGCAAGATAAACATCAGGGTCATCACTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAAT  >  NC_000913/3218534‑3218799                                                                                                                                    |                                                                                                                                       ccAAAGCATGGCCGCCAGCAAAATAAACATCAGGGTCATCACTCCACGCGCCCGCCAGCGAAGGGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCg                                                                                                                                 <  2:407602/139‑1 (MQ=255)   aaaGCATGGCCGCCAGCAAGATAAACATCAGGGTCATCACTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGc                                                                                                                               <  2:647243/139‑1 (MQ=255)     aGCATGGCCGCCAGCAAGATAAACATCAGGGTCATCACTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGt                                                                                                                             >  2:220069/1‑139 (MQ=255)        aTGGCCGCCAGCAAGATAAACATCAGGGTCATCACTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCa                                                                                                                          <  2:958413/139‑1 (MQ=255)        aTGGCCGCCAGCAAGATAAACATCAGGGTCATCACTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCa                                                                                                                          <  1:220069/139‑1 (MQ=255)        aTGGCCGCCAGCAAGATAAACATCAGGGTCATCACTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCa                                                                                                                          <  2:394939/139‑1 (MQ=255)          ggCCGCCAGCAAGATAAACATCAGGGTCATCACTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAAc                                                                                                                        >  1:532025/1‑139 (MQ=255)          ggCCGCCAGCAAGATAAACATCAGGGTCATCACTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAAc                                                                                                                        <  2:532025/139‑1 (MQ=255)           gccgccAGCAAGATAAACATCAGGGTCATCACTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACg                                                                                                                       <  2:1114651/139‑1 (MQ=255)                                 gggTCATCACTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCAccgc                                                                                                 >  2:252927/1‑139 (MQ=255)                                 gggTCATCACTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCAccgc                                                                                                 >  1:470650/1‑139 (MQ=255)                                  ggTCATCACTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCAccgcc                                                                                                >  2:477758/1‑139 (MQ=255)                                     catcaCTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCcgcg                                                                                             <  1:163919/139‑1 (MQ=255)                                       tcaCTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGAt                                                                                           <  1:241236/139‑1 (MQ=255)                                            ccACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCtt                                                                                      >  1:667918/1‑139 (MQ=255)                                              aCGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCa                                                                                    >  1:159238/1‑139 (MQ=255)                                               cgcgCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCAt                                                                                   >  2:611474/1‑139 (MQ=255)                                                   cccgccAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGt                                                                               >  2:23682/1‑139 (MQ=255)                                                     cgccAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCg                                                                             >  2:471504/1‑139 (MQ=255)                                                            gAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGt                                                                      <  1:252927/139‑1 (MQ=255)                                                                          ggCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACatat                                                        <  2:470650/139‑1 (MQ=255)                                                                             aGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCa                                                                                                                          >  1:298068/1‑70 (MQ=255)                                                                             aGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCa                                                                                                                          <  2:298068/70‑1 (MQ=255)                                                                               tttACCCAGCCAGCCTTTACGCACCACCCGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGCATCGACATATAGCCa                                                   <  2:13060/139‑1 (MQ=255)                                                                                  aCCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCcgcg                                                                                             >  2:824777/1‑94 (MQ=255)                                                                                  aCCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCcgcg                                                                                             <  1:824777/94‑1 (MQ=255)                                                                                  aCCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTg                                                >  2:928699/1‑139 (MQ=255)                                                                                        ccagccTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCg                                                                                                                                 <  1:170272/52‑1 (MQ=255)                                                                                        ccagccTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCg                                                                                                                                 >  2:170272/1‑52 (MQ=255)                                                                                                tACGCACCACCAGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTcgcg                                  <  2:428589/139‑1 (MQ=255)                                                                                                      ccaccaGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCa                                                                                                                          <  1:120369/45‑1 (MQ=255)                                                                                                      ccaccaGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCa                                                                                                                          >  2:120369/1‑45 (MQ=255)                                                                                                      ccaccaGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTAc                                                                                                               <  2:991745/56‑1 (MQ=255)                                                                                                      ccaccaGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTAc                                                                                                               >  1:991745/1‑56 (MQ=255)                                                                                                      ccaccaGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCAt                            <  1:928699/139‑1 (MQ=255)                                                                                                       caccaGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACg                                                                                                                       <  1:517528/47‑1 (MQ=255)                                                                                                       caccaGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACg                                                                                                                       >  2:517528/1‑47 (MQ=255)                                                                                                       caccaGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCgcgc                                 <  2:1042740/133‑1 (MQ=255)                                                                                                       caccaGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCgcgc                                 >  1:1042740/1‑133 (MQ=255)                                                                                                          caGGCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGt                        >  2:580400/1‑139 (MQ=255)                                                                                                            ggCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTAc                      >  2:395459/1‑139 (MQ=255)                                                                                                             gCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCa                                                                                                                          <  2:7542/38‑1 (MQ=255)                                                                                                             gCCTTTATGAATACGCTTACTGGCACGTCCGGCGTTCa                                                                                                                          >  1:7542/1‑38 (MQ=255)                                                                                                                  tATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCggg                                                                       <  1:696111/84‑1 (MQ=255)                                                                                                                  tATGAATACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCggg                                                                       >  2:696111/1‑84 (MQ=255)                                                                                                                       aTACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGccc           >  2:1108768/1‑139 (MQ=255)                                                                                                                        tACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCg          <  1:566356/139‑1 (MQ=255)                                                                                                                        tACGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCg          <  1:471504/139‑1 (MQ=255)                                                                                                                         aCGCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACGGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGc         <  1:804482/139‑1 (MQ=255)                                                                                                                           gCTTACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGAGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTGTCCCTAGCGCACCATCGGTACCGCATTGGCCCGCAc       >  2:703163/1‑139 (MQ=255)                                                                                                                              tACTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCa    >  1:870972/1‑139 (MQ=255)                                                                                                                                cTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATc                                                               >  2:646606/1‑78 (MQ=255)                                                                                                                                cTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATc                                                               <  1:646606/78‑1 (MQ=255)                                                                                                                                cTGGCACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCaat  >  2:423438/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                       CCAAAGCATGGCCGCCAGCAAGATAAACATCAGGGTCATCACTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAAT  >  NC_000913/3218534‑3218799

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