Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,742,962	G→A	D194N (GAT→AAT)	yfiN →	putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,742,962	0	G	A	100.0%	163.1 / NA	51	D194N (GAT→AAT)	yfiN	putative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/26);  total (25/26)

CGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCAGCTTCGTTTA  >  NC_000913/2742826‑2743096                                                                                                                                         |                                                                                                                                       cgACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAAt                                                                                                                                      <  1:1216120/139‑1 (MQ=255)                ccATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGAtgtg                                                                                                                       >  2:606715/1‑139 (MQ=255)                 cATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGc                                                                                                                      <  1:879991/139‑1 (MQ=255)                          tGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAAcc                                                                                                             <  1:514576/139‑1 (MQ=255)                              tttCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTa                                                                                                                                >  2:480226/1‑116 (MQ=255)                              tttCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTa                                                                                                                                <  1:480226/116‑1 (MQ=255)                              tttCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCaa                                                                                                         >  2:174587/1‑139 (MQ=255)                                       ccGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcg                                                                                                >  2:285145/1‑139 (MQ=255)                                               aCCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGttt                                                                                        <  1:530689/139‑1 (MQ=255)                                                ccGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACatgatg                                                                                                                         <  1:552035/105‑1 (MQ=255)                                                ccGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACatgatg                                                                                                                         >  2:552035/1‑105 (MQ=255)                                                    ttGTATTCTGCTGGCATCAGGCATCTCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTcc                                                                                                                 <  2:33727/109‑1 (MQ=255)                                                    ttGTATTCTGCTGGCATCAGGCATCTCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTcc                                                                                                                 >  1:33727/1‑109 (MQ=255)                                                    ttGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGgaa                                                                                   >  2:1315867/1‑139 (MQ=255)                                                         ttCTGCTGGCATCAGGCATCGAAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACg                                                                              >  2:369687/1‑139 (MQ=255)                                                              ctgGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCg                                                                                                                                  >  2:252499/1‑82 (MQ=255)                                                              ctgGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCg                                                                                                                                  <  1:252499/82‑1 (MQ=255)                                                              ctgGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCg                                                                         <  1:1315867/139‑1 (MQ=255)                                                                   aTCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAg                                                                    >  1:356133/1‑139 (MQ=255)                                                                                         aCCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCg                                                                         <  2:365919/112‑1 (MQ=255)                                                                                         aCCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCg                                                                         >  1:365919/1‑112 (MQ=255)                                                                                         aCCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACt                                              >  1:511155/1‑139 (MQ=255)                                                                                            cGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCa                                           >  1:103720/1‑139 (MQ=255)                                                                                                      acaATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTctgct                                 <  2:118169/139‑1 (MQ=255)                                                                                                         aTGGCCTGGTAGAAGCGCTGAAAAATATCACCAAt                                                                                                                                      <  2:540025/35‑1 (MQ=255)                                                                                                         aTGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCtctc                                                  >  1:309137/1‑119 (MQ=255)                                                                                                         aTGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCtctc                                                  <  2:309137/119‑1 (MQ=255)                                                                                                         aTGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGa                              <  2:1071078/139‑1 (MQ=255)                                                                                                         aTGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGa                              >  1:1401861/1‑139 (MQ=255)                                                                                                         aTGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGa                              <  2:1401861/139‑1 (MQ=255)                                                                                                          tGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACatgatg                                                                                                                         >  1:1253319/1‑47 (MQ=255)                                                                                                          tGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACatgatg                                                                                                                         <  2:1253319/47‑1 (MQ=255)                                                                                                             ccTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg                                                                                             >  1:649060/1‑72 (MQ=255)                                                                                                             ccTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg                                                                                             <  2:649060/72‑1 (MQ=255)                                                                                                             ccTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATg                            >  1:346928/1‑137 (MQ=255)                                                                                                             ccTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATg                            <  2:346928/137‑1 (MQ=255)                                                                                                             ccTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATg                            >  2:640952/1‑137 (MQ=255)                                                                                                             ccTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATg                            <  1:640952/137‑1 (MQ=255)                                                                                                             ccTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGaa                          >  1:96940/1‑139 (MQ=255)                                                                                                             ccTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGaa                          >  1:1234396/1‑139 (MQ=255)                                                                                                               tGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGGTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAAt                        >  1:1028442/1‑139 (MQ=255)                                                                                                                      aGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGt                                      <  2:179996/118‑1 (MQ=255)                                                                                                                      aGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGt                                      >  1:179996/1‑118 (MQ=255)                                                                                                                      aGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAagag                 >  1:281366/1‑139 (MQ=255)                                                                                                                         gcTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTgg              <  2:356133/139‑1 (MQ=255)                                                                                                                         gcTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTgg              <  2:96940/139‑1 (MQ=255)                                                                                                                           tGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAg                                                                    <  2:1291492/83‑1 (MQ=255)                                                                                                                           tGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAg                                                                    >  1:1291492/1‑83 (MQ=255)                                                                                                                           tGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCa            <  2:103720/139‑1 (MQ=255)                                                                                                                               aaaTATCACCAATGTCGTACATGATGTGCGGTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCAGCtt        <  2:686060/139‑1 (MQ=255)                                                                                                                                     cACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCAGCTTCGTTTa  <  1:606715/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       CGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCAGCTTCGTTTA  >  NC_000913/2742826‑2743096

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