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%	189.1 / NA	58	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 (27/31);  total (27/31)

TTCAGTTCATCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTA  >  NC_000913/3218232‑3218493                                                                                                                                      |                                                                                                                                 ttcagttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGc                                                                                                                             <  1:148894/139‑1 (MQ=255)   cagttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTg                                                                                                                           <  1:949265/139‑1 (MQ=255)   cagttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTg                                                                                                                           <  1:871832/139‑1 (MQ=255)   cagttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTg                                                                                                                           <  1:1445723/139‑1 (MQ=255)      ttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCTATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCg                                                                                                                              <  1:1285173/133‑1 (MQ=255)      ttcaTCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCg                                                                                                                              >  2:1285173/1‑133 (MQ=255)           cGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAg                                                                                                                   >  1:1019278/1‑139 (MQ=255)            ggTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGa                                                                                                                  <  1:1310843/139‑1 (MQ=255)                ccTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                              <  1:215137/139‑1 (MQ=255)                           cGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTa                                                                                                   >  2:973229/1‑139 (MQ=255)                           cGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTa                                                                                                   >  2:74086/1‑139 (MQ=255)                             tctcACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACg                                                                                                 >  2:259979/1‑139 (MQ=255)                               tcACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCaa                                                                                                           >  1:603662/1‑127 (MQ=255)                               tcACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCaa                                                                                                           <  2:603662/127‑1 (MQ=255)                                        tttCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCaa                                                                                                           <  2:555106/118‑1 (MQ=255)                                        tttCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCaa                                                                                                           >  1:555106/1‑118 (MQ=255)                                         ttCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGc                                                                                     <  1:973229/139‑1 (MQ=255)                                          tCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGcc                                                                                    >  1:947065/1‑139 (MQ=255)                                                    gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCg                                                                                                                              >  2:944345/1‑87 (MQ=255)                                                    gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTg                                                                          >  2:566731/1‑139 (MQ=255)                                                    gACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATCTCAGCCCCAGCTTATCg                                                                                                                              <  1:944345/87‑1 (MQ=255)                                                               tGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCa                                                               <  1:74086/139‑1 (MQ=255)                                                                 agaCATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAAc                                                             <  2:1019278/139‑1 (MQ=255)                                                                      tCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                              >  1:840758/1‑85 (MQ=255)                                                                      tCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                              <  2:840758/85‑1 (MQ=255)                                                                      tCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACAtttt                                                        >  2:1335720/1‑139 (MQ=255)                                                                      tCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACAtttt                                                        >  2:1429670/1‑139 (MQ=255)                                                                       cGTTAATTAGCCAACGGCACATCAGGGCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAaca                                                                                                         <  1:897130/89‑1 (MQ=255)                                                                       cGTTAATTAGCCAACGGCACATCAGGCAAAGTGGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAaca                                                                                                         >  2:897130/1‑89 (MQ=255)                                                                           aaTTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCa                                                                                                            >  1:1235414/1‑82 (MQ=255)                                                                           aaTTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCa                                                                                                            <  2:1235414/82‑1 (MQ=255)                                                                                    aCGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGcaca                                          <  1:1429670/139‑1 (MQ=255)                                                                                               ggCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTc                               >  1:1101415/1‑139 (MQ=255)                                                                                                  cAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCaaa                            >  1:502896/1‑139 (MQ=255)                                                                                                       tGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTc                                                                                                             >  1:304849/1‑53 (MQ=255)                                                                                                       tGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTc                                                                                                             <  2:304849/53‑1 (MQ=255)                                                                                                       tGCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAgg                       <  1:1335720/139‑1 (MQ=255)                                                                                                        gCCCTACCGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGgc                      <  1:1509740/139‑1 (MQ=255)                                                                                                              ccGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAaca                                                                                                         >  1:1238284/1‑50 (MQ=255)                                                                                                              ccGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAaca                                                                                                         <  2:1238284/50‑1 (MQ=255)                                                                                                               cGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGcc               <  2:947065/139‑1 (MQ=255)                                                                                                               cGCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGcc               <  2:895147/139‑1 (MQ=255)                                                                                                                gCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCa                                                               <  2:272792/90‑1 (MQ=255)                                                                                                                gCACGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCa                                                               >  1:272792/1‑90 (MQ=255)                                                                                                                   cGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAgg                       >  1:684008/1‑127 (MQ=255)                                                                                                                   cGTAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAgg                       <  2:684008/127‑1 (MQ=255)                                                                                                                     tAATGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAAt         <  1:566731/139‑1 (MQ=255)                                                                                                                        tGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                              >  2:1334856/1‑35 (MQ=255)                                                                                                                        tGTCAGCCCCAGCTTATCGCTGCGATTCAGATGCt                                                                                                              <  1:1334856/35‑1 (MQ=255)                                                                                                                            aGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCtt                                                                               <  1:98548/62‑1 (MQ=255)                                                                                                                            aGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCtt                                                                               >  2:98548/1‑62 (MQ=255)                                                                                                                            aGCCCCAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTa  >  2:590284/1‑139 (MQ=255)                                                                                                                              ccccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCa                                                               <  1:248712/76‑1 (MQ=255)                                                                                                                              ccccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCa                                                               >  2:248712/1‑76 (MQ=255)                                                                                                                              ccccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGca                                            <  2:977364/95‑1 (MQ=255)                                                                                                                              ccccAGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGca                                            >  1:977364/1‑95 (MQ=255)                                                                                                                                  aGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGt                                                                           >  2:764567/1‑60 (MQ=255)                                                                                                                                  aGCTTATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGt                                                                           <  1:764567/60‑1 (MQ=255)                                                                                                                                      |                                                                                                                                 TTCAGTTCATCGGTGCCTTTCGCCAGCGTCTCACTGCCATTTCTGACACTGGACACCTGGCTTGAGACATCGTTAATTAGCCAACGGCACATCAGGCCAAGTTGCCCTACCGCACGTAATGTCAGCCCCAGCTCATCGCTGCGATTCAGATGCTCAACACTATTACGTTCTCCGGTCGCCACCTTCAGTGCCTGATGGGCAACATTTTCTATCGGGCGCACAATCTGCCATTCAAAACAGGCGCTTGCCAACAATACCACTA  >  NC_000913/3218232‑3218493

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