Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,626,944	C→T	L665L (CTG→CTA)	pdeF ←	cyclic di‑GMP phosphodiesterase PdeF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,626,944	0	C	T	100.0%	116.3 / NA	34	L665L (CTG→CTA)	pdeF	cyclic di‑GMP phosphodiesterase PdeF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/18);  total (16/18)

GATCTCTTCGTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTAC  >  NC_000913/2626809‑2627080                                                                                                                                        |                                                                                                                                         gATCTCTTCGTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGAt                                                                                                                                      <  2:1885757/140‑1 (MQ=255)            tttCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTcaca                                                                                                                           >  1:2038244/1‑140 (MQ=255)                    tACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAg                                                                                                                   <  2:2222087/140‑1 (MQ=255)                    tACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAg                                                                                                                   <  1:2200288/140‑1 (MQ=255)                     aCTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGc                                                                                                                  >  1:2113468/1‑140 (MQ=255)                         tGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCcgcg                                                                                                              <  2:1857597/140‑1 (MQ=255)                                    gCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCg                                                                                                   <  2:1718024/140‑1 (MQ=255)                                     cATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGt                                                                                                  >  1:1596081/1‑140 (MQ=255)                                     cATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGt                                                                                                  >  1:2264019/1‑140 (MQ=255)                                                           tGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAAt                                                                            <  2:1688047/140‑1 (MQ=255)                                                               aaaTCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCa                                                                        >  2:2220273/1‑140 (MQ=255)                                                               aaaTCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCa                                                                        >  2:1701507/1‑140 (MQ=255)                                                               aaaTCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCa                                                                        <  1:1796387/140‑1 (MQ=255)                                                                     aTGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTgg                                                                  >  2:1800092/1‑140 (MQ=255)                                                                                  aTAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCt                                                     >  2:1603897/1‑140 (MQ=255)                                                                                  aTAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCt                                                     <  2:1878206/140‑1 (MQ=255)                                                                                      tCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGc                                                                 <  2:1669037/124‑1 (MQ=255)                                                                                      tCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGc                                                                 >  1:1669037/1‑124 (MQ=255)                                                                                            cTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGAtgctg                                           >  1:1926123/1‑140 (MQ=255)                                                                                             tATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGAtgctgc                                          <  1:1800092/140‑1 (MQ=255)                                                                                               tttGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCaa                                        >  2:1794599/1‑140 (MQ=255)                                                                                                ttGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCGTCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAg                                       >  1:1668122/1‑140 (MQ=255)                                                                                                 tGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAgg                                      <  1:2220273/140‑1 (MQ=255)                                                                                                         tatTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTgc                              >  1:2217969/1‑140 (MQ=255)                                                                                                           ttGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTgcgc                            <  2:1563192/140‑1 (MQ=255)                                                                                                           ttGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTgcgc                            <  2:2264019/140‑1 (MQ=255)                                                                                                           ttGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTgcgc                            <  1:1603897/140‑1 (MQ=255)                                                                                                            tGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGcc                           >  1:1597503/1‑140 (MQ=255)                                                                                                               ggATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGA‑TTGCGCCTGCt                       <  2:2038244/140‑1 (MQ=255)                                                                                                                gATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCt                       <  1:1701507/140‑1 (MQ=255)                                                                                                                   aaaaCTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTaa                    <  2:1992672/140‑1 (MQ=255)                                                                                                                   aaaaCTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTaa                    >  1:1784870/1‑140 (MQ=255)                                                                                                                   aaaaCTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTaa                    >  2:2276055/1‑140 (MQ=255)                                                                                                                                     aaGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTAc  <  2:2113468/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                         GATCTCTTCGTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTAC  >  NC_000913/2626809‑2627080

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