Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,684,604	G→A	Q339* (CAG→TAG)	yhjK ←	cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,684,604	0	G	A	100.0%	207.8 / NA	59	Q339* (CAG→TAG)	yhjK	cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/26);  total (33/26)

GATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGC  >  NC_000913/3684508‑3684695                                                                                                 |                                                                                            gaTCGCCGTCGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCTGGCGCTCGCTCATGATAGTGAGCACTTGCTAAc                                                                                           <  1:1118819/99‑1 (MQ=255) gaTCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAAcc                                                                                          >  2:952289/1‑100 (MQ=255) gaTCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAAcc                                                                                          <  2:417527/100‑1 (MQ=255)   tCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTa                                                                                        <  1:77152/100‑1 (MQ=255)    ngCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTaa                                                                                       >  1:1328746/2‑100 (MQ=255)          aGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGAt                                                                                 >  1:949927/1‑100 (MQ=255)            aaCATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTg                                                                               >  2:671080/1‑100 (MQ=255)             aCATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGc                                                                              <  2:278453/100‑1 (MQ=255)             aCATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGc                                                                              >  2:242163/1‑100 (MQ=255)             aCATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGc                                                                              >  1:563361/1‑100 (MQ=255)              cATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCg                                                                             >  1:574818/1‑100 (MQ=255)              cATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCg                                                                             <  1:473213/100‑1 (MQ=255)                tCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTg                                                                           <  1:952285/100‑1 (MQ=255)                      cgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACgg                                                                     >  1:80801/1‑100 (MQ=255)                         caATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTc                                                                  >  1:863133/1‑100 (MQ=255)                          aaTGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTcc                                                                 >  2:1072149/1‑100 (MQ=255)                            tGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTg                                                               <  1:169159/100‑1 (MQ=255)                            tGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTg                                                               >  2:1159508/1‑100 (MQ=255)                             gCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTg                                                               <  2:514061/99‑1 (MQ=255)                                     gCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGtt                                                      >  1:659093/1‑100 (MQ=255)                                        gacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGc                                                   <  1:1521293/100‑1 (MQ=255)                                        gacggaCTTGAATCCGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGttc                                                     <  1:1015253/98‑2 (MQ=255)                                         acggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCa                                                  <  2:718594/100‑1 (MQ=255)                                           ggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAAt                                                <  2:244508/100‑1 (MQ=255)                                             aGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACATCGTTGGCAATGa                                              >  2:69372/1‑100 (MQ=255)                                             aGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGa                                              >  1:556966/1‑100 (MQ=255)                                              gTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGAc                                             <  1:1300030/100‑1 (MQ=255)                                               ttGAATACGTTCAATCGGAAGGCGCTCGCTTATGATAGTGAGTACTTGCTAACCTATTGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACa                                            <  1:1481438/100‑1 (MQ=255)                                                tGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAg                                           <  2:672518/100‑1 (MQ=255)                                                tGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAg                                           >  2:1260780/1‑100 (MQ=255)                                                tGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTAAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAg                                           >  1:1045821/1‑100 (MQ=255)                                                 gAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGc                                          >  2:1257104/1‑100 (MQ=255)                                                 gAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGCTTCCTGTACACCGTTGGCAATGACAGc                                          <  1:283827/100‑1 (MQ=255)                                                    tACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCaaa                                       >  2:656055/1‑100 (MQ=255)                                                       gTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTc                                    <  2:177446/100‑1 (MQ=255)                                                       gTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTc                                    >  2:25088/1‑100 (MQ=255)                                                       gTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTc                                    <  2:1342458/100‑1 (MQ=255)                                                        ttCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCa                                   <  2:516791/100‑1 (MQ=255)                                                          cAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATa                                 <  2:1302267/100‑1 (MQ=255)                                                            aTCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCGAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAAc                               >  1:221550/1‑100 (MQ=255)                                                               ggCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGc                            >  2:1566512/1‑100 (MQ=255)                                                               ggCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGc                            >  1:1201917/1‑100 (MQ=255)                                                               ggCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGc                            >  1:591565/1‑100 (MQ=255)                                                                   ggCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAAt                        >  1:963848/1‑100 (MQ=255)                                                                     cgctcgctCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCt                      >  1:151903/1‑100 (MQ=255)                                                                      gctcgctcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTg                     >  1:608908/1‑100 (MQ=255)                                                                       ctcgctcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTgc                    >  2:1540231/1‑100 (MQ=255)                                                                         cgctcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTgcgc                  >  2:230241/1‑100 (MQ=255)                                                                          gctcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAATGTCATAACCGCTAATCTGcgcg                 <  1:990922/100‑1 (MQ=255)                                                                           ctcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGcgcg                 <  1:1169791/99‑1 (MQ=255)                                                                            tcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAg               <  2:890032/100‑1 (MQ=255)                                                                            tcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAg               <  2:1211274/100‑1 (MQ=255)                                                                             cATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCCTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGATAATCTGCGCGAgg              >  2:1198459/1‑100 (MQ=255)                                                                               tGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGAt            >  2:643954/1‑100 (MQ=255)                                                                               tGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGAt            >  1:690374/1‑100 (MQ=255)                                                                                   aGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCAt         >  1:1362658/1‑99 (MQ=255)                                                                                      gAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGt     >  1:860555/1‑100 (MQ=255)                                                                                        gCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTgg   <  2:258127/100‑1 (MQ=255)                                                                                         cACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGc  <  2:949930/100‑1 (MQ=255)                                                                                                 |                                                                                            GATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGC  >  NC_000913/3684508‑3684695

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑