Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,516,313	T→C	Y507C (TAC→TGC)	yfeA ←	putative diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,516,313	0	T	C	100.0%	170.9 / NA	46	Y507C (TAC→TGC)	yfeA	putative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/22);  total (24/22)

AAAACGCGCGCTAAGGTTAAACTGAGCAATAAGGGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTT  >  NC_000913/2516180‑2516459                                                                                                                                      |                                                                                                                                                   aaaaCGCGCGCTAAGGTTAAACTGAGCAATAAGGGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCAccc                                                                                                                                   >  1:582052/1‑151 (MQ=255)         cgcTAAGGTTAAACTGAGCAATAAGGGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTaa                                                                                                                           <  1:739251/151‑1 (MQ=255)          gcTAAGGTTAAACTGAGCAATAAGGGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCCGAGCAGCAAATCACCCTGATCTAAc                                                                                                                          >  2:742995/1‑151 (MQ=255)             aaGGTTAAACTGAGCAATAAGGGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGcc                                                                                                                       >  2:350870/1‑151 (MQ=255)             aaGGTTAAACTGAGCAATAAGGGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGcc                                                                                                                       <  2:236348/151‑1 (MQ=255)                   aaaCTGAGCAATAAGGGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCg                                                                                                                  <  1:312511/150‑1 (MQ=255)                   aaaCTGAGCAATAAGGGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCg                                                                                                                  <  2:842929/150‑1 (MQ=255)                       tGAGCAATAAGGGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATgg                                                                                                             >  2:744448/1‑151 (MQ=255)                        gAGCAATAAGGGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATg                                                                                                              <  1:178963/149‑1 (MQ=255)                           cAATAAGGGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTAtt                                                                                                         >  2:360390/1‑151 (MQ=255)                               aaGGGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAg                                                                                                      <  1:210857/150‑1 (MQ=255)                                  gggCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAAt                                                                                                  >  2:382609/1‑151 (MQ=255)                                    gCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAg                                                                                                <  2:307024/151‑1 (MQ=255)                                      aGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGc                                                                                               >  2:79254/1‑150 (MQ=255)                                         aaCTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCt                                                                                           >  2:54731/1‑151 (MQ=255)                                         aaCTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCt                                                                                           >  2:296262/1‑151 (MQ=255)                                         aaCTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCt                                                                                           <  1:728625/151‑1 (MQ=255)                                          aCTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTg                                                                                          <  2:312783/151‑1 (MQ=255)                                             tATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGttt                                                                                       >  1:56783/1‑151 (MQ=255)                                             tATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGt                                                                                         <  1:480160/149‑1 (MQ=255)                                                 cGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTg                                                                                   >  1:108517/1‑151 (MQ=255)                                                  ggggTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTgg                                                                                  >  1:612361/1‑151 (MQ=255)                                                      tCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTg                                                                              >  1:44287/1‑151 (MQ=255)                                                       cATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTg                                                                              <  1:268578/150‑1 (MQ=255)                                                        aTAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAACGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGcc                                                                            >  1:812406/1‑151 (MQ=255)                                                          aaTGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGccc                                                                           <  2:108495/150‑1 (MQ=255)                                                                gccgTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAAcc                                                                    >  1:180928/1‑151 (MQ=255)                                                                     cATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCAtctc                                                               >  2:272324/1‑151 (MQ=255)                                                                     cATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCAtctc                                                               >  1:383821/1‑151 (MQ=255)                                                                           tCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCctctct                                                         >  1:847378/1‑151 (MQ=255)                                                                           tCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCctctc                                                          <  2:39339/150‑1 (MQ=255)                                                                              gTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCt                                                      <  2:806147/151‑1 (MQ=255)                                                                                      aGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAgcg                                              <  1:504387/151‑1 (MQ=255)                                                                                         aTCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGcc                                            >  1:300687/1‑150 (MQ=255)                                                                                           ctcATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGc                                         >  2:435702/1‑151 (MQ=255)                                                                                           ctcATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGc                                         >  2:174577/1‑151 (MQ=255)                                                                                             catcatAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCa                                           <  1:762227/147‑1 (MQ=255)                                                                                             catcatAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAg                                          <  1:418154/148‑1 (MQ=255)                                                                                                   aaCCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCTCCAGCACatgat                                  <  2:161915/150‑1 (MQ=255)                                                                                                        tCACCCTCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGAtgtg                            >  1:904680/1‑151 (MQ=255)                                                                                                        tCACCCGCTTTGTTGCGAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACACCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGAtgt                             <  2:644735/150‑1 (MQ=255)                                                                                                                        gAATTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCgc            >  2:421133/1‑151 (MQ=255)                                                                                                                          aTTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCgcca          <  1:482978/151‑1 (MQ=255)                                                                                                                          aTTGGCTGGGCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCgcca          <  1:435749/151‑1 (MQ=255)                                                                                                                                 gggCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCcagccagc    <  1:906168/150‑1 (MQ=255)                                                                                                                                  ggCGCAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCtt  >  1:130991/1‑151 (MQ=255)                                                                                                                                      |                                                                                                                                                   AAAACGCGCGCTAAGGTTAAACTGAGCAATAAGGGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTT  >  NC_000913/2516180‑2516459

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