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%	108.4 / NA	34	intergenic (‑52/‑366)	aer/patA	aerotaxis receptor/putrescine aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/15);  total (19/15)

GCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTTCACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCAATTATTGCCCTGAAAACGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAA  >  NC_000913/3218886‑3219358                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       ccTTTTTGCCTATCCGGGTGACGCCCCATGTTGTGCGGGTGCCTTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTAGCGTCCGCCAGCGGGGTATTTGGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCCGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGcccc                                                                                                                                                                                                                                <  2:346688/250‑1 (MQ=255)    tttGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCg                                                                                                                                                                                                                               <  1:513765/249‑1 (MQ=255)     ttGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCt                                                                                                                                                                                                                             <  1:44562/250‑1 (MQ=255)            atCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTAt                                                                                                                                                                                                                     >  1:294546/1‑251 (MQ=255)               cGGGTGACGCACCATGTTGTGCCGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCAGGAGTGATATAGCGTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGTCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTGCGCGTCCTGTGCAGGTTATATGTTTGTTAACTTCTTGTCAGAGGTTATATCGGCCCCGCTGCGGTTATCtt                                                                                                                                                                                                                  <  2:43392/251‑1 (MQ=255)                           cATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAAt                                                                                                                                                                                                      <  1:316636/251‑1 (MQ=255)                            aTGTTGTGCGGCTGCCATTGTAACTCTTGCAAGGTATAGCCGCGCACCAGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGAACGGAAGTGGACATCAGAGAGGGATCGTACGCAAGCGGGGAGTTTGGCTGGGTGACATAAGGCTGAGAAGACATGTTCGCGTACTGTGCAGGATATATGGGTGTTGACTTCTTGTCAGAGTTTATATAGGCCCCGCTTCGGTTGTCTTTAACCGATTactt                                                                                                                                                                                                     >  2:63207/1‑248 (MQ=255)                                          cccTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCa                                                                                                                                                                                       <  1:340680/251‑1 (MQ=255)                                             ttGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAAtt                                                                                                                                                                                    >  1:314375/1‑251 (MQ=255)                                               gtaacatgtGCAAGGTATAGCCGTTCACCTGCACAAAAGTGTCATTAGCATGAGTTATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTTGTATCTTCCGCCAGCGGGGTATTTTGCTGGGTGACATTCGTATGAGAAGACATGTCCGCATCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTTTATCCGCCCCGCTGCGGTTATCATTAACTGATTAATTTGATTTAGATCGCAATTTg                                                                                                                                                                                  <  1:213890/243‑1 (MQ=255)                                                       tGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTa                                                                                                                                                                           >  1:350826/1‑250 (MQ=255)                                                                     gCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGTTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAAtt                                                                                                                                                            <  1:316500/251‑1 (MQ=255)                                                                                 aaaaGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGTTTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTaa                                                                                                                                                 >  2:212807/1‑250 (MQ=255)                                                                                        tcagtAGCATGAGTGATATAGCTTTGCAGATCTGTAGTGGTCATCAGAGTGGTATTGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACtt                                                                                                                                          <  2:541488/246‑1 (MQ=255)                                                                                           ttAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTTCAc                                                                                                                                      >  1:170077/1‑251 (MQ=255)                                                                                               cATGAGTGATAAAGATTTGACGATCGGTAGTGGACAGCAGAGTGGGATCGTCCGCCAGAGGGGTATTTTGCTGGGTGAAATACGGATGAGCAGACAAGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTATTGTCAGAGTTTATATCGGCACAGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGTTCGCAATTTGCGATTTAAACAAAAATATAAGTCCTTAATTTAAAATACtttacctatt                                                                                                                                  >  2:389270/1‑244 (MQ=255)                                                                                                aTGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTTCACTCTGt                                                                                                                                 <  1:262943/251‑1 (MQ=255)                                                                                                     tGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCACTTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTTAACTCTGTTACt                                                                                                                             >  2:200043/1‑250 (MQ=255)                                                                                                      aagataacttttaaGATCGGTAGTGGACATCAGAGTGGTTTCGTCCTCCAGCGGGGTATTTTGCTGGGTGACATACAGATGAGCAGACATGTTCGCGCCCTGGGCGGGTTATATGGTTGTTAACTTCTTGTAATAGGTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGCTTTAGTGCACACTTTGAGATTTAAACACAAATCTAATTCCTTTATTTAACAAACTTTCACTCTGTTACtat                                                                                                                           >  2:358691/14‑251 (MQ=255)                                                                                                                   aGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTTCACTATGTTACTATAAGAAACCGTTa                                                                                                               >  1:302469/1‑250 (MQ=255)                                                                                                                     aTCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTTCACTCTGTTACTATACGACAACGTTAAtt                                                                                                            >  1:184751/1‑251 (MQ=255)                                                                                                                         gtggTGGACAGCAGAGTGGTATTTTCCGCCAGTGGGGTATTTTGCTTGTTTACATACGGATGAGAAGACATGTTCGCATCATGTGCAGGTGATATGGTTGTTAACTATTTGTCAGAGTTTATATCTGCCCCTCTGCGGTTAACTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATGCCTTTATTTAAAATACTTTCTCTCTGTTACTATACGAAAACGTTAATTATca                                                                                                        <  2:203880/248‑2 (MQ=255)                                                                                                                               gcaatacgtgttgTATCGTCCTCCCGCGTGGTATTTGGCTGGGATACAGACGGGTGAGAAGACACGTTTGTGCCCTGGGCAGGTTATATTTTTGTTAACTTCTTGTCAGAGTGTATAGCGGCTCATCTGCGGTTATCTTTAACCGATTAATTTGATGTTGATCGCAATTTGCGATTGAAACTCAAATCTAATTCCTTTATTTAAAATGCTTTCTCCCTGTTTTTATATGAAAACGTTAAGTATCATGcgca                                                                                                  <  2:245436/239‑4 (MQ=255)                                                                                                                                      aagaggTATCGTCCGCCAGCGGGGAATTTTGCTGGGAGACATACGGGTGATCAGACATGTTCGCGTCCTGTGCAGGGTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTACTTTTATTTAGATCGCAATTTGCGATTTAAACACAACTCTAATTCCTTGATTAAAAATACTTTCACTCTGTTACTATACGAAACCGTTAATTATCTTGCCCAACAATc                                                                                            >  2:125013/5‑250 (MQ=255)                                                                                                                                         tGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTTCACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGc                                                                                        >  1:109572/1‑251 (MQ=255)                                                                                                                                             aTCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTTCACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCAatt                                                                                    >  1:496123/1‑251 (MQ=255)                                                                                                                                                gTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTTCACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCAattatt                                                                                 >  1:3284/1‑251 (MQ=255)                                                                                                                                                          ggggTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTTCACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCAATTATTGCCCTGaaa                                                                        >  1:93736/1‑250 (MQ=255)                                                                                                                                                            ggTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTTCACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAACAATCAGGCACTTATTGCCCTGAAAACg                                                                     >  1:42708/1‑251 (MQ=255)                                                                                                                                                                                            acttGTTCGCGTCCTTTGCAGGTGATATGGTTGTTAACTTCTTTTCAGAGTTTATATCGGCCCTGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTTATTTAAAATACTTTCCCTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCAATTATTGCCCTGAAAACGTGCATTTGCGCAGCAATCATCAAATCCATAcc                                     <  2:140373/248‑1 (MQ=255)                                                                                                                                                                                                                 gTTATATGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTTCACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCAATTATTGCCCTGAAAACGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAAt                 >  1:354856/1‑250 (MQ=255)                                                                                                                                                                                                                    gtatGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGATATCTTTAACCGATTAATTTGATTTAGATCGCCATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTTCACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCAATTATTGCCCTGAAAACGTGCATTTGCGCAGCAATCATCAAATCCATACCAGACAAAAACCGGGCAAAATaac              >  2:404672/2‑250 (MQ=255)                                                                                                                                                                                                                     tatGGTTGTTAACTTCTTGTCAGAGTTTATATCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTTCACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCAATTATTGCCCTGAAAACGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATaaca             <  2:294546/250‑1 (MQ=255)                                                                                                                                                                                                                                tcttcctGTTAGAGTTTATATCGGCCACGTTTAGGTTATCTTTAACCGATTAAGTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTTATTTAAAATTCTTTCACTCTTTTACTATCCGAATACGTTAATTATCTCGCCCAAAAATCAGGCAATGATTGCCCTGAAAACGTGCATTTGCGCAGCAATCATCAAATCCATACCCGGCAAAAACCGTGCAAAATAAACACAAAGGTTaa  <  2:43982/249‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       GCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTTCACTCTGTTACTATACGAAAACGTTAATTATCTTGCCCAAAAATCAGGCAATTATTGCCCTGAAAACGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAA  >  NC_000913/3218886‑3219358

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