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

GACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCAGCTTCGTTTA  >  NC_000913/2742827‑2743096                                                                                                                                        |                                                                                                                                       gACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATg                                                                                                                                     <  1:30178/139‑1 (MQ=255)    aGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCg                                                                                                                                  <  2:1049485/139‑1 (MQ=255)               ccATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGAtgtg                                                                                                                       >  2:520718/1‑139 (MQ=255)               ccATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGAtgtg                                                                                                                       >  2:211497/1‑139 (MQ=255)                         tGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAAcc                                                                                                             >  2:119346/1‑139 (MQ=255)                              ttCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTcc                                                                                                                 >  2:607290/1‑130 (MQ=255)                              ttCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTcc                                                                                                                 <  1:607290/130‑1 (MQ=255)                               tcgctcgcCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTAc                                                                                                                               <  2:514681/115‑1 (MQ=255)                               tcgctcgcCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTAc                                                                                                                               >  1:514681/1‑115 (MQ=255)                                 gctcgcCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGAtgt                                                                                                                        <  2:1405523/120‑1 (MQ=255)                                 gctcgcCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGAtgt                                                                                                                        >  1:1405523/1‑120 (MQ=255)                                   tcgcCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTc                                                                                                   <  2:449330/139‑1 (MQ=255)                                      ccGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcg                                                                                                <  2:765022/139‑1 (MQ=255)                                       cGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcga                                                                                               <  2:275995/139‑1 (MQ=255)                                         tACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg                                                                                             <  2:493293/139‑1 (MQ=255)                                                cGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGttt                                                                                        <  2:121102/137‑1 (MQ=255)                                                cGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGttt                                                                                        >  1:121102/1‑137 (MQ=255)                                                   ttGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGttt                                                                                        <  1:1226519/134‑1 (MQ=255)                                                   ttGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGttt                                                                                        >  2:1226519/1‑134 (MQ=255)                                                                gCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTg                                                                      >  1:1526838/1‑139 (MQ=255)                                                                       gCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAAttt                                                                                                      >  2:811312/1‑100 (MQ=255)                                                                       gCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAAttt                                                                                                      <  1:811312/100‑1 (MQ=255)                                                                               aTTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCg                                                       <  2:1505578/139‑1 (MQ=255)                                                                                ttACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTcc                                                                                                                 >  1:1164947/1‑80 (MQ=255)                                                                                ttACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTcc                                                                                                                 <  2:1164947/80‑1 (MQ=255)                                                                                  aCCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGctc                                                    <  1:211497/139‑1 (MQ=255)                                                                                       cACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGAc                                               >  2:591208/1‑139 (MQ=255)                                                                                        aCCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACatgatg                                                                                                                         >  1:1327305/1‑64 (MQ=255)                                                                                        aCCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACatgatg                                                                                                                         <  2:1327305/64‑1 (MQ=255)                                                                                            gCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCaa                                          <  2:128442/139‑1 (MQ=255)                                                                                              cATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACa                                        <  1:591208/139‑1 (MQ=255)                                                                                                       aaTGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCg                                                                         <  2:1400621/97‑1 (MQ=255)                                                                                                       aaTGGCCTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCg                                                                         >  1:1400621/1‑97 (MQ=255)                                                                                                            ccTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTcc                                                                                                                 <  2:238463/52‑1 (MQ=255)                                                                                                            ccTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTcc                                                                                                                 >  1:238463/1‑52 (MQ=255)                                                                                                             cTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGaaa                         <  2:1526838/139‑1 (MQ=255)                                                                                                             cTGGTAGAAGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGaaa                         <  1:119346/139‑1 (MQ=255)                                                                                                                     aGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACtt                                             <  1:664923/111‑1 (MQ=255)                                                                                                                     aGCGCTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACtt                                             >  2:664923/1‑111 (MQ=255)                                                                                                                      gcgcTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg                                                                                             <  2:791023/62‑1 (MQ=255)                                                                                                                      gcgcTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg                                                                                             >  1:791023/1‑62 (MQ=255)                                                                                                                        gcTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTCTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTgg              <  1:520718/139‑1 (MQ=255)                                                                                                                         cTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg                                                                                             >  1:729193/1‑59 (MQ=255)                                                                                                                         cTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCcgacg                                                                                             <  2:729193/59‑1 (MQ=255)                                                                                                                         cTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACg                                                                              >  2:806591/1‑74 (MQ=255)                                                                                                                         cTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACg                                                                              <  1:806591/74‑1 (MQ=255)                                                                                                                         cTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACa                                        >  1:1661167/1‑112 (MQ=255)                                                                                                                         cTGAAAAATATCACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACa                                        <  2:1661167/112‑1 (MQ=255)                                                                                                                                    cACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCAGCTTCGTTTa  >  1:615437/1‑139 (MQ=255)                                                                                                                                    cACCAATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCAGCTTCGTTTa  <  1:905145/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                       GACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAGCGCTGAAAAATATCACCGATGTCGTACATGATGTGCGTTCCAACCGCAATTTTTCCCGACGAGTTTCGGAAGAACGTATCGCTGAGTTTCACCGCTTCGCTCTCGACTTCAACAGTCTGCTGGATGAAATGGAAGAGTGGCAGCTTCGTTTA  >  NC_000913/2742827‑2743096

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