Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,027,332	T→A	E229V (GAA→GTA)	dgcQ ←	diguanylate cyclase DgcQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,027,332	0	T	A	100.0%	28.7 / NA	12	E229V (GAA→GTA)	dgcQ	diguanylate cyclase DgcQ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (6/6);  total (6/6)

CAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATC  >  NC_000913/2027111‑2027476                                                                                                                                                                                                                              |                                                                                                                                                 cAAAAATATTCCCTGTGTGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGAGACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCGGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACcgcn                                                                                                                     <  1:63290/251‑2 (MQ=255)                 tGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATATAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGTGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTtct                                                                                                      <  1:23405/250‑1 (MQ=255)                                ggAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGa                                                                                                            >  1:196137/1‑229 (MQ=255)                                                                                   gAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCACCAGCACGCCATACCAGTAGTTCTTACTGTCTACCGGAACACTGACGGTCACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGAGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGGTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGcgcg                                   >  1:358145/1‑251 (MQ=255)                                                                                              cGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTAATACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAAcg                                                    <  1:356775/223‑1 (MQ=255)                                                                                                             aaaaTTGCTGCATAGTACGCACGGGAATTCTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCGGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGTCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGGAAACAAGGTCGGCTGCGTCGAAACg         <  1:290602/251‑1 (MQ=255)                                                                                                                    ctgcATAGTACGCACGGGAATACTCATCCCAAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTAACGCGCGGTGACGATTTTTTCGTTGTGCATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATc  >  1:147669/1‑251 (MQ=255)                                                                                                                                          cTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACATGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAAcc                                                                                                                                >  1:196182/1‑103 (MQ=255)                                                                                                                                            cATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGCCGGTAACCTGCGGATCAGCCTTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCa                                                                                       >  1:199996/1‑142 (MQ=255)                                                                                                                                                          cATACNNGTNNNTANNACTGTCTACCGGAACACTGACGGNAACCTGCGGTTCAGTATTGCTGGCGTGTANCNGNTGCGNNGTGNACCAGCGTAc                                                                                                                         >  2:256945/1‑94 (MQ=255)                                                                                                                                                          cATACCAGTAGGTAGGACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTAGTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTAc                                                                                                                         <  1:256945/94‑1 (MQ=255)                                                                                                                                                                        nnaCTGTCNNCCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTtc                                                                                                       <  2:16106/96‑1 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                 CAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATC  >  NC_000913/2027111‑2027476

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