Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,993,027	C→A	G625G (GGC→GGA)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,993,027	0	C	A	100.0%	157.2 / NA	45	G625G (GGC→GGA)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/29);  total (16/29)

GCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCA  >  NC_000913/3992897‑3993157                                                                                                                                   |                                                                                                                                   gCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCAt                                                                                                                           >  2:280713/1‑140 (MQ=255) gCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCAt                                                                                                                           <  2:436936/140‑1 (MQ=255)      ggCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGa                                                                                                                      <  2:325485/140‑1 (MQ=255)       gCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGAc                                                                                                                     >  2:611101/1‑140 (MQ=255)           gaATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAATATGCATCAGGAcgccg                                                                                                                 <  2:570051/140‑1 (MQ=255)           gaATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGAcgccg                                                                                                                 <  2:766357/140‑1 (MQ=255)           gaATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGAcgccg                                                                                                                 <  1:281553/140‑1 (MQ=255)            aaTCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGAcgccgc                                                                                                                >  1:722953/1‑140 (MQ=255)             aTCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCa                                                                                                               <  1:89225/140‑1 (MQ=255)                         tGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGc                                                                                                   >  2:143952/1‑140 (MQ=255)                          ggtaggtaGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCg                                                                                                  <  2:560107/140‑1 (MQ=255)                            taggtagCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTg                                                                                                >  1:308236/1‑140 (MQ=255)                             aggtagCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTgg                                                                                               <  1:442431/140‑1 (MQ=255)                             aggtagCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTgg                                                                                               <  1:280713/140‑1 (MQ=255)                             aggtagCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTgg                                                                                               <  2:556635/140‑1 (MQ=255)                                tagCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAg                                                                                            <  1:76319/140‑1 (MQ=255)                                    gTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGtctt                                                                                        <  2:246065/140‑1 (MQ=255)                                    gTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGtctt                                                                                        <  1:752634/140‑1 (MQ=255)                                    gTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGtctt                                                                                        <  1:504086/140‑1 (MQ=255)                                    gTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGtctt                                                                                        <  2:795186/140‑1 (MQ=255)                                     ttGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGtcttc                                                                                       >  2:601842/1‑140 (MQ=255)                                         ccTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGtt                                                                                   >  2:110695/1‑140 (MQ=255)                                             ctgTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAg                                                                               >  1:285261/1‑140 (MQ=255)                                               gTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGcc                                                                             <  1:131456/140‑1 (MQ=255)                                                         tCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTgcg                                                                   <  2:722953/140‑1 (MQ=255)                                                          cGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAAATATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTgcgc                                                                  >  2:776097/1‑140 (MQ=255)                                                            tGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCgg                                                                >  2:259245/1‑140 (MQ=255)                                                            tGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCgg                                                                <  2:364555/140‑1 (MQ=255)                                                                  gAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAAt                                                          <  1:716618/140‑1 (MQ=255)                                                                       gCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTa                                                     <  1:96011/140‑1 (MQ=255)                                                                       gCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTa                                                     >  1:453484/1‑140 (MQ=255)                                                                                gCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGAACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTgc                                            >  2:258858/1‑140 (MQ=255)                                                                                       aaCGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACt                                     <  2:565951/140‑1 (MQ=255)                                                                                              aGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCt                              >  1:608751/1‑140 (MQ=255)                                                                                               gCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTg                             <  2:720665/140‑1 (MQ=255)                                                                                                cAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGa                            <  1:110695/140‑1 (MQ=255)                                                                                                cAATCGATGATCGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGa                            <  2:659838/140‑1 (MQ=255)                                                                                                           cGAAGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCACCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAAt                 <  2:289513/140‑1 (MQ=255)                                                                                                             aaGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTg               >  2:194205/1‑140 (MQ=255)                                                                                                             aaGCCCTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTg               >  1:409371/1‑140 (MQ=255)                                                                                                                 ccTGAAAACTATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTc           >  2:737072/1‑140 (MQ=255)                                                                                                                          tATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCa  <  1:377046/140‑1 (MQ=255)                                                                                                                          tATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCa  <  1:295665/140‑1 (MQ=255)                                                                                                                          tATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCa  <  1:143952/140‑1 (MQ=255)                                                                                                                          tATTCTCGGAAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCa  <  2:800774/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                   GCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCA  >  NC_000913/3992897‑3993157

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑