Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,407,859	A→G	100%	intergenic (‑64/‑191)	ydaM ← / → ydaN	diguanylate cyclase, csgD regulator/putative Zn(II) transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,407,859	0	A	G	100.0%	281.4 / ‑6.7	84	intergenic (‑64/‑191)	ydaM/ydaN	diguanylate cyclase, csgD regulator/putative Zn(II) transporter
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (37/44);  total (38/46)

AAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAAT  >  NC_000913/1407764‑1407954                                                                                                |                                                                                                aaGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAg                                                                                             <  1:1234312/100‑1 (MQ=255) aaGTCCAGGATATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGTAAAAAAAAAAAAAAAAAGCCACGTCATATGGTCAg                                                                                             <  2:1691805/100‑1 (MQ=255)  aGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGa                                                                                            >  1:1040888/1‑100 (MQ=255)  aGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGa                                                                                            >  2:185028/1‑100 (MQ=255)   ntCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTTTTTTTTTTTTATCGCAAAAAAGCCACGTCATATGGTCAGAt                                                                                           >  2:1972553/2‑100 (MQ=255)   gTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGAt                                                                                           >  2:398759/1‑100 (MQ=255)      naGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAAt                                                                                        >  1:896638/2‑100 (MQ=255)      cAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATa                                                                                          <  2:2393330/98‑1 (MQ=255)       aGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAAtt                                                                                       >  1:78362/1‑100 (MQ=255)        gggTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCCATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTc                                                                                      >  1:528716/1‑100 (MQ=255)         ggTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCt                                                                                     <  1:2397443/100‑1 (MQ=255)         ggTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCt                                                                                     <  1:2397376/100‑1 (MQ=255)           tATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGa                                                                                   <  2:2141939/100‑1 (MQ=255)           tATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGa                                                                                   >  1:427314/1‑100 (MQ=255)           tATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGa                                                                                   <  1:406105/100‑1 (MQ=255)             ttGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATa                                                                                 <  1:2110583/100‑1 (MQ=255)               gAAGTTGTGCGTAATCATTGAGATGCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCCAAAAAGCCACGTCATATGGTCAGATAATTCTGATAAt                                                                               >  2:1130424/1‑100 (MQ=255)               gAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAAt                                                                               <  1:371111/100‑1 (MQ=255)                aaGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATg                                                                              >  2:2101602/1‑100 (MQ=255)                   ttGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATa                                                                           <  2:2077295/100‑1 (MQ=255)                   ttGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATa                                                                           <  2:1674439/100‑1 (MQ=255)                    tgtgCGTAATCATTGAGATCCCATTATTTTTAGCAATTTTTGTTGCTATACAATAAAAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATa                                                                           <  2:146761/99‑1 (MQ=255)                     gtgCGTAATCATTGAGATCCGGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGa                                                                         >  2:1561805/1‑100 (MQ=255)                      tgCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGAc                                                                        >  2:1064188/1‑100 (MQ=255)                      tgCGTAATCAGTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCTTATGGTCAGATAATTCTGATAATGATAGAc                                                                        <  2:1690120/100‑1 (MQ=255)                        cGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGc                                                                      <  1:1713557/100‑1 (MQ=255)                        cGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGc                                                                      <  1:1713614/100‑1 (MQ=255)                        cGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGc                                                                      >  1:1998460/1‑100 (MQ=255)                         atAATCATTAAGATCCCATTATTTTTATCAATATTTGTTACTATCCGATCAAAAAAAATCCACATAATATAAACAAATAATACTAATAATAATAAAca                                                                       >  2:992198/2‑97 (MQ=255)                           aaTCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTAt                                                                   <  2:325552/100‑1 (MQ=255)                           aaTCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTAt                                                                   >  1:751924/1‑100 (MQ=255)                             tCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTAttt                                                                 <  1:930308/100‑1 (MQ=255)                              cATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAAGGATATACGCTATTTa                                                                >  2:2143598/1‑100 (MQ=255)                               aTTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTa                                                                <  2:135018/99‑1 (MQ=255)                                 tGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAaca                                                             <  1:1242406/100‑1 (MQ=255)                                 tGAGATCCCGTTATTTTTATCAATTTTTGTTGCTAAAAAATAAAAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAaca                                                             <  2:1834393/100‑1 (MQ=255)                                  gagaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATCATTCTGATAATGATAGACGCTATTTAacac                                                            <  1:515542/100‑1 (MQ=255)                                  gagaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAacac                                                            >  2:2412408/1‑100 (MQ=255)                                  gagaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAacac                                                            >  1:1497195/1‑100 (MQ=255)                                    gaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACtt                                                          >  1:2311180/1‑100 (MQ=255)                                       cccGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTcac                                                       >  1:2531304/1‑100 (MQ=255)                                         cGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCacac                                                     <  2:1715056/100‑1 (MQ=255)                                         cGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCacac                                                     >  2:1631349/1‑100 (MQ=255)                                         cGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCacac                                                     <  2:2461429/100‑1 (MQ=255)                                              tttttATCAATTTTTGTTTCAAAAAAAAAAAAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGttt                                                <  2:2286483/100‑1 (MQ=255)                                              tttttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTTAGATAATTCTGATAATGCTAGACGCTATTTAACACTTCACACGGttt                                                >  2:981548/1‑100 (MQ=255)                                                tttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGt                                              <  1:1961926/100‑1 (MQ=255)                                                 ttATCAATTTTTGTTTCTATACGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGta                                             <  2:249364/100‑1 (MQ=255)                                                  tATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGtat                                            >  1:291858/1‑100 (MQ=255)                                                  tATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGtat                                            <  1:713827/100‑1 (MQ=255)                                                     gtggcttttttgCGATCGGATAGCAACAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACg                                         >  1:710274/10‑100 (MQ=25)                                                       aTTTTCGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAAATCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGa                                       <  1:2375949/100‑1 (MQ=255)                                                           ttgttgCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAg                                   <  2:1267020/100‑1 (MQ=255)                                                             gttgCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCa                                 >  2:644730/1‑100 (MQ=255)                                                             gttgCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCa                                 <  2:1894956/100‑1 (MQ=255)                                                              ttgCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTTATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCa                                 <  1:2410180/99‑1 (MQ=255)                                                               tgCTATCGGATCGCAAAAAAGCGACGTCATATTCTAAGATAATTCTGATAATGATAGACACTAATTAACACTTCACACGGTTTGTATACGGAAAAGCAtt                               <  2:777491/100‑1 (MQ=255)                                                               tgCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtt                               >  1:2288067/1‑100 (MQ=255)                                                                gCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAttt                              >  2:81920/1‑100 (MQ=255)                                                                gCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAttt                              >  1:114779/1‑100 (MQ=255)                                                                gCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAttt                              >  2:1079964/1‑100 (MQ=255)                                                                  tATCCGATCGCAAAAAAGCAACGTCATATGGTTAGATAATTCTGATAATGATAGACGCTATTTTATACTTCACACGGTTTGTATACGGAAAAGCATTTTg                            <  2:1102865/100‑1 (MQ=255)                                                                   aTCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGc                           >  1:1968444/1‑100 (MQ=255)                                                                   aTCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGc                           >  1:2207900/1‑100 (MQ=255)                                                                      cGATCGCAAAAAAGCCACGTCATATGTTCATATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCtt                         <  2:1129008/99‑1 (MQ=255)                                                                       gATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCtttt                       <  1:1448071/100‑1 (MQ=255)                                                                       gATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCtttt                       >  1:1926501/1‑100 (MQ=255)                                                                        aTCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCttttt                      <  2:427310/100‑1 (MQ=255)                                                                         tCGCAAAAAAGGCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTGCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTg                     <  1:284957/100‑1 (MQ=255)                                                                          cgaaaaaaaGCCGCGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTGAACACTTCGCACGGTTTGTATACGGAAAAGCATTTTGCTTTTTg                     <  2:2317248/96‑1 (MQ=255)                                                                          cGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGt                    >  1:1535502/1‑100 (MQ=255)                                                                            cAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTAt                  <  2:701823/100‑1 (MQ=255)                                                                              aaaaaGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTc                >  1:1603134/1‑100 (MQ=255)                                                                                 aaGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAAt             <  2:548394/100‑1 (MQ=255)                                                                                   gCCACGTCATATGGTCAGATAATTATGATAATGATAGATGCTATTTAACACTTCACCCTGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAAttt           <  2:1395909/100‑1 (MQ=255)                                                                                     cACGTCATATGGTCAGATAATTCTGATAATGATAGGCGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAg         >  2:1638583/1‑100 (MQ=255)                                                                                     cACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAg         <  2:2066622/100‑1 (MQ=255)                                                                                      aCGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGa        <  2:816169/100‑1 (MQ=255)                                                                                       cGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGAc       <  2:291854/100‑1 (MQ=255)                                                                                         tCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAg     <  1:451565/100‑1 (MQ=255)                                                                                         tCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAg     >  1:54801/1‑100 (MQ=255)                                                                                           atatGGTTAAATAATTCTGATAATGATAGATTCTATTTAATACTTTATATGGTTTGTATACGGAAAAGTATTTTGCTTTTTGTATTCAATTTAGat       <  2:799445/96‑2 (MQ=255)                                                                                            tatGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAAt  >  1:1948086/1‑100 (MQ=255)                                                                                            tatGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAAt  >  2:2540581/1‑100 (MQ=255)                                                                                                |                                                                                                AAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAAT  >  NC_000913/1407764‑1407954

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