Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,993,461	A→G	H770R (CAC→CGC)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,993,461	0	A	G	100.0%	184.8 / NA	51	H770R (CAC→CGC)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/26);  total (25/26)

ATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGA  >  NC_000913/3993326‑3993596                                                                                                                                        |                                                                                                                                        aTCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATcgccac                                                                                                                                     >  1:1105960/1‑140 (MQ=255) aTCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAAAGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATcgccac                                                                                                                                     >  1:260811/1‑140 (MQ=255)    aCGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGc                                                                                                                                  <  2:603801/140‑1 (MQ=255)      gTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGa                                                                                                                                <  2:370986/140‑1 (MQ=255)      gTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGa                                                                                                                                <  1:85596/140‑1 (MQ=255)           aTTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGgca                                                                                                                           <  2:828851/140‑1 (MQ=255)                 ggcggTGGGGGGCGGCTTTGCCAGCGAAGGGGTCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTTATATCTACATTCTCGACGAAAGCAACCGGGGTGAGGGATATCGCCACTGCGAAGGCAGCaaaa                                                                                                                     >  2:746453/1‑139 (MQ=255)                     gtggtgGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAgagga                                                                                                                 <  1:1281262/140‑1 (MQ=255)                     gtggtgGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAAACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGGATATCGCCACTGCGAAGGCAGCAAAgaggg                                                                                                                 >  2:791413/1‑139 (MQ=255)                      tggtggACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAgaggag                                                                                                                <  2:145052/140‑1 (MQ=255)                        gtggACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCt                                                                                                              <  1:1317876/140‑1 (MQ=255)                         tggACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTg                                                                                                             >  1:1090342/1‑140 (MQ=255)                               gCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGt                                                                                                       >  2:1108239/1‑140 (MQ=255)                                cTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTg                                                                                                      >  1:933109/1‑140 (MQ=255)                                cTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTg                                                                                                      >  2:432437/1‑140 (MQ=255)                                  ttGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGAc                                                                                                    <  2:1362863/140‑1 (MQ=255)                                  ttGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGAc                                                                                                    >  1:1327532/1‑140 (MQ=255)                                   tGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACg                                                                                                   <  2:1105960/140‑1 (MQ=255)                                   tGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACg                                                                                                   <  2:1053747/140‑1 (MQ=255)                                        gCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGt                                                                                              <  1:1380656/140‑1 (MQ=255)                                                     ccAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGt                                                                                 <  2:257635/140‑1 (MQ=255)                                                         ttCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATc                                                                             >  1:39019/1‑140 (MQ=255)                                                         ttCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATc                                                                             >  1:447930/1‑140 (MQ=255)                                                                 cGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGAcc                                                                     <  1:432437/140‑1 (MQ=255)                                                                   aagaaACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGt                                                                   >  2:884935/1‑140 (MQ=255)                                                                   aagaaACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGt                                                                   >  2:605816/1‑140 (MQ=255)                                                                                      tGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCt                                                <  2:260811/140‑1 (MQ=255)                                                                                        gCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtc                                              >  2:1373423/1‑140 (MQ=255)                                                                                         cTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtca                                             <  2:950724/140‑1 (MQ=255)                                                                                             aaTATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCaa                                         >  2:562499/1‑140 (MQ=255)                                                                                                  cTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCa                                    >  2:900927/1‑140 (MQ=255)                                                                                                  cTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCa                                    >  1:345154/1‑140 (MQ=255)                                                                                                  cTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCa                                    >  2:994296/1‑140 (MQ=255)                                                                                                   tACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtcat                                            <  2:1301344/131‑1 (MQ=255)                                                                                                   tACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtcat                                            >  1:1301344/1‑131 (MQ=255)                                                                                                   tACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCaa                                   >  1:208130/1‑140 (MQ=255)                                                                                                   tACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCaa                                   >  1:284782/1‑140 (MQ=255)                                                                                                    aCATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAAc                                  >  1:582761/1‑140 (MQ=255)                                                                                                    aCATTCTCGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAAc                                  <  1:605816/140‑1 (MQ=255)                                                                                                        tctcGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGc                              >  2:305576/1‑140 (MQ=255)                                                                                                        tctcGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGc                              >  1:165185/1‑140 (MQ=255)                                                                                                        tctcGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGc                              >  2:690577/1‑140 (MQ=255)                                                                                                         ctcGACGAAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGcc                             <  1:884935/140‑1 (MQ=255)                                                                                                           cgacgaAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGc                           <  2:1090342/140‑1 (MQ=255)                                                                                                           cgacgaAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGc                           <  2:582761/140‑1 (MQ=255)                                                                                                           cgacgaAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGc                           <  1:746453/140‑1 (MQ=255)                                                                                                               gaAAGCAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGtt                       >  1:1393227/1‑140 (MQ=255)                                                                                                                    cAACCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATc                  >  1:733653/1‑140 (MQ=255)                                                                                                                     aaCCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCa                 <  1:228945/140‑1 (MQ=255)                                                                                                                     aaCCGGGTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCa                 <  1:305576/140‑1 (MQ=255)                                                                                                                           gTTGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATtgt           >  1:1378487/1‑140 (MQ=255)                                                                                                                            ttGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATtgtg          >  1:87901/1‑140 (MQ=255)                                                                                                                            ttGAGGTATATCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATtgtg          <  2:208130/140‑1 (MQ=255)                                                                                                                                    tatCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGa  <  1:791413/140‑1 (MQ=255)                                                                                                                                    tatCGCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGa  <  2:733653/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        ATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGA  >  NC_000913/3993326‑3993596

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