Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,345,482	G→A	Q421* (CAG→TAG)	gmr ←	cyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of RNase II stability

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,345,482	0	G	A	100.0%	124.2 / NA	32	Q421* (CAG→TAG)	gmr	cyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of RNase II stability
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/18);  total (14/18)

GCCCTGACTCTTCGGCGTAGGAAATGAAGTCCAACGGTGGAATCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTGATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTGGGGTAAAAACGCAAAATTGTCCTCGTCCGCCTTCCTTCGCTGTGTACATTGCTGTGTCGGCGTGACGAATAAT  >  NC_000913/1345344‑1345626                                                                                                                                           |                                                                                                                                                 gCCCTGACTCTTCGGCGTAGGAAATGAAGTCCAACGGTGGAATCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAg                                                                                                                                      <  1:320739/151‑1 (MQ=255)  cccTGACTCTTCGGCGTAGGAAATGAAGTCCAACGGTGGAATCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGt                                                                                                                                     <  2:430672/151‑1 (MQ=255)            tCGGCGTAGGAAATGAAGTCCAACGGTGGAATCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTAc                                                                                                                              <  1:82644/148‑1 (MQ=255)            tCGGCGTAGGAAATGAAGTCCAACGGTGGAATCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGc                                                                                                                            >  1:524242/1‑150 (MQ=255)                       aaTGAAGTCCAACGGTGGAATCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATcc                                                                                                                <  2:526163/151‑1 (MQ=255)                        aTGAAGTCCAACGGTGGAATCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATcca                                                                                                               <  1:423144/151‑1 (MQ=255)                             gTCCAACGGTGGAATCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATccagccag                                                                                                          <  2:178312/151‑1 (MQ=255)                                 aaCGGTGGAATCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATccagccag                                                                                                          >  1:17067/1‑147 (MQ=255)                                 aaCGGTGGAATCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATccagccag                                                                                                          >  2:151386/1‑147 (MQ=255)                                        gAATCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCaaa                                                                                               <  2:410943/151‑1 (MQ=255)                                        gAATCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCaa                                                                                                >  1:58590/1‑150 (MQ=255)                                          aTCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAAt                                                                                              >  2:539130/1‑150 (MQ=255)                                          aTCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATa                                                                                             <  2:624625/151‑1 (MQ=255)                                           tCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATAc                                                                                            >  1:236345/1‑151 (MQ=255)                                               cccACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGt                                                                                        <  1:343056/151‑1 (MQ=255)                                                  aCGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTg                                                                                      >  1:310708/1‑150 (MQ=255)                                                  aCGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGa                                                                                     <  1:86803/151‑1 (MQ=255)                                                   cGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGAt                                                                                    <  2:152746/151‑1 (MQ=255)                                                            gACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTg                                                                           <  1:592486/151‑1 (MQ=255)                                                            gACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTg                                                                           <  2:185201/151‑1 (MQ=255)                                                                              gCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTGGGGTAAAAACGCAAAAtt                                                         >  1:618556/1‑151 (MQ=255)                                                                                ttCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTGGGGTAAAAACGCAAAATTGt                                                       >  1:22301/1‑151 (MQ=255)                                                                                 tCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTGGGGTAAAAACGCAAAATTg                                                        <  2:536862/149‑1 (MQ=255)                                                                                    aGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTGGGGTAAAAACGCAAAATTGTCCTc                                                   >  2:423685/1‑151 (MQ=255)                                                                                      aCTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTGGGGTAAAAACGCAAAATTGTCCTCGt                                                 <  2:387176/151‑1 (MQ=255)                                                                                         gcgcACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTGGGGTAAAAACGCAAAATTGTCCTCGTcc                                               <  1:380842/150‑1 (MQ=255)                                                                                         gcgcACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTGGGGTAAAAACGCAAAATTGTCCTCGTCCg                                              >  2:284145/1‑151 (MQ=255)                                                                                                             tGATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTGGGGTAAAAACGCAAAATTGTCCTCGTCCGCCTTCCTTCGCTGTGTACAt                          >  1:160273/1‑151 (MQ=255)                                                                                                              gATTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTGGGGTAAAAACGCAAAATTGTCCTCGTCCGCCTTCCTTCGCTGTGTACAtt                         <  2:236340/151‑1 (MQ=255)                                                                                                               aTTTTCGGTTGATAGTGAATAACCAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTGGGGTAAAAACGCAAAATTGTCCTCGTCCGCCTTCCTTCGCTGTGTACATTg                        >  2:551824/1‑151 (MQ=255)                                                                                                                                     ccAACTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTGGGGTAAAAACGCAAAATTGTCCTCGTCCGCCTTCCTTCGCTGTGTACATTGCTGTGTCGGCGTGACGaataat  <  1:436846/151‑1 (MQ=255)                                                                                                                                       aaCTAATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTGGGGTAAAAACGCAAAATTGTCCTCGTCCGCCTTCCTTCGCTGTGTACATTGCTGTGTCGGCGTGACGaataat  >  1:99967/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                                 GCCCTGACTCTTCGGCGTAGGAAATGAAGTCCAACGGTGGAATCAACCCACGTTCAGGTGACTGCCAACGTACTAGTGCTTCCAGACTGCGCACTTCGCCACGCCAGGTGATTTTCGGTTGATAGTGAATAACCAACTGATCGTTTTCCAGTGCTTTACGCAAGTTGGTATCCAGCCAGAGATATTCAAATACCCGTTGATTCATTTCTGGGGTAAAAACGCAAAATTGTCCTCGTCCGCCTTCCTTCGCTGTGTACATTGCTGTGTCGGCGTGACGAATAAT  >  NC_000913/1345344‑1345626

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