Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	325,394	G→A	intergenic (‑30/+183)	ykgH ← / ← betA	uncharacterized protein YkgH/choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	325,394	0	G	A	100.0%	115.8 / NA	35	intergenic (‑30/+183)	ykgH/betA	uncharacterized protein YkgH/choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/14);  total (21/14)

AAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTC  >  NC_000913/325149‑325634                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 aaaaTTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAg                                                                                                                                                                                                                                             <  1:499569/251‑1 (MQ=255)             gCTATCACAGCTAACAAGAATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCCCTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTTTCGCAAATAAGGTTGTATTAAt                                                                                                                                                                                                                                  >  2:381305/1‑250 (MQ=255)                aTCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGGAAATCACTTTTATATTTTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACTCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGt                                                                                                                                                                                                                               <  2:422964/250‑1 (MQ=255)                         aaCAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAg                                                                                                                                                                                                                      <  2:181795/250‑1 (MQ=255)                            aaGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATCATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCAGTATCAGGAAt                                                                                                                                                                                                                  >  1:20706/1‑251 (MQ=255)                              gTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGt                                                                                                                                                                                                                >  2:378867/1‑251 (MQ=255)                              gTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGt                                                                                                                                                                                                                <  2:18788/251‑1 (MQ=255)                                 taCCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCg                                                                                                                                                                                                             >  2:231638/1‑251 (MQ=255)                                                                gTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGcatcat                                                                                                                                                                               <  1:94181/250‑1 (MQ=255)                                                                   aaTCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAAt                                                                                                                                                                           <  1:112061/251‑1 (MQ=255)                                                                                                aCAGAAACAACAGCACCATAGCCATACATATAGCCATACTAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCAAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCACGCATCATGAATACTGGATATGAAGCATGAGAGATACCTc                                                                                                                                               >  2:509371/1‑250 (MQ=255)                                                                                                aCAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTc                                                                                                                                               >  2:175434/1‑250 (MQ=255)                                                                                                  aGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAgt                                                                                                                                            <  1:112970/251‑1 (MQ=255)                                                                                                     aacaacAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGGATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTTCCTCAGTGtt                                                                                                                                         <  2:361158/251‑1 (MQ=255)                                                                                                      acaacaGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGttt                                                                                                                                        >  1:300587/1‑251 (MQ=255)                                                                                                             cACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTAACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTtatata                                                                                                                                   >  1:42562/1‑249 (MQ=255)                                                                                                                             ataGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCctct                                                                                                                 >  2:260133/1‑251 (MQ=255)                                                                                                                                    aaataaAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAAc                                                                                                          <  2:78929/250‑1 (MQ=255)                                                                                                                                     aataaAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCACTCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATACAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGATCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGCTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACg                                                                                                         >  1:485585/1‑251 (MQ=255)                                                                                                                                                               tcttcAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCt                                                                                <  1:378867/250‑1 (MQ=255)                                                                                                                                                                  tcAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTg                                                                               <  2:270111/248‑1 (MQ=255)                                                                                                                                                                      ataaaataaaATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGt                                                                        >  2:491986/1‑251 (MQ=255)                                                                                                                                                                            taaaatCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACAt                                                                  >  1:47739/1‑251 (MQ=255)                                                                                                                                                                                 tctcAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACg                                                              >  1:173372/1‑250 (MQ=255)                                                                                                                                                                                        cAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAttt                                                      >  1:389520/1‑251 (MQ=255)                                                                                                                                                                                        cAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAtt                                                       <  1:123735/250‑1 (MQ=255)                                                                                                                                                                                             cGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCg                                                 >  1:252087/1‑251 (MQ=255)                                                                                                                                                                                             cGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCg                                                 <  1:443371/251‑1 (MQ=255)                                                                                                                                                                                                      gTTTGATTTGTTAACGCATCATATATCCAGAGAATAAACTCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGATCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCATCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATAACGACTCATTTTTTCGCTCTCACCg                                        >  2:361928/1‑251 (MQ=255)                                                                                                                                                                                                           aTTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATc                                   >  1:303859/1‑251 (MQ=255)                                                                                                                                                                                                              tGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATccc                                 <  1:231638/250‑1 (MQ=255)                                                                                                                                                                                                                         aataTATCCAGAGAATAAAATCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCACCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACaa                      >  1:260863/1‑250 (MQ=255)                                                                                                                                                                                                                          ctatatCCAGAGCATAAAATCGATTGCAATTAAGGCTGTATTCATATTCTGTATCCGGAATGTTCGGGATAAATATAAGAAAAAAACCCGCATCATGAATACTGGATATGAAGCATGCGAGATACCTAAGAGTTTATATAAGGATTCGGTCCCCCTCCCTGGAACGGTACATCCCAATCTGATTGATTCCTGCGTTAGTTCACATCACGACTCATTTATTCGCTCTCACCGGCATACCATTTGCCACaaa                     >  2:393690/2‑250 (MQ=255)                                                                                                                                                                                                                                            aTCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAACAAGCCCGCCTCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACCAAATATCCCGCCGTGctcctc  >  2:572957/1‑251 (MQ=255)                                                                                                                                                                                                                                            aTCTGTCGCAAATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGctcctc  >  1:360947/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 AAAATTATTATTGCTATCACAGCTAACAAGTATACCTGGCGTTGTCTGGTTTATTTTTAGGCTGTAAATCACTTTTATATTGTCTTTATTGATATACAGAAACAACAGCACCATAGCCATACATATAGCCATAATAAAAAGGATAACACGAATCTTTTTCTTCAGATAAAATAAAATCTCAATCAGATCGATATCCTGTTTGATTTGTTCACGCATAATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTC  >  NC_000913/325149‑325634

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