Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,684,903	T→C	N239S (AAC→AGC)	pdeK ←	putative c‑di‑GMP phosphodiesterase PdeK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,684,903	0	T	C	100.0%	172.0 / NA	50	N239S (AAC→AGC)	pdeK	putative c‑di‑GMP phosphodiesterase PdeK
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/23);  total (27/23)

GTATCACGCAGGGTTTCACAGGTGATGATCATCAGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGTTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGAT  >  NC_000913/3684771‑3685033                                                                                                                                     |                                                                                                                                   gTATCACGCAGGGTTTCTCAGGCGATGATCAACAGTGCGTTGGTTTGTGTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCgctctgttc                                                                                                                             <  2:236093/140‑1 (MQ=255)   aTCACGCAGGGTTTCACAGGTGATGATCATCAGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGttctt                                                                                                                           <  2:135106/140‑1 (MQ=255)     cACGCAGGGTTTCACAGGTGATGATCATCAGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCa                                                                                                                         >  2:31228/1‑140 (MQ=255)       cGCAGGGTTTCACAGGTGATGATCATCAGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATa                                                                                                                       >  2:480740/1‑140 (MQ=255)         cAGGGTTTCACAGGTGATGATCATCAGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAAt                                                                                                                     >  2:406628/1‑140 (MQ=255)          aGGGTTTCACAGGTGATGATCATCAGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATg                                                                                                                    >  1:476982/1‑140 (MQ=255)            ggTTTCACAGGTGATGATCATCAGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGgc                                                                                                                  >  1:152735/1‑140 (MQ=255)            ggTTTCACAGGTGATGATCATCAGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGCAAGCGCAACGCGTTGTCGCTCAGTTCTTCATAATGgc                                                                                                                  >  1:170075/1‑140 (MQ=255)            ggTTTCACAGGTGATCATCATCTGCACGGCGGTTTGTATACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACGGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGgc                                                                                                                  <  1:212444/140‑1 (MQ=255)                      gTGATGATCATCAGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCaa                                                                                                        >  2:555762/1‑140 (MQ=255)                         atgatCATCAGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTg                                                                                                     <  2:537968/140‑1 (MQ=255)                             tcatcaGCGCGGTGGTTTGTTTTCGCGCGACAACCTGCTACAGCATCTGCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTgg                                                                                                 <  2:226305/140‑1 (MQ=255)                              catcaGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGt                                                                                                >  1:600293/1‑140 (MQ=255)                               atcaGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGtt                                                                                               <  1:31228/140‑1 (MQ=255)                               atcaGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGtt                                                                                               <  1:406628/140‑1 (MQ=255)                               atcaGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGtt                                                                                               <  1:529617/140‑1 (MQ=255)                               atcaGCGCGGTGGTTCGTTTACTCGGGACAACCTGCTACAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGtt                                                                                               <  1:227167/140‑1 (MQ=255)                                tcaGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGgcgc                                                                                                                >  1:634797/1‑122 (MQ=255)                                tcaGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGgcgc                                                                                                                <  2:634797/122‑1 (MQ=255)                                tcaGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTg                                                                                              >  2:399239/1‑140 (MQ=255)                                tcaGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCGCTGTTCTACATGATGGCGCTGGAGCAATTGCTGGTTg                                                                                              >  1:170504/1‑140 (MQ=255)                                 caGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGTATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGa                                                                                             >  2:487253/1‑140 (MQ=255)                                  aGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAg                                                                                            >  2:429934/1‑140 (MQ=255)                                   gcgcGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTTTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAgg                                                                                           >  1:303200/1‑140 (MQ=255)                                    cgcgGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGt                                                                                          >  2:139766/1‑140 (MQ=255)                                    cgcgGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGGTCTTCATAATGACGGTGCAGCAATTGCTGGTTGAGGt                                                                                          >  1:204694/1‑140 (MQ=255)                                         tggtTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTaa                                                                                     <  1:344995/140‑1 (MQ=255)                                          ggtTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAAc                                                                                    <  1:139766/140‑1 (MQ=255)                                              tgtttACGCGCGACAACCTGCTCCAGTATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTgcg                                                                                <  1:487253/140‑1 (MQ=255)                                                 ttACGCGCGACAACCTGCTCCACCAACTCGATCAGCACGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCAc                                                                             <  1:197655/140‑1 (MQ=255)                                                     gcgcgACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAAc                                                                         >  2:365078/1‑140 (MQ=255)                                                                 gCTCCAGCATCTCCATCCGCAATGCTTTGTTCGGCCAATTCGACACCGTGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCa                                                             <  2:170075/140‑1 (MQ=255)                                                                     cAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGACTTCATCGt                                                         >  1:335739/1‑140 (MQ=255)                                                                        cATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCACATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCt                                                      >  2:420975/1‑140 (MQ=255)                                                                         aTCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTg                                                     >  2:402092/1‑140 (MQ=255)                                                                            tcCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATg                                                  >  1:27050/1‑140 (MQ=255)                                                                                tCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGa                                              <  1:480740/140‑1 (MQ=255)                                                                                tCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGa                                              <  2:165248/140‑1 (MQ=255)                                                                                  aGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGAcg                                            >  2:638712/1‑140 (MQ=255)                                                                                   gCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGAcgc                                           <  2:465272/140‑1 (MQ=255)                                                                                     aaGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCgg                                         <  1:638712/140‑1 (MQ=255)                                                                                        gCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAAAATACCGATTTCAACGTCCTGATGCAGACGCGGTaa                                      >  2:296894/1‑140 (MQ=255)                                                                                         cTTTGTTCGGCAAATCAGACACCGGGAAGCGGATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAAt                                     <  2:303200/140‑1 (MQ=255)                                                                                            tGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGcc                                  >  1:68310/1‑140 (MQ=255)                                                                                                 ggCAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCATCAACATACGGATGTCATCGTCCTGAAGCAGGCGCGGTAATGCCAGTTg                             >  2:194370/1‑140 (MQ=255)                                                                                                   cAAATCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGt                           <  2:152735/140‑1 (MQ=255)                                                                                                     aaTCCGACACCGGGAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGa                         <  1:420975/140‑1 (MQ=255)                                                                                                                 ggAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTcc             >  2:27004/1‑140 (MQ=255)                                                                                                                 ggAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTcc             >  1:519953/1‑140 (MQ=255)                                                                                                                 ggAAGCGCATCGCATTTTCGCTCTGCTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTcc             <  1:319891/140‑1 (MQ=255)                                                                                                                  gAAGCGCATCGCATTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCt            >  2:430564/1‑140 (MQ=255)                                                                                                                            gcatTTTCGCTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGAt  <  2:519953/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                   GTATCACGCAGGGTTTCACAGGTGATGATCATCAGCGCGGTGGTTTGTTTACGCGCGACAACCTGCTCCAGCATCTCCATCAGCAAGGCTTTGTTCGGCAAATCCGACACCGGGAAGCGCATCGCATTTTCGTTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGAT  >  NC_000913/3684771‑3685033

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