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

AATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCAGCTTCGTTTA  >  NC_000913/2742835‑2743096                                                                                                                                |                                                                                                                                       aaTCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACat                                                                                                                             <  1:1128857/139‑1 (MQ=255)  aTCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACatg                                                                                                                            <  1:630945/139‑1 (MQ=255)        cATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGc                                                                                                                      >  2:1308414/1‑139 (MQ=255)                        cgctcgcCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAAttt                                                                                                      >  1:1197900/1‑139 (MQ=255)                                      aCCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGttt                                                                                        <  2:1197900/139‑1 (MQ=255)                                           ttGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACatgatg                                                                                                                         <  1:721548/101‑1 (MQ=255)                                           ttGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACatgatg                                                                                                                         >  2:721548/1‑101 (MQ=255)                                                  ctgctgGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTa                                                                            <  1:1308414/139‑1 (MQ=255)                                                    gctgGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATc                                                                          >  1:1095382/1‑139 (MQ=255)                                                     ctgGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCg                                                                         >  1:501758/1‑139 (MQ=255)                                                     ctgGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCg                                                                         >  1:320923/1‑139 (MQ=255)                                                      tgGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGc                                                                        >  2:773084/1‑139 (MQ=255)                                                             aGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGttt                                                                 >  1:234813/1‑139 (MQ=255)                                                             aGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGttt                                                                 <  2:234813/139‑1 (MQ=255)                                                               gCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAAcc                                                                                                             >  1:122644/1‑93 (MQ=255)                                                               gCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAAcc                                                                                                             <  2:122644/93‑1 (MQ=255)                                                                  tCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCg                                                            >  2:440265/1‑139 (MQ=255)                                                                   cGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGgaaga                                                                                 >  2:127790/1‑117 (MQ=255)                                                                   cGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTACGgaaga                                                                                 <  1:127790/117‑1 (MQ=255)                                                                          aCCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGctc                                                    <  1:440265/139‑1 (MQ=255)                                                                             cTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGt                                                                                                                                    >  1:559027/1‑56 (MQ=255)                                                                             cTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGt                                                                                                                                    <  2:559027/56‑1 (MQ=255)                                                                                aCCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgac                                                                                              <  1:12378/91‑1 (MQ=255)                                                                                aCCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgac                                                                                              >  2:12378/1‑91 (MQ=255)                                                                                                 tGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGctc                                                    <  2:712054/116‑1 (MQ=255)                                                                                                 tGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGctc                                                    >  1:712054/1‑116 (MQ=255)                                                                                                      tGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGAc                                               >  2:861431/1‑116 (MQ=255)                                                                                                      tGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGAc                                               <  1:861431/116‑1 (MQ=255)                                                                                                         tAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGc                                                                        >  2:1583207/1‑88 (MQ=255)                                                                                                         tAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGc                                                                        <  1:1583207/88‑1 (MQ=255)                                                                                                         tAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATgg                      <  2:563209/138‑1 (MQ=255)                                                                                                         tAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATgg                      >  1:563209/1‑138 (MQ=255)                                                                                                           gAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCg                                                                         >  1:830882/1‑85 (MQ=255)                                                                                                           gAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCg                                                                         <  2:830882/85‑1 (MQ=255)                                                                                                             aGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAagag                 >  2:149072/1‑139 (MQ=255)                                                                                                                  tGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGa                     <  2:283012/130‑1 (MQ=255)                                                                                                                  tGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGa                     >  1:283012/1‑130 (MQ=255)                                                                                                                  tGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCa            >  1:991127/1‑139 (MQ=255)                                                                                                                  tGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCa            >  1:925138/1‑139 (MQ=255)                                                                                                                            cACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCAGCTTCGTTTa  <  1:149072/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                       AATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCAGCTTCGTTTA  >  NC_000913/2742835‑2743096

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