Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,093,358	(G)6→7	coding (483/1359 nt)	dgcT →	putative diguanylate cyclase DgcT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,093,352	1	.	G	100.0%	128.3 / NA	36	F159L (TTT→TTG)	dgcT	putative diguanylate cyclase DgcT
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (15/21);  total (15/21)

CTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTT‑GGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGT  >  NC_000913/1093217‑1093485                                                                                                                                           |                                                                                                                                      cTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTT‑gggg                                                                                                                                   <  2:402258/140‑1 (MQ=255) cTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTT‑gggg                                                                                                                                   <  2:367534/140‑1 (MQ=255) cTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTT‑gggg                                                                                                                                   <  1:672841/140‑1 (MQ=255)     tattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTT‑gggggggc                                                                                                                               >  1:299091/1‑138 (MQ=255)     tattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTT‑gggggggc                                                                                                                               >  1:85548/1‑138 (MQ=255)      attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGAGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTT‑gggggggcc                                                                                                                              <  1:589088/140‑4 (MQ=255)          gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATa                                                                                                                           >  2:688347/1‑140 (MQ=255)            ttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATa                                                                                                                           >  2:273697/1‑140 (MQ=255)            ttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGTCGATAGt                                                                                                                         >  2:716757/1‑140 (MQ=255)            ttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGt                                                                                                                         <  2:740562/140‑1 (MQ=255)              cGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTT‑gggg                                                                                                                                   <  1:156020/127‑1 (MQ=255)              cGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTT‑gggg                                                                                                                                   >  2:156020/1‑127 (MQ=255)                tCAGATGAGTTTGTGTATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCt                                                                                                                     <  2:299091/140‑1 (MQ=255)                tCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCt                                                                                                                     <  2:85548/140‑1 (MQ=255)                tCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCt                                                                                                                     <  2:330795/140‑1 (MQ=255)                tCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCt                                                                                                                     <  1:688347/140‑1 (MQ=255)                    aTGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTcaca                                                                                                                 <  1:768159/140‑1 (MQ=255)                     tGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTcac                                                                                                                  <  1:273697/140‑1 (MQ=255)                           tgtgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTc                                                                                                          <  2:622531/140‑1 (MQ=255)                           tgtgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTc                                                                                                          <  2:280880/140‑1 (MQ=255)                             tgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTGCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGTCGATAGTTGCTCACATACTTTCTa                                                                                                        <  1:716757/140‑1 (MQ=255)                                            cTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTAt                                                                                                >  2:664045/1‑133 (MQ=255)                                            cTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTAt                                                                                                <  1:664045/133‑1 (MQ=255)                                            cTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCt                                                                                         >  2:386022/1‑140 (MQ=255)                                            cTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCt                                                                                         >  1:231443/1‑140 (MQ=255)                                            cTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCt                                                                                         <  2:231443/140‑1 (MQ=255)                                                cATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACaa                                                                                     >  1:618745/1‑140 (MQ=255)                                                     gTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCt                                                                                         <  1:1189/131‑1 (MQ=255)                                                     gTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCt                                                                                         >  2:1189/1‑131 (MQ=255)                                                              gTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCg                                                                       >  2:280840/1‑140 (MQ=255)                                                               ttAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGa                                                                      >  1:826210/1‑140 (MQ=255)                                                               ttAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGa                                                                      >  2:385101/1‑140 (MQ=255)                                                                    gAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTaa                                                                 >  1:609426/1‑140 (MQ=255)                                                                       aaCACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCa                                                              >  1:696805/1‑140 (MQ=255)                                                                             aaaCAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGaaa                                                        <  2:826210/140‑1 (MQ=255)                                                                                    gcAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCa                                                 >  2:81985/1‑140 (MQ=255)                                                                                         ttatttTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAggtgg                                            <  1:310972/140‑1 (MQ=255)                                                                                          tatttTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAggtggt                                           <  1:81985/140‑1 (MQ=255)                                                                                                ccAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGaa                                     <  1:280840/140‑1 (MQ=255)                                                                                                ccAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGaa                                     <  1:261811/140‑1 (MQ=255)                                                                                                  aaaaaaaTGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAg                                   >  2:123768/1‑140 (MQ=255)                                                                                                    aaaaaTGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGcc                                 <  1:123768/140‑1 (MQ=255)                                                                                                     aaaaTGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCt                                <  2:696805/140‑1 (MQ=255)                                                                                                     aaaaTGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCt                                <  2:635114/140‑1 (MQ=255)                                                                                                                             tgttttttgtttttgGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGAtttt          >  2:375425/1‑140 (MQ=255)                                                                                                                                    gtttttg‑gGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGt  <  2:658836/138‑1 (MQ=255)                                                                                                                                    gtttttg‑gGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGt  <  1:385101/138‑1 (MQ=255)                                                                                                                                           |                                                                                                                                      CTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT‑‑TGTGCATAAGTTTGTTTTTTGTTTTT‑GGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGT  >  NC_000913/1093217‑1093485

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: ‑