Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,027,398	C→T	W207* (TGG→TAG)	yedQ ←	putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,027,398	0	C	T	100.0%	119.5 / NA	35	W207* (TGG→TAG)	yedQ	putative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/18);  total (17/18)

CATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAG  >  NC_000913/2027264‑2027517                                                                                                                                       |                                                                                                                        cATACCAGTAGTTATTACTGTATACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAggg                                                                                                                     <  1:231403/139‑1 (MQ=255)           gtTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACa                                                                                                           >  2:866989/1‑139 (MQ=255)            ttattaCTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAATCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACAt                                                                                                          <  2:363746/139‑1 (MQ=255)              attaCTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAg                                                                                                        <  1:866989/139‑1 (MQ=255)                taCTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGcc                                                                                                      >  2:80535/1‑139 (MQ=255)                     tCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACAt                                                                                                          <  1:159023/130‑1 (MQ=255)                     tCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACAt                                                                                                          >  2:159023/1‑130 (MQ=255)                     tCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGtaat                                                                                                 >  2:449353/1‑139 (MQ=255)                         ccGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAAcg                                                                                             >  2:308497/1‑139 (MQ=255)                         ccGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAAcg                                                                                             >  2:10951/1‑139 (MQ=255)                                      ggTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTg                                                                                <  2:353799/139‑1 (MQ=255)                                          aCCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGcgcg                                                                            >  2:659420/1‑139 (MQ=255)                                                            cTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAggg                                                                                                                     >  2:3734/1‑80 (MQ=255)                                                            cTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAggg                                                                                                                     <  1:3734/80‑1 (MQ=255)                                                            cTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCg                                                          >  1:581741/1‑139 (MQ=255)                                                                 gtgtTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGaa                                                     >  1:546026/1‑139 (MQ=255)                                                                    ttCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGAc                                                                                                            >  1:300933/1‑81 (MQ=255)                                                                    ttCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGAc                                                                                                            <  2:300933/81‑1 (MQ=255)                                                                             cGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCg                                         <  1:263451/139‑1 (MQ=255)                                                                              gAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAAcgcg                                                                                           <  1:633432/88‑1 (MQ=255)                                                                              gAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAAcgcg                                                                                           >  2:633432/1‑88 (MQ=255)                                                                                 gTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCAc                                     <  1:274750/139‑1 (MQ=255)                                                                                   gAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCAcgc                                   <  1:659420/139‑1 (MQ=255)                                                                                       cAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCgaga                               <  2:546026/139‑1 (MQ=255)                                                                                       cAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCgaga                               <  1:10951/139‑1 (MQ=255)                                                                                                 cgcgGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAAcgcg                                                                                           <  2:877960/69‑1 (MQ=255)                                                                                                 cgcgGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAAcgcg                                                                                           >  1:877960/1‑69 (MQ=255)                                                                                                   cgGTGACGATTTTCTCGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCt                   >  2:843711/1‑139 (MQ=255)                                                                                                                ctcGTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATcga      >  1:144103/1‑139 (MQ=255)                                                                                                                   gTTGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATcgacga   <  2:581741/139‑1 (MQ=255)                                                                                                                    ttGTGAATGGCCAATAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACgag  <  1:80535/139‑1 (MQ=255)                                                                                                                              ccAATAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATc                       >  2:3047/1‑108 (MQ=255)                                                                                                                              ccAATAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATc                       <  1:3047/108‑1 (MQ=255)                                                                                                                                 aTAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACAtt                                                                                 <  2:200392/47‑1 (MQ=255)                                                                                                                                 aTAAACTAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACAtt                                                                                 >  1:200392/1‑47 (MQ=255)                                                                                                                                       |                                                                                                                        CATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAG  >  NC_000913/2027264‑2027517

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: ‑