Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,683,630	C→T	*663* (TAG→TAA)	yhjK ←	cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,683,630	0	C	T	100.0%	73.7 / NA	24	*663* (TAG→TAA)	yhjK	cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/12);  total (12/12)

CTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTG  >  NC_000913/3683502‑3683739                                                                                                                                 |                                                                                                              cTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                     <  2:678590/139‑1 (MQ=255)           cacCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAct                                                                                           <  2:524621/139‑1 (MQ=255)            acCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAActc                                                                                          >  1:802415/1‑139 (MQ=255)                aGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCtt                                                                                      >  1:71894/1‑139 (MQ=255)                  acaacaTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCg                                                                                    <  2:71894/139‑1 (MQ=255)                          acaGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAct                                                                                           <  1:338384/124‑1 (MQ=255)                          acaGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAct                                                                                           >  2:338384/1‑124 (MQ=255)                                        ggATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGc                                                              >  2:926001/1‑139 (MQ=255)                                                   gAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAgg                                                   <  2:802415/139‑1 (MQ=255)                                                   gAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACCTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAgg                                                   <  1:632933/139‑1 (MQ=255)                                                         aaaaGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGccc                                             <  2:886473/139‑1 (MQ=255)                                                               ttAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCa                                       <  1:926001/139‑1 (MQ=255)                                                                                 cTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCtttt                     >  2:295132/1‑139 (MQ=255)                                                                                                    tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTggg                                         <  2:343818/100‑1 (MQ=255)                                                                                                    tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTggg                                         >  1:343818/1‑100 (MQ=255)                                                                                                    tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGc                             <  1:679131/112‑1 (MQ=255)                                                                                                    tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGc                             >  2:679131/1‑112 (MQ=255)                                                                                                    tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTtgtg  >  2:295798/1‑139 (MQ=255)                                                                                                    tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTtgtg  >  2:543904/1‑139 (MQ=255)                                                                                                    tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTtgtg  >  1:416801/1‑139 (MQ=255)                                                                                                              tGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGaa                                                                                  >  1:708964/1‑49 (MQ=255)                                                                                                              tGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGaa                                                                                  <  2:708964/49‑1 (MQ=255)                                                                                                              tGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGccagccag            <  2:314800/119‑1 (MQ=255)                                                                                                              tGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGccagccag            >  1:314800/1‑119 (MQ=255)                                                                                                                                 |                                                                                                              CTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTG  >  NC_000913/3683502‑3683739

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