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%	180.5 / NA	56	D194N (GAT→AAT)	yfiN	putative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/29);  total (27/29)

GACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCAGCTTC  >  NC_000913/2742827‑2743091                                                                                                                                        |                                                                                                                                  gACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATg                                                                                                                                >  2:774052/1‑139 (MQ=255) gACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATg                                                                                                                                >  1:1292189/1‑139 (MQ=255)  aCAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGt                                                                                                                               <  1:662133/139‑1 (MQ=255)               ccATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGAtgtg                                                                                                                  >  2:806941/1‑139 (MQ=255)                cATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGc                                                                                                                 >  1:602493/1‑139 (MQ=255)                        cTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAAc                                                                                                         <  2:480099/139‑1 (MQ=255)                           gTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGc                                                                                                      >  1:187275/1‑139 (MQ=255)                              ttCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAAt                                                                                                   >  1:1075258/1‑139 (MQ=255)                                   tcgcCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTc                                                                                              <  1:774052/139‑1 (MQ=255)                                         tACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg                                                                                        <  1:1256285/139‑1 (MQ=255)                                         tACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg                                                                                        <  2:1292189/139‑1 (MQ=255)                                               ccGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTc                                                                                  >  1:677678/1‑139 (MQ=255)                                                cGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACGCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcga                                                                                          >  2:892889/1‑130 (MQ=255)                                                cGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcga                                                                                          <  1:892889/130‑1 (MQ=255)                                                 ggTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCgg                                                                                <  2:1279293/139‑1 (MQ=255)                                                   ttGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGgaa                                                                              <  1:1263558/139‑1 (MQ=255)                                                        ttCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACg                                                                         <  2:602493/139‑1 (MQ=255)                                                             ctgGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg                                                                                        <  2:85513/119‑1 (MQ=255)                                                             ctgGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg                                                                                        >  1:85513/1‑119 (MQ=255)                                                             ctgGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCg                                                                    >  1:279999/1‑139 (MQ=255)                                                             ctgGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCg                                                                    >  2:1170992/1‑139 (MQ=255)                                                               ggCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCt                                                                  <  2:279999/139‑1 (MQ=255)                                                               ggCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCt                                                                  <  2:187275/139‑1 (MQ=255)                                                                gCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg                                                                                        >  1:119678/1‑116 (MQ=255)                                                                gCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg                                                                                        <  2:119678/116‑1 (MQ=255)                                                                  aTCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg                                                                                        <  1:181670/114‑1 (MQ=255)                                                                  aTCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg                                                                                        >  2:181670/1‑114 (MQ=255)                                                                     aGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACatgatg                                                                                                                    <  1:66132/83‑1 (MQ=255)                                                                     aGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACatgatg                                                                                                                    >  2:66132/1‑83 (MQ=255)                                                                          tCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCg                                                       >  2:134380/1‑139 (MQ=255)                                                                           cGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGc                                                      >  1:43812/1‑139 (MQ=255)                                                                                 tACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGct                                                <  2:1075258/139‑1 (MQ=255)                                                                                                tttGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGt                                 <  1:1017184/139‑1 (MQ=255)                                                                                                   gCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTctg                              <  1:1204084/139‑1 (MQ=255)                                                                                                    cacaATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTctgc                             >  2:1208170/1‑139 (MQ=255)                                                                                                    cacaATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTctgc                             >  2:437445/1‑139 (MQ=255)                                                                                                       aaTGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGgaaga                                                                            <  1:1260403/89‑1 (MQ=255)                                                                                                       aaTGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGgaaga                                                                            >  2:1260403/1‑89 (MQ=255)                                                                                                        aTGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGgaa                                                                              >  1:1036334/1‑86 (MQ=255)                                                                                                        aTGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGgaa                                                                              <  2:1036334/86‑1 (MQ=255)                                                                                                            ccTGGTAGTAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGaa                     >  1:117141/1‑139 (MQ=255)                                                                                                            ccTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGa                      >  2:378408/1‑138 (MQ=255)                                                                                                            ccTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGa                      <  1:378408/138‑1 (MQ=255)                                                                                                              tGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGa                         >  1:286054/1‑133 (MQ=255)                                                                                                              tGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGa                         <  2:286054/133‑1 (MQ=255)                                                                                                              tGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAAt                   >  1:512054/1‑139 (MQ=255)                                                                                                                     aGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAagag            >  1:144805/1‑139 (MQ=255)                                                                                                                       cgcTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTg          >  1:152738/1‑139 (MQ=255)                                                                                                                         cTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGc        <  1:437445/139‑1 (MQ=255)                                                                                                                         cTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGc        <  1:1208170/139‑1 (MQ=255)                                                                                                                          tGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcga                                                                                          <  2:141775/56‑1 (MQ=255)                                                                                                                          tGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcga                                                                                          >  1:141775/1‑56 (MQ=255)                                                                                                                          tGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCa       <  2:144805/139‑1 (MQ=255)                                                                                                                          tGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCa       <  1:134380/139‑1 (MQ=255)                                                                                                                             aaaaTATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCAGCt    <  2:152738/139‑1 (MQ=255)                                                                                                                               aaTATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCAGCTTc  <  2:833487/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  GACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCAGCTTC  >  NC_000913/2742827‑2743091

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