Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,352,281	A→G	Y249H (TAC→CAC)	arcB ←	sensory histidine kinase ArcB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,352,281	0	A	G	100.0%	124.2 / NA	37	Y249H (TAC→CAC)	arcB	sensory histidine kinase ArcB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/21);  total (16/21)

GACTGATGGTGGAGATAAACGTCGTTTTGTCGCGGCTGGCCCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTGACACATTATGACGGAACACTTTTTCATCGGTTTCAATGACTTTTGCGGCGGCTTCCGGTGAG  >  NC_000913/3352152‑3352414                                                                                                                                  |                                                                                                                                      gACTGATGGTGGAGATAAACGTCGTTTTGTCGCGGCTGGCCCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTa                                                                                                                             <  2:559553/140‑1 (MQ=255)   cTGATGGTGGAGATAAACGTCGTTTTGTCGCGGCTGGCCCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACg                                                                                                                           <  1:435893/140‑1 (MQ=255)      aTGGTGGAGATAAACGTCGTTTTGTCGCGGCTGGCCCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGAt                                                                                                                        >  1:652586/1‑140 (MQ=255)      aTGGTGGAGATAAACGTCGTTTTGTCGCGGCTGGCCCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGAt                                                                                                                        >  2:62249/1‑140 (MQ=255)                tAAACGTCGTTTTGTCGCGGCTGGCCCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAg                                                                                                              >  2:161178/1‑140 (MQ=255)                tAAACGTCGTTTTGTCGCGGCTGGCCCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAg                                                                                                              >  1:799874/1‑140 (MQ=255)                   aCGTCGTTTTGTCGCGGCTGGCCCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCt                                                                                                           <  2:327653/140‑1 (MQ=255)                    cgtcgtTTTGTCGCGGCTGGCCCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCtt                                                                                                          >  2:64351/1‑140 (MQ=255)                         ttttGTCGCGGCTGGCCCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTgcgc                                                                                                     <  1:161178/140‑1 (MQ=255)                               cgcgGCTGGCCCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTcc                                                                                               >  1:415942/1‑140 (MQ=255)                                   gCTGGCCCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGt                                                                                           <  2:674977/140‑1 (MQ=255)                                   gCTGGCCCGTTCAAGCGCATCCTGACACCCCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGt                                                                                           <  1:214745/140‑1 (MQ=255)                                             tCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAACCCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCa                                                                                 <  2:73965/140‑1 (MQ=255)                                             tCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCa                                                                                 <  1:462058/140‑1 (MQ=255)                                             tCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCa                                                                                 <  1:489098/140‑1 (MQ=255)                                             tCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCa                                                                                 <  1:62249/140‑1 (MQ=255)                                             tCAAGCGCATCCTGATACCGCTGACGCTCAGTAATGTTGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCa                                                                                 <  1:259889/140‑1 (MQ=255)                                                       ccTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAg                                                                       <  1:110308/140‑1 (MQ=255)                                                         tGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGt                                                                     <  2:652586/140‑1 (MQ=255)                                                              ccGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTg                                                                <  2:689926/140‑1 (MQ=255)                                                               cGCTTACGCTCGGTAGTGTCGCGGCTAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTGa                                                               <  1:286183/140‑1 (MQ=255)                                                                gCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTGac                                                              <  2:799874/140‑1 (MQ=255)                                                                 cTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTGaca                                                             >  1:165704/1‑140 (MQ=255)                                                                         cGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTCTACGGATTTCAAAGCAGGCTTTGCGGCCGTCCGGATAAGCCAACCACTGTTTATAGCTCAGTGACACATTATGa                                                     >  2:286183/1‑140 (MQ=255)                                                                            tAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTGACACATTATGACgg                                                  <  2:415942/140‑1 (MQ=255)                                                                                      aCCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTGACACATTATGACGGAACACttttt                                        >  1:23960/1‑140 (MQ=255)                                                                                            gCCCATCAAACCGTGACTTTTACCCACGCGGTCGTAGTGCGGCACTTTTCGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTGACACATTATGACGGAACACTTTTTCATCgg                                  <  2:301096/140‑1 (MQ=255)                                                                                                 tCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTGAGTCAAACTTTATGACGGAACCATTTTTCACCGGTTTCa                             >  2:399754/1‑140 (MQ=255)                                                                                                 tCAAACCGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAATGACACATTATGACGGAACACATTTTCATCGGTTTCa                             >  1:349133/1‑140 (MQ=255)                                                                                                       cGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTGACACATTATGACGGAACACTTTTTCAt                                     <  1:725198/126‑1 (MQ=255)                                                                                                       cGTGACGTTTACCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTGACACATTATGACGGAACACTTTTTCAt                                     >  2:725198/1‑126 (MQ=255)                                                                                                                 aCCCACGCGGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTGACACATTATGACGGAACACTTTTTCATCGGTTTCAATGACTTTTGcggcgg             >  1:150870/1‑140 (MQ=255)                                                                                                                       gcgGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATGGGTCAGTGACACATTATGACGGAACACTTTTTCATCGGTTTCAATGACTTTTGCGGCGGCTTCCg       >  2:29080/1‑140 (MQ=255)                                                                                                                       gcgGTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTGACACATTATGACGGAACACTTTTTCATCGGTTTCAATGACTTTTGCGGCGGCTTCCg       >  1:98884/1‑140 (MQ=255)                                                                                                                         ggTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTGACACATTATGACGGAACACTTTTTCATCGGTTTCAATGACTTTTGCGGCGGCTTCCGGt     <  2:23960/140‑1 (MQ=255)                                                                                                                          gTCGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTGACACATTATGACGGAACACTTTTTCATAGGTTTCTATGACTTTTACGGCGGCTTCCGGTg    >  1:340414/1‑140 (MQ=255)                                                                                                                            cGTAGTGCGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTGACACATTATGACGGAACACTTTTTCATCGGTTTCAATGACTTTTGCGGCGGCTTCCGGTGAg  <  2:150870/140‑1 (MQ=255)                                                                                                                            cGTAGTGCGGCACTTTACGGATCTCAAAGCAAGCTTTGCGCCCGTCCGCGTAATCCAGCCACTGTTCATAGGTCAGTGACACATTATGACGGAACACTTTTTCATAGGTTTCAATGACTTTTGCGGCGGCTTCCGGTGAg  <  2:340414/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                      GACTGATGGTGGAGATAAACGTCGTTTTGTCGCGGCTGGCCCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACGGATTTCAAAGCAGGCTTTGCGCCCGTCCGGGTAATCCAGCCACTGTTCATAGGTCAGTGACACATTATGACGGAACACTTTTTCATCGGTTTCAATGACTTTTGCGGCGGCTTCCGGTGAG  >  NC_000913/3352152‑3352414

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