Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,992,624	T→G	100%	L491R (CTG→CGG)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,992,624	0	T	G	100.0%	246.0 / NA	64	L491R (CTG→CGG)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (32/32);  total (32/32)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGC  >  NC_000913/3992486‑3992762                                                                                                                                           |                                                                                                                                           cGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTg                                                                                                                                   <  4:118518/148‑1 (MQ=255)   gAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACggg                                                                                                                                          <  3:315716/139‑1 (MQ=255)   gAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACggg                                                                                                                                          >  4:315716/1‑139 (MQ=255)              gACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAAc                                                                                                                      >  2:489009/1‑148 (MQ=255)                ccTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACgg                                                                                                                    <  1:489009/148‑1 (MQ=255)                                gccgggccggGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCa                                                                                                                               >  3:1106644/1‑121 (MQ=255)                                gccgggccggGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCa                                                                                                                               <  4:1106644/121‑1 (MQ=255)                                      ccgggcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGt                                                                                              >  1:883952/1‑148 (MQ=255)                                            tAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGTATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGtat                                                                                        <  2:883952/148‑1 (MQ=255)                                              aCCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTg                                                                                                                                   >  2:752448/1‑103 (MQ=255)                                              aCCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTa                                                                                                                        >  4:701974/1‑114 (MQ=255)                                              aCCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTa                                                                                                                        <  3:701974/114‑1 (MQ=255)                                              aCCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCAATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTg                                                                                                                                   <  1:752448/103‑1 (MQ=255)                                               ctgTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCctg                                                                                                                <  4:424480/119‑1 (MQ=255)                                               ccGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCctg                                                                                                                >  3:424480/1‑121 (MQ=255)                                                  ttCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGtt                                                                                             >  4:317901/1‑137 (MQ=255)                                                  ttCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGtt                                                                                             <  3:317901/137‑1 (MQ=255)                                                   tCAGGTTGGTATCTTTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggt                                                                                                                                         <  1:336108/92‑1 (MQ=255)                                                   tCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggt                                                                                                                                         >  2:336108/1‑92 (MQ=255)                                                   tCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggt                                                                                                                                         <  3:1020241/92‑1 (MQ=255)                                                   tCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggt                                                                                                                                         >  4:1020241/1‑92 (MQ=255)                                                        ttGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTAtaataa                                                                                   >  4:861529/1‑141 (MQ=255)                                                        ttGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaa                                                                                   <  3:861529/141‑1 (MQ=255)                                                            tATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATg                                                                                                                             >  4:420255/1‑95 (MQ=255)                                                            tATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATg                                                                                                                             <  3:420255/95‑1 (MQ=255)                                                            tATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGtt                                                                                             >  4:1093534/1‑127 (MQ=255)                                                            tATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGtt                                                                                             <  3:1093534/127‑1 (MQ=255)                                                              tCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACg                                                                           >  2:1004380/1‑143 (MQ=255)                                                              tCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACg                                                                           <  1:1004380/143‑1 (MQ=255)                                                               cTGTATAACCGCGCGCCAAATGTTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGtt                                                                                             >  1:686476/1‑124 (MQ=255)                                                               cTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTaa                                                                                                                       >  3:1067520/1‑98 (MQ=255)                                                               cTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTaa                                                                                                                       <  4:1067520/98‑1 (MQ=255)                                                               cTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGtt                                                                                             <  2:686476/124‑1 (MQ=255)                                                                  tataACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCctgc                                                                                                               <  2:64890/103‑1 (MQ=255)                                                                  tataACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCctgc                                                                                                               >  1:64890/1‑103 (MQ=255)                                                                        cgcgcgCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAgg                                                                                >  2:563481/1‑128 (MQ=255)                                                                        cgcgcgCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAgg                                                                                <  1:563481/128‑1 (MQ=255)                                                                                aaaTATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGc                                                                                                                                <  1:335247/72‑1 (MQ=255)                                                                                aaaTATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGc                                                                                                                                >  2:473436/1‑72 (MQ=255)                                                                                aaaTATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGc                                                                                                                                >  2:335247/1‑72 (MQ=255)                                                                                aaaTATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGc                                                                                                                                <  1:473436/72‑1 (MQ=255)                                                                                aaaTATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGt                                                                                          >  3:298942/1‑110 (MQ=255)                                                                                aaaTATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGt                                                                                          <  4:298942/110‑1 (MQ=255)                                                                                 aaTATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCa                                                                                                   >  4:1116161/1‑100 (MQ=255)                                                                                 aaTATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCa                                                                                                   <  3:1116161/100‑1 (MQ=255)                                                                                          tCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTcg                                                                                                       <  3:761422/87‑1 (MQ=255)                                                                                          tCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTcg                                                                                                       >  4:761422/1‑87 (MQ=255)                                                                                              tcatcaGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATg                                                                                                                             <  3:632700/61‑1 (MQ=255)                                                                                              tcatcaGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATg                                                                                                                             >  4:632700/1‑61 (MQ=255)                                                                                              tcatcaGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGtt                                                                                             <  4:877495/93‑1 (MQ=255)                                                                                              tcatcaGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGtt                                                                                             >  3:877495/1‑93 (MQ=255)                                                                                              tcatcaGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGc                                      >  4:869429/1‑148 (MQ=255)                                                                                                   agccatcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGt                                                                                                                           <  2:894903/58‑1 (MQ=255)                                                                                                   agccatcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGt                                                                                                                           >  1:894903/1‑58 (MQ=255)                                                                                                         cagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCg                                                                                               <  1:598004/80‑1 (MQ=255)                                                                                                         cagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCg                                                                                               >  2:598004/1‑80 (MQ=255)                                                                                                                   aaTATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTa                                                         >  4:383698/1‑108 (MQ=255)                                                                                                                   aaTATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTa                                                         <  3:383698/108‑1 (MQ=255)                                                                                                                    atatAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCg                <  3:869429/148‑1 (MQ=255)                                                                                                                      ataACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCt              >  3:243512/1‑148 (MQ=255)                                                                                                                        aaCCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTaa                                                                             >  1:783481/1‑83 (MQ=255)                                                                                                                        aaCCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTaa                                                                             <  2:783481/83‑1 (MQ=255)                                                                                                                                  tGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGgag                                              <  2:625505/104‑1 (MQ=255)                                                                                                                                  tGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGgag                                              >  1:625505/1‑104 (MQ=255)                                                                                                                                  tGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTa       >  1:160602/1‑143 (MQ=255)                                                                                                                                  tGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTa       <  2:160602/143‑1 (MQ=255)                                                                                                                                  tGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGc  <  4:243512/148‑1 (MQ=255)                                                                                                                                  tGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGc  >  3:441314/1‑148 (MQ=255)                                                                                                                                   gAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGAt                                            <  3:908857/105‑1 (MQ=255)                                                                                                                                   gAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGAt                                            >  4:908857/1‑105 (MQ=255)                                                                                                                                           |                                                                                                                                           CGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGC  >  NC_000913/3992486‑3992762

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