Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,143,684	(G)6→7	coding (419/3318 nt)	dgcE →	putative diguanylate cyclase DgcE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,143,678	1	.	G	100.0%	86.4 / NA	26	L138W (TTG→TGG)	dgcE	putative diguanylate cyclase DgcE
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (13/13);  total (13/13)

GGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑GGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACT  >  NC_000913/2143543‑2143806                                                                                                                                         |                                                                                                                                 ggACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑gggg                                                                                                                              >  2:65147/1‑140 (MQ=255)   acgacCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑gggggg                                                                                                                            <  1:771209/140‑1 (MQ=255)      acCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑gggggggtg                                                                                                                         <  1:125076/140‑4 (MQ=255)        cATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGCCGGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑gggggggggtt                                                                                                                       >  2:299461/1‑135 (MQ=255)                        gCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGa                                                                                                        >  2:74870/1‑140 (MQ=255)                          cGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCctgctg                                                                                                         >  1:143739/1‑137 (MQ=255)                          cGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCctgctg                                                                                                         <  2:143739/137‑1 (MQ=255)                                 ggggCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGgag                                                                                               <  2:266499/140‑1 (MQ=255)                                        tgctgcTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGc                                                                                                     >  2:310083/1‑127 (MQ=255)                                        tgctgcTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGc                                                                                                     <  1:310083/127‑1 (MQ=255)                                        tgctgcTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATcc                                                                                        >  2:313337/1‑140 (MQ=255)                                             ctACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCa                                                                                   >  1:860970/1‑140 (MQ=255)                                                       tgctgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCAttttt                                                                         <  1:74870/140‑1 (MQ=255)                                                           gcCGTGGTATAACCCCTTGCAAAATCTGGCTGTTTGGCTGCGTCTGGCACTCTGCAGCGCCATTGTTCCGCCTCTCTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGata                                                                     <  2:806249/140‑1 (MQ=255)                                                           gcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGata                                                                     <  1:65147/140‑1 (MQ=255)                                                             cGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATg                                                                   >  2:548552/1‑140 (MQ=255)                                                                            tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATc                                                    >  2:630140/1‑140 (MQ=255)                                                                            tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATc                                                    >  2:661997/1‑140 (MQ=255)                                                                            tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATc                                                    >  2:289273/1‑140 (MQ=255)                                                                            tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATc                                                    >  1:463288/1‑140 (MQ=255)                                                                                  aTCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCgg                                              <  2:463288/140‑1 (MQ=255)                                                                                             tGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATACGGGTACTGTCAGAATCCATCGGCGCTCTGGCAc                                   <  2:860970/140‑1 (MQ=255)                                                                                                 tGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGt                               >  2:640178/1‑140 (MQ=255)                                                                                                     tCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCg                           >  2:827327/1‑140 (MQ=255)                                                                                                         gCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTgg                       <  1:313337/140‑1 (MQ=255)                                                                                                                  aGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTAt              <  1:630140/140‑1 (MQ=255)                                                                                                                   gcgcCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTAtt             >  2:89518/1‑140 (MQ=255)                                                                                                                        aTTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAAc        <  1:548552/140‑1 (MQ=255)                                                                                                                         ttGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAAcc       <  1:289273/140‑1 (MQ=255)                                                                                                                              ccGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACt  <  1:299461/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                 GGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑GGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACT  >  NC_000913/2143543‑2143806

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