Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,991,903	C→T	L251F (CTC→TTC)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,991,903	0	C	T	100.0%	102.6 / NA	30	L251F (CTC→TTC)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/12);  total (18/12)

GGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCT  >  NC_000913/3991663‑3992137                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           ggCGAAAATTGTGGCTCCACCAAGCATATACTGCTGCTTGACGAATTTAATCGTACCGCCGTGCGTCGCGCCGGTAAGCGTATTCTGTGGAATATGGCGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGttctgcacgag                                                                                                                                                                                                                                  <  2:195874/251‑8 (MQ=255)     aaGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATc                                                                                                                                                                                                                              >  1:237723/1‑251 (MQ=255)                                  ctgctTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAac                                                                                                                                                                                                 >  2:77378/1‑251 (MQ=255)                                     ctTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACAct                                                                                                                                                                                              >  2:150384/1‑251 (MQ=255)                                     ctTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATATCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACAct                                                                                                                                                                                              >  2:1697/1‑251 (MQ=255)                                       tGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACAct                                                                                                                                                                                              <  2:154818/249‑1 (MQ=255)                                                     gTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAt                                                                                                                                                                              >  1:69649/1‑251 (MQ=255)                                                      tACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACCCTGCTGCTGGAAGCCTAtt                                                                                                                                                                             >  2:198214/1‑251 (MQ=255)                                                          gccgTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTg                                                                                                                                                                         >  1:136526/1‑251 (MQ=255)                                                                                       gtgGAATATGGAGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCAAGTCTGCTGGCg                                                                                                                                             >  1:69512/1‑250 (MQ=255)                                                                                           aaTATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGaa                                                                                                                                           >  2:121298/1‑248 (MQ=255)                                                                                              atGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGatat                                                                                                                                     <  1:101272/251‑1 (MQ=255)                                                                                                    ccGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAAc                                                                                                                                >  1:161348/1‑250 (MQ=255)                                                                                                          gacgaAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGttt                                                                                                                         >  1:203354/1‑251 (MQ=255)                                                                                                                                    gtttatgaCGCTTTACGCGCAGGTCGGGCTGACGCCAAATGAATGGCTGGATCGCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTTCTTTGGTGACAGCCTTTGGCAGTTCTACAAGAGTATTGATTCCCCATACAAAGCGTTACGTAGATCTCGGGTTCTGGATTCCTGTTCCTGGTGATACCCGAACCCAGGTCTGCTGGCGAAAGATATCACACAGCGAGTTCGCGATGGCGAGATTGTATCGTTTg                                                                                               <  2:248855/247‑1 (MQ=255)                                                                                                                                       gatgaCGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGGTTGGtc                                                                                            <  1:89894/251‑1 (MQ=255)                                                                                                                                            cGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAt                                                                                       >  1:171668/1‑251 (MQ=255)                                                                                                                                             gCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTAGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATc                                                                                      <  2:68781/251‑1 (MQ=255)                                                                                                                                                  aCGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATa                                                                                  <  1:201692/250‑1 (MQ=255)                                                                                                                                                     cgcAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGc                                                                              >  2:164345/1‑251 (MQ=255)                                                                                                                                                                  gACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGTTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAAGACCCGAACCCACGTCTGCTGGCGAAAGATGTCAAACAGTGTTTGCACGACTGCTAGATTTTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGa                                                                   <  1:228628/249‑1 (MQ=255)                                                                                                                                                                                 gCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATAcc                                                  >  1:144838/1‑251 (MQ=255)                                                                                                                                                                                      gtctcGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACg                                             <  2:135473/250‑1 (MQ=255)                                                                                                                                                                                        ctcGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAACCAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCAATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGCCGGc                                           >  2:190258/1‑251 (MQ=255)                                                                                                                                                                                                     gCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGAtttt                              <  1:54427/251‑1 (MQ=255)                                                                                                                                                                                                               tgctgAAGAGTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTgg                    >  1:41733/1‑251 (MQ=255)                                                                                                                                                                                                                       agTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGAGTGAAGATTTTACCCGTCTGGATTTAGTa            <  2:237723/251‑1 (MQ=255)                                                                                                                                                                                                                        gTACTTTGGTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTAc           >  2:167112/1‑251 (MQ=255)                                                                                                                                                                                                                               ggTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCgctg    <  1:133826/251‑1 (MQ=255)                                                                                                                                                                                                                               ggTGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCgct     <  1:34616/250‑1 (MQ=255)                                                                                                                                                                                                                                 tGCCAGCCTTTGGCAGTTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCgctgct  >  2:24785/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           GGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCT  >  NC_000913/3991663‑3992137

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