Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,872,031	G→A	intergenic (+171/‑10)	cdgI → / → dgcJ	putative c‑di‑GMP binding protein CdgI/putative diguanylate cyclase DgcJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,872,031	0	G	A	100.0%	211.5 / NA	61	intergenic (+171/‑10)	cdgI/dgcJ	putative c‑di‑GMP binding protein CdgI/putative diguanylate cyclase DgcJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/27);  total (34/27)

TTTATTACTATTACCCCTAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCG  >  NC_000913/1871909‑1872167                                                                                                                           |                                                                                                                                         tttATTACTATTACCCCTAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATATGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAAtt                                                                                                                         >  1:1138881/1‑140 (MQ=255)    attactattacCCCTAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCa                                                                                                                      <  2:709795/140‑1 (MQ=255)         tattacCCCTAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCa                                                                                                                      <  2:446477/135‑1 (MQ=255)         tattacCCCTAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCa                                                                                                                      >  1:446477/1‑135 (MQ=255)         tattacCCCTAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATa                                                                                                                 >  2:721076/1‑140 (MQ=255)           ttacCCCTAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCa                                                                                                                      >  2:884303/1‑133 (MQ=255)           ttacCCCTAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCa                                                                                                                      <  1:884303/133‑1 (MQ=255)           ttacCCCTAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGa                                                                                                               >  1:943300/1‑140 (MQ=255)           ttacCCCTAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGa                                                                                                               >  1:520838/1‑140 (MQ=255)             acCCCTAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAAt                                                                                                             >  1:275101/1‑140 (MQ=255)             acCCCTAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAAt                                                                                                             <  2:520838/140‑1 (MQ=255)                 cTAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTc                                                                                                         <  1:1106906/140‑1 (MQ=255)                  tAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAAtt                                                                                                                         <  1:121815/123‑1 (MQ=255)                  tAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAAtt                                                                                                                         >  2:121815/1‑123 (MQ=255)                     aTTTCCCCTCCCCAGCGGCGCAGAGGTGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGc                                                                                                     >  2:217311/1‑140 (MQ=255)                         cccctccccAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCAttt                                                                                                 >  2:933494/1‑140 (MQ=255)                         cccctccccAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCAttt                                                                                                 >  2:395839/1‑140 (MQ=255)                              ccccAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATcg                                                                                            >  2:299540/1‑140 (MQ=255)                               cccAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATcgc                                                                                           >  1:550035/1‑140 (MQ=255)                                ccAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATcgcc                                                                                          <  2:550035/140‑1 (MQ=255)                                ccAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATcgcc                                                                                          <  1:180137/140‑1 (MQ=255)                                  aGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATcgccgc                                                                                        <  2:943300/140‑1 (MQ=255)                                       cgcAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAgtg                                                                                   >  1:172796/1‑140 (MQ=255)                                                gAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGc                                                                          <  1:395839/140‑1 (MQ=255)                                                 aGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCt                                                                         >  1:851452/1‑140 (MQ=255)                                                      aaTTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACAt                                                                    >  1:667923/1‑140 (MQ=255)                                                        ttAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATct                                                                  >  1:1085963/1‑140 (MQ=255)                                                          aGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAGTGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATcttc                                                                >  2:1142211/1‑140 (MQ=255)                                                           gcgcCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATcttcc                                                               >  1:593443/1‑140 (MQ=255)                                                            cgcCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATcttcct                                                              >  2:8961/1‑140 (MQ=255)                                                             gcCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATcttcctt                                                             >  1:866895/1‑140 (MQ=255)                                                                ccTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTtccttcct                                                          <  1:721076/140‑1 (MQ=255)                                                                ccTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTtccttcct                                                          >  1:858459/1‑140 (MQ=255)                                                                ccTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTtccttcct                                                          <  1:247050/140‑1 (MQ=255)                                                                      ccAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTGACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTAtttt                                                    >  2:1027926/1‑140 (MQ=255)                                                                      ccAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTAtttt                                                    >  2:591139/1‑140 (MQ=255)                                                                         ggCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTtccttcct                                                          >  1:621719/1‑131 (MQ=255)                                                                         ggCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTtccttcct                                                          <  2:621719/131‑1 (MQ=255)                                                                         ggCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGa                                                 >  2:89648/1‑140 (MQ=255)                                                                              cAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTg                                            >  2:1092938/1‑140 (MQ=255)                                                                              cAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTg                                            >  2:100092/1‑140 (MQ=255)                                                                              cAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTg                                            >  1:495248/1‑140 (MQ=255)                                                                                   gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGcc                                       <  2:172796/140‑1 (MQ=255)                                                                                   gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGcc                                       <  1:100092/140‑1 (MQ=255)                                                                                        ttATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGa                                  <  2:275101/140‑1 (MQ=255)                                                                                         tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGAc                                 <  2:1085963/140‑1 (MQ=255)                                                                                                   tctcGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCTCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGa                       >  2:1010602/1‑140 (MQ=255)                                                                                                   tctcGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACagag                             >  2:305699/1‑134 (MQ=255)                                                                                                   tctcGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACagag                             <  1:305699/134‑1 (MQ=255)                                                                                                    ctcGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAg                                     >  2:512668/1‑125 (MQ=255)                                                                                                    ctcGTTTCACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAg                                     <  1:512668/125‑1 (MQ=255)                                                                                                           cACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTAt               <  1:1027926/140‑1 (MQ=255)                                                                                                           cACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTAt               <  2:161503/140‑1 (MQ=255)                                                                                                           cACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTAt               <  2:866895/140‑1 (MQ=255)                                                                                                           cACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTAt               <  2:667923/140‑1 (MQ=255)                                                                                                           cACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTAt               <  2:593443/140‑1 (MQ=255)                                                                                                           cACTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTAt               <  2:495248/140‑1 (MQ=255)                                                                                                             cTAACCGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCt             <  2:851452/140‑1 (MQ=255)                                                                                                                  cGAAGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCt        <  2:858459/140‑1 (MQ=255)                                                                                                                     aGGAGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCtata     >  1:69222/1‑140 (MQ=255)                                                                                                                        aGTACCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCg  >  1:165032/1‑140 (MQ=255)                                                                                                                           |                                                                                                                                         TTTATTACTATTACCCCTAAATTTCCCCTCCCCAGCGGCGCAGAGATGAGTATAATTAGCGCCCCTGTGCCAGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCG  >  NC_000913/1871909‑1872167

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