Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,993,531	T→C	R793R (CGT→CGC)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,993,531	0	T	C	100.0%	93.8 / NA	32	R793R (CGT→CGC)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/17);  total (15/17)

CATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGA  >  NC_000913/3993286‑3993752                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                              cATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTAc                                                                                                                                                                                                                          <  2:321637/251‑1 (MQ=255)    aaCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTa                                                                                                                                                                                                                       <  1:368900/251‑1 (MQ=255)     aCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTAc                                                                                                                                                                                                                      <  2:241934/251‑1 (MQ=255)     aCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTAc                                                                                                                                                                                                                      <  1:235141/251‑1 (MQ=255)                 gCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTtc                                                                                                                                                                                                          <  1:3231/251‑1 (MQ=255)                       cAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGAGAAGGGATCATCCAGGTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGACAGCAACCGGGGTGAGGGATATCACCACTGCGAAAGCAGCAAAGAGGAGCCGGGACGTG‑ACGTCAGTCGCCTCTACTCGTCATCGCATGACCGCTTTACCTTCGGCTCAAGGTTCATCaa                                                                                                                                                                                                     >  1:39561/1‑250 (MQ=255)                       cAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGGGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGGTGAGGGATATCACCACTGCGAAGGGAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCCCAAGCTTCATCAAc                                                                                                                                                                                                    >  1:112888/1‑251 (MQ=255)                         gTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACtt                                                                                                                                                                                                  <  2:373343/251‑1 (MQ=255)                                  gAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGcc                                                                                                                                                                                         <  1:38419/251‑1 (MQ=255)                                     aCCAATCACGTCAAATTACCGGCGGTGGGGGACGGCTTTGCCAGCGAAGGGGTCATCCAGTTCTTTTTCGAAGAAACGCAAAACGAGAATGGGTTTTATTTCTACATTCTCGACGAAAACAACCGGGGTGAGGGATATCACCACTGCGAAGGCAGCCAAGAGGAGCCGGTACGTG‑ACGTCCGTCGCTTCTACTCCTCATCGCCTGACCGCTTTACCTACGGCCCAAGCTTCATCCACTTCAACCTGccccc                                                                                                                                                                                      >  2:314537/1‑248 (MQ=255)                                      ccAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGGACATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGGTGAGGGATATCCCCACTGCGAAGGCAGCCAAGAGGAGCTGGTACGTG‑GCGTCAGTCGCTTCTACTCGTCATCGCCTGACCGCCTTACCTACGGCTCAAGCTTCCTCAACTTCCACCTGCCGc                                                                                                                                                                                       >  1:151096/1‑249 (MQ=255)                                       cAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCGACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGt                                                                                                                                                                                    <  2:125158/251‑1 (MQ=255)                                                              tggtggACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATAGCTACATTCTCGAAGACCGCAACCGGGTTGAGGAATATAACCACTGCGAAGGCCGCAAAGAGGAACTGGTACTTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGTCTCAAGATTCATCAACATCAACCTGCCGCAGTTCTATCAGATTGTGCAGGTTGAt                                                                                                                                                             >  2:313049/1‑251 (MQ=255)                                                                  ggACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACAGTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTc                                                                                                                                                         <  2:291793/251‑1 (MQ=255)                                                                               aGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAt                                                                                                                                            <  1:314537/251‑1 (MQ=255)                                                                                       gATCATCCAGTTCTTTTTCGCAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAttccgttcc                                                                                                                                    >  2:320247/1‑251 (MQ=255)                                                                                       gATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAttccgttcc                                                                                                                                    >  1:40209/1‑251 (MQ=255)                                                                                                     ttttCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTa                                                                                                                        <  2:24363/249‑1 (MQ=255)                                                                                                          gaagaaACGCAAGACGCGAATGGCTTTCATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTACGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTa                                                                                                                  >  2:183012/1‑250 (MQ=255)                                                                                                                        cGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTCTATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGcc                                                                                                    >  2:75810/1‑250 (MQ=255)                                                                                                                               ggCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGa                                                                                            <  2:323141/251‑1 (MQ=255)                                                                                                                               ggCTTTAATATACACATTCTTGCCGAAAGCACCCGCGTTGAGGTATATCAACACAGCGCAGGCAGCAAAGACGAGCTGGCACCTG‑AAGTCAGCCGCTTATACTAGTCATCGCATGAACGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAgg                                                                                             >  2:37401/1‑250 (MQ=255)                                                                                                                                    tAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACg                                                                                       >  2:200542/1‑251 (MQ=255)                                                                                                                                      atatCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAACATCTATCGGTACCATGCCGCCTGCCAATCAGGATCACGAt                                                                                     >  2:31570/1‑251 (MQ=255)                                                                                                                                              aTTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCt                                                                             >  1:197015/1‑251 (MQ=255)                                                                                                                                                                 gggTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCGCGATACGCCGCTAGTACAGCAATATGTTTCg                                                          <  2:151096/251‑1 (MQ=255)                                                                                                                                                                            tatCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAc                                                                                   >  1:335557/1‑215 (MQ=255)                                                                                                                                                                            tatCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAc                                                                                   <  2:335557/215‑1 (MQ=255)                                                                                                                                                                                              agcaAAGGGGAGCTGGTAC‑TGTACTTCAGTCGTTTCTATGCGTCAACGCATGGCCGCTTTACCTACGGATCAAGCGTCATCAACTTCAACCTGCCGCAATTCTATCAGATTGGGAAGTGTGATGGTCGTGAACGGGTTATTCCGTTCCGCACAAAATCTATCGGTAACATGCAGCCTGCCAATCAGGATCACGATACGCCGCGATTACAGCAATATTGTTCGTGATGAACGTTCCTGAAAGCTAGGCTTAt                             <  2:248085/251‑1 (MQ=255)                                                                                                                                                                                               gcaAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATc                            <  2:112888/251‑1 (MQ=255)                                                                                                                                                                                                 aaaGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCg                          >  2:291254/1‑251 (MQ=255)                                                                                                                                                                                                                          aGTCGCTTCTACTCGTCAGCTCAGGACCGCTTTACCTTCGGCGCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACCGGTGATTCCGTGCCGCACAAAATCTATGGGTAACATGCTGCCTGCCAATCAGGATCACGATACGCCGCTATTCCAGCAATATTTTTCGTGATTGACGTGCCGGAATGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAactga  <  1:191016/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                              CATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTG‑ACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGA  >  NC_000913/3993286‑3993752

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