Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,567,285	T→C	intergenic (‑145/+219)	dosC ← / ← yddW	diguanylate cyclase, cold‑ and stationary phase‑induced oxygen‑dependent biofilm regulator/liprotein, glycosyl hydrolase homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,567,285	0	T	C	100.0%	188.1 / NA	53	intergenic (‑145/+219)	dosC/yddW	diguanylate cyclase, cold‑ and stationary phase‑induced oxygen‑dependent biofilm regulator/liprotein, glycosyl hydrolase homolog
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/26);  total (27/26)

CTGGGAAATCAGACGGAATCAAAATGAAACGCAACGTGCGAGATCGACTAACTGCACCATATTCTCCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCTACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGT  >  NC_000913/1567160‑1567419                                                                                                                              |                                                                                                                                       cTGGGAAATCAGACGGAATCAAAATGAAACGCAACGTGCGAGATCGACTAACTGCACCATATTCTCCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCaga                                                                                                                           <  1:34808/139‑1 (MQ=255)  tGGGAAATCAGACGGAATCAAAATGAAACGCAACGTGCGAGATCGACTAACTGCACCATATTCTCCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCagaa                                                                                                                          <  2:1033001/139‑1 (MQ=255)    ggAAATCAGACGGAATCAAAATGAAACGCAACGTGCGAGATCGACTAACTGCACCATATTCTCCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCagaaga                                                                                                                        >  2:189753/1‑139 (MQ=255)    ggAAATCAGACGGAATCAAAATGAAACGCAACGTGCGAGATCGACTAACTGCACCATATTCTCCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGAAAGTTTCAATCTACGCCATTGTTCTTCagaaga                                                                                                                        >  2:1135378/1‑139 (MQ=255)                           aaaCGCAACGTGCGAGATCGACTAACTGCACCATATTCTCCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTAtttt                                                                                                 >  2:32375/1‑139 (MQ=255)                               gCAACGTGCGAGATCGACTAACTGCACCATATTCTCCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCgg                                                                                             <  1:189753/139‑1 (MQ=255)                                               aCTAACTGCACCATATTCTCCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTg                                                                                                                                    >  2:443803/1‑84 (MQ=255)                                               aCTAACTGCACCATATTCTCCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTg                                                                                                                                    <  1:443803/84‑1 (MQ=255)                                                cTAACTGCACCATATTCTCCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTa                                                                            >  1:1114718/1‑139 (MQ=255)                                                       cACCATATTCTCCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTg                                                                                                                                    <  2:64308/76‑1 (MQ=255)                                                       cACCATATTCTCCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTg                                                                                                                                    >  1:64308/1‑76 (MQ=255)                                                         ccATATTCTCCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTttc                                                                   <  1:1129055/139‑1 (MQ=255)                                                               tctcCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTttcatattc                                                             >  2:1482658/1‑139 (MQ=255)                                                                 tcCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCag                                                                                                                            <  2:1130898/74‑1 (MQ=255)                                                                 tcCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTttcatattcat                                                           >  2:346561/1‑139 (MQ=255)                                                                 tcCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTAAATCTACGCCACTGTTCTTCag                                                                                                                            >  1:1130898/1‑74 (MQ=255)                                                                   cTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTg                                                                                           <  2:194040/105‑1 (MQ=255)                                                                   cTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTg                                                                                           >  1:194040/1‑105 (MQ=255)                                                                      aaaTATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCTTACAtt                                                      >  2:598489/1‑139 (MQ=255)                                                                           tGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCAt                                                                                                             >  2:1477613/1‑79 (MQ=255)                                                                           tGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCAt                                                                                                             <  1:1477613/79‑1 (MQ=255)                                                                           tGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTttat                                                 <  1:32375/139‑1 (MQ=255)                                                                            gAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTttatt                                                <  1:1482658/139‑1 (MQ=255)                                                                                atatACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAggg                                            >  2:1300157/1‑139 (MQ=255)                                                                                        aaaaGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACa                                                                                        <  2:430453/87‑1 (MQ=255)                                                                                        aaaaGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACa                                                                                        >  1:430453/1‑87 (MQ=255)                                                                                        aaaaGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATgg                                    <  1:1300157/139‑1 (MQ=255)                                                                                        aaaaGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATgg                                    <  1:598489/139‑1 (MQ=255)                                                                                          aaGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGttttt                                                                    <  1:71159/105‑1 (MQ=255)                                                                                          aaGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGttttt                                                                    >  2:71159/1‑105 (MQ=255)                                                                                          aaGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACAt                                                       <  2:393931/118‑1 (MQ=255)                                                                                          aaGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACAt                                                       >  1:393931/1‑118 (MQ=255)                                                                                          aaGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTttatt                                                <  2:346865/125‑1 (MQ=255)                                                                                          aaGAAATAAGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTttatt                                                >  1:346865/1‑125 (MQ=255)                                                                                                  aGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCag                                                                                                                            >  1:57792/1‑41 (MQ=255)                                                                                                  aGCGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCag                                                                                                                            <  2:57792/41‑1 (MQ=255)                                                                                                    cGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGa                                                                                          >  1:822756/1‑73 (MQ=255)                                                                                                    cGATTTAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGa                                                                                          <  2:822756/73‑1 (MQ=255)                                                                                                     gATTTAGGACAGTTTCAATCTACGCCACTGttctt                                                                                                                               <  2:121695/35‑1 (MQ=255)                                                                                                     gATTTAGGACAGTTTCAATCTACGCCACTGttctt                                                                                                                               >  1:121695/1‑35 (MQ=255)                                                                                                       tttAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGta                                                                                                                     >  2:1256305/1‑43 (MQ=255)                                                                                                       tttAGGACAGTTTCAATCTACGCCACTGTTCTTCAGAAGAGta                                                                                                                     <  1:1256305/43‑1 (MQ=255)                                                                                                                 tttCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTg                                                                                           >  1:910093/1‑59 (MQ=255)                                                                                                                 tttCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTg                                                                                           <  2:910093/59‑1 (MQ=255)                                                                                                                 tttCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGAc                                                                                         <  1:925842/61‑1 (MQ=255)                                                                                                                 tttCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGAc                                                                                         >  2:925842/1‑61 (MQ=255)                                                                                                                 tttCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTaa                            >  2:605815/1‑122 (MQ=255)                                                                                                                 tttCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTaa                            <  1:605815/122‑1 (MQ=255)                                                                                                                 tttCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAAttt           >  1:862089/1‑139 (MQ=255)                                                                                                                 tttCAATCTACGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAAttt           >  1:1129028/1‑139 (MQ=255)                                                                                                                          aCGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTtaagt  <  2:1129028/139‑1 (MQ=255)                                                                                                                           cGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTttcata                                                                <  2:1377651/76‑1 (MQ=255)                                                                                                                           cGCCACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTttcata                                                                >  1:1377651/1‑76 (MQ=255)                                                                                                                              |                                                                                                                                       CTGGGAAATCAGACGGAATCAAAATGAAACGCAACGTGCGAGATCGACTAACTGCACCATATTCTCCTGAAATATGAAGATATACTGAAAAGAAATAAGCGATTTAGGACAGTTTCAATCTACGCTACTGTTCTTCAGAAGAGTATAGCCCATCGTAATTATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGT  >  NC_000913/1567160‑1567419

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