Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,219,128	G→A	intergenic (‑52/‑366)	aer ← / → patA	aerotaxis receptor/putrescine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,219,128	0	G	A	100.0%	77.6 / NA	24	intergenic (‑52/‑366)	aer/patA	aerotaxis receptor/putrescine aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/10);  total (14/10)

CGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT‑CACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCAATTATTGCCCTGAAAACGTGCATTTGCGCAGCAATCATCAAATCCAT  >  NC_000913/3218885‑3219320                                                                                                                                                                                                                                                    |                                                                                                                                                                                                  cgcTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCATTGTAACTCTTGCAAGGGATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGAGACATGCGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGccc                                                                                                                                                                                            >  2:324716/1‑251 (MQ=255)                         cACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGAtt                                                                                                                                                                    <  2:132766/251‑1 (MQ=255)                           ccATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATCTCGCTTTGCCGAACGGTAGTGGACATCAGAGTGGTCTCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACAGGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTa                                                                                                                                                                   >  1:401420/1‑250 (MQ=255)                                                gTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTg                                                                                                                                             <  1:74221/251‑1 (MQ=255)                                                             ggTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGCTTAATTTGATTTAGATCGCAATTTGCGATTTAAacac                                                                                                                                 >  2:270049/1‑250 (MQ=255)                                                              gTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGGGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAacaaa                                                                                                                                >  2:14476/1‑248 (MQ=255)                                                                 taGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATc                                                                                                                            >  1:184362/1‑251 (MQ=255)                                                                                       tgtCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTCTATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCATTGATTTAAAATACt                                                                                                      >  1:148372/1‑251 (MQ=255)                                                                                                  tgactGATATAGCTTTGCAGATCGGAAGTGGACATCAGCGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCACCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT‑CACTCTGtt                                                                                          >  2:274705/5‑251 (MQ=255)                                                                                                   gAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT‑CACTCTGTTa                                                                                         <  1:316990/251‑1 (MQ=255)                                                                                                      tGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT‑CACTCTGTTACt                                                                                       >  1:329880/1‑250 (MQ=255)                                                                                                       gATATAGCTTTTCAGATCGTTTGTGGTCAGCAGATTGGTTGCGTCCGCCAGCGGTGTATTTTGCTTGGTTACATTCGGATGAGAAGACATGTTCGCGTCCTGTGCATGTTATATGGTTGTTAACTGCTTGTCAGTGTTTATATCGTCCCGGCTGCGGTTATCTTTATCCGATTAATTTGATGTATATCTCAATGTGCGATTTAACCACAAATCTAATTCCTTGATGTAAA‑TGCTTTTCTCTCTGTTACta                                                                                      <  2:256264/250‑1 (MQ=255)                                                                                                                                     cAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT‑CACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAAt                                                       >  1:73355/1‑251 (MQ=255)                                                                                                                                      agagTGGTAGCGTCCGCCAGCGGGGTATTTTGCTGGTTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT‑CACTCTGTTACTATACGAAAACGTTAATTATCGTGCCCAAAAATc                                                      <  2:184362/251‑1 (MQ=255)                                                                                                                                            gTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGCCATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT‑CACTCTGTTACTATACGAAAACg                                                                            >  2:228500/1‑223 (MQ=255)                                                                                                                                            gTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT‑CACTCTGTTACTATACGAAAACg                                                                            <  1:228500/223‑1 (MQ=255)                                                                                                                                                          cGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGGTTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT‑CACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCAATTATTGCCCTGaaa                                  <  1:167017/251‑1 (MQ=255)                                                                                                                                                             ggTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCATTGATTTAAAATACTTT‑CACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCGATTATTGCCCTGAAAACg                               >  1:192611/1‑251 (MQ=255)                                                                                                                                                                   ttGCTGGGTGACATGCGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTATCCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT‑CACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCAATTATTGCCCTGAAAACGTGCAt                          <  1:43958/250‑1 (MQ=255)                                                                                                                                                                      cTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT‑CACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCAATTATTGCCCTGAAAACGTGCATTTgc                      >  2:130893/1‑251 (MQ=255)                                                                                                                                                                            gccATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT‑CACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCAATTATTGCCCTGAAAACGTGCATTTGCgcagc                 <  2:104242/248‑1 (MQ=255)                                                                                                                                                                                           agaCATGATCGCGTCCTGTGCAGATTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT‑CACTCTGTTACTATACGAAAACGTTAATTATCTTGCCAAAAAATCAGGCAATTATTGCCCTGAAAACGTGCATTTGCGCAGCAATCATCAAATCCAt  >  2:269210/1‑250 (MQ=255)                                                                                                                                                                                            gaCATGTTCGCGTCCTGTGCATGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATGTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATTCTTT‑CACTCTg                                                                                            <  2:28343/159‑1 (MQ=255)                                                                                                                                                                                            gaCATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGCGTTTATATCGGCCCCGCTGCGGTTATCTTTCACCGCTTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT‑CACTCTg                                                                                            >  1:28343/1‑159 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                  CGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT‑CACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCAATTATTGCCCTGAAAACGTGCATTTGCGCAGCAATCATCAAATCCAT  >  NC_000913/3218885‑3219320

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