Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,991,678	T→C	S176P (TCC→CCC)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,991,678	0	T	C	100.0%	191.0 / NA	54	S176P (TCC→CCC)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/24);  total (30/24)

GCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGC  >  NC_000913/3991532‑3991822                                                                                                                                                   |                                                                                                                                                 gCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCAc                                                                                                                                              <  2:494919/151‑1 (MQ=255)  cGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCAcc                                                                                                                                             <  2:357651/151‑1 (MQ=255)  cGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCACCCTGGGGGGCGAAGATTGTGGCCCCAcc                                                                                                                                             >  1:305164/1‑151 (MQ=255)    aaGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCcccaccca                                                                                                                                           <  2:61982/151‑1 (MQ=255)     agagCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAg                                                                                                                                          <  1:612849/151‑1 (MQ=255)        gcgcCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCat                                                                                                                                       <  1:401541/151‑1 (MQ=255)         cgcCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCATCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCata                                                                                                                                      <  1:19212/151‑1 (MQ=255)          gcCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCatat                                                                                                                                     <  2:736090/151‑1 (MQ=255)           ccAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCatat                                                                                                                                     >  1:409290/1‑150 (MQ=255)            cAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATAc                                                                                                                                   <  2:599007/151‑1 (MQ=255)                    aCAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCtt                                                                                                                           >  2:726189/1‑151 (MQ=255)                     cAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTg                                                                                                                          <  1:775153/151‑1 (MQ=255)                          tAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGa                                                                                                                      >  2:135978/1‑150 (MQ=255)                           aaaTGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAAt                                                                                                                    <  2:578834/151‑1 (MQ=255)                            aaTGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAAtt                                                                                                                   <  1:540638/151‑1 (MQ=255)                             aTGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAAttt                                                                                                                  >  2:651381/1‑151 (MQ=255)                              tGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAAtttt                                                                                                                 >  1:733556/1‑151 (MQ=255)                               gTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTa                                                                                                                <  2:682849/151‑1 (MQ=255)                                     tgctgGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTAc                                                                                                          >  2:469513/1‑151 (MQ=255)                                         ggAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTAccgcc                                                                                                      >  1:62799/1‑151 (MQ=255)                                                        cTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGt                                                                                       >  1:684330/1‑151 (MQ=255)                                                              gggTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGt                                                                                 >  2:554574/1‑151 (MQ=255)                                                               ggTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTa                                                                                <  1:702062/151‑1 (MQ=255)                                                                      aaGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCtgtg                                                                         <  2:506423/151‑1 (MQ=255)                                                                          cAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAat                                                                     <  2:209071/151‑1 (MQ=255)                                                                          cAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAat                                                                     >  1:555259/1‑151 (MQ=255)                                                                              ttcttcCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATgg                                                                 >  2:179079/1‑151 (MQ=255)                                                                               tcttcCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATgg                                                                 >  2:570160/1‑150 (MQ=255)                                                                                 ttcCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGc                                                              >  2:136518/1‑151 (MQ=255)                                                                                  tcCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGcc                                                             >  2:573071/1‑151 (MQ=255)                                                                                   ccTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCg                                                            <  1:696926/151‑1 (MQ=255)                                                                                   ccTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCg                                                            >  2:536751/1‑151 (MQ=255)                                                                                      gattgatGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGc                                                         >  1:583171/1‑151 (MQ=255)                                                                                          gatgaAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGcgacg                                                     <  2:74957/151‑1 (MQ=255)                                                                                            tgaAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGcgacga                                                    >  1:501467/1‑150 (MQ=255)                                                                                                 aCCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAagaaga                                              <  1:762531/151‑1 (MQ=255)                                                                                                    gCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGc                                            >  1:612652/1‑150 (MQ=255)                                                                                                       tCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCAt                                          <  1:170533/149‑1 (MQ=255)                                                                                                         cGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTacg                                      >  2:440352/1‑151 (MQ=255)                                                                                                            cATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTacgacg                                   >  2:521371/1‑151 (MQ=255)                                                                                                            cATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTacgacg                                   <  2:684404/151‑1 (MQ=255)                                                                                                             aTAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTAcgacga                                  <  2:269302/151‑1 (MQ=255)                                                                                                               aaTGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACt                                <  2:270117/151‑1 (MQ=255)                                                                                                                aTGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTAcgacga                                  >  2:153027/1‑148 (MQ=255)                                                                                                                  gAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTAtg                             >  2:332394/1‑151 (MQ=255)                                                                                                                  gAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTAtg                             >  1:362535/1‑151 (MQ=255)                                                                                                                  gAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTAtg                             >  1:538084/1‑151 (MQ=255)                                                                                                                    aaGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGGATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTAtgtg                           <  2:713071/151‑1 (MQ=255)                                                                                                                     aGCGGCAGCCTGGGGGGCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTAtgtg                           >  2:649520/1‑150 (MQ=255)                                                                                                                                ggggggCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCttt                 <  1:504067/149‑1 (MQ=255)                                                                                                                                ggggggCGAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTAc               <  2:612709/151‑1 (MQ=255)                                                                                                                                       gAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAgg        >  1:427770/1‑151 (MQ=255)                                                                                                                                       gAAGATTGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAgg        >  1:312903/1‑151 (MQ=255)                                                                                                                                            ttGTGGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTg   >  1:178103/1‑151 (MQ=255)                                                                                                                                             tgtgGCCCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGc  >  2:709858/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGC  >  NC_000913/3991532‑3991822

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