Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	404,305	A→G	E201E (GAA→GAG)	dgcC →	diguanylate cyclase DgcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	404,305	0	A	G	100.0%	147.7 / NA	43	E201E (GAA→GAG)	dgcC	diguanylate cyclase DgcC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/17);  total (26/17)

GTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGA  >  NC_000913/404168‑404438                                                                                                                                          |                                                                                                                                      gTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTAGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCa                                                                                                                                     >  1:121999/1‑140 (MQ=255)            gCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGt                                                                                                                          <  2:974099/140‑1 (MQ=255)              tGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTg                                                                                                                        <  2:866633/140‑1 (MQ=255)                aCGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCa                                                                                                                      <  1:130180/140‑1 (MQ=255)                      aTTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCAt                                                                                                                <  1:982717/140‑1 (MQ=255)                             tgtCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAGCTGGCGGAGCATAAACGCAGGATGCAGGTCATGAGTAcc                                                                                                         >  2:263024/1‑140 (MQ=255)                             tgtCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTAcc                                                                                                         >  1:521661/1‑140 (MQ=255)                                 gTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCcgcg                                                                                                     <  2:974924/140‑1 (MQ=255)                                   tCAATAGTGCGCCGCTGGAATGTTGACTCACCCTTCCTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGAc                                                                                                   <  1:411803/140‑1 (MQ=255)                                        aGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACGGGCGGAGCATAAGCTCGCGTTGCAGATCATCGGTACCCGCGACGGCAt                                                                                              >  2:492484/1‑140 (MQ=255)                                         gTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATg                                                                                             >  1:611101/1‑140 (MQ=255)                                          tGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGa                                                                                            >  1:13764/1‑140 (MQ=255)                                            cgccgcTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGTAGGTCAGGAGTACACGCGACGGCATGACg                                                                                          >  1:300268/1‑140 (MQ=255)                                            cgccgcTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACg                                                                                          >  2:1001708/1‑140 (MQ=255)                                            cgccgcTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACg                                                                                          >  1:796211/1‑140 (MQ=255)                                             gccgcTGGAATGGTGGCTCTCCCTTCCCATTCTTGCCATCTATCCTCTTCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACgg                                                                                         <  2:300268/140‑1 (MQ=255)                                                           ggCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGa                                                                           <  1:925088/140‑1 (MQ=255)                                                                tcCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCa                                                                      <  1:239583/140‑1 (MQ=255)                                                                  ccTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCAtt                                                                    >  1:876831/1‑140 (MQ=255)                                                                   cTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTg                                                                   <  2:120257/140‑1 (MQ=255)                                                                      cccATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGa                                                                <  2:796211/140‑1 (MQ=255)                                                                      cccATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGGCAGGAGTACCCACGAAGACGTGACGGGCGTGTAGAACCGACGTCATTGGGa                                                                >  2:451909/1‑140 (MQ=255)                                                                        cATTATTGTCATTTAGCCTCTGCTGTTTGACTGGGTCCGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGaaa                                                              <  1:451909/140‑1 (MQ=255)                                                                          ttattGTCATTTATCGTCTGCTATTAGGCCGGGTCAGTTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACt                                                            <  1:474773/140‑1 (MQ=255)                                                                          ttattGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGTGCGGGTGTAACCGTCGTCATTGGGAAACt                                                            >  1:372848/1‑140 (MQ=255)                                                                              tGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGt                                                        >  2:974376/1‑140 (MQ=255)                                                                              tGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCAGCGACGGCATGACGGGCGTGTATAATCGACTTCATTGGGAAACTATGt                                                        >  2:333071/1‑140 (MQ=255)                                                                                tCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTa                                                      >  2:593318/1‑140 (MQ=255)                                                                                  aTTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACg                                                    >  1:519282/1‑140 (MQ=255)                                                                                       tCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATg                                               >  2:854190/1‑140 (MQ=255)                                                                                       tCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATg                                               >  1:253637/1‑140 (MQ=255)                                                                                       tCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATg                                               <  2:521661/140‑1 (MQ=255)                                                                                              ctgTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTGATTGGGAAACTATGGTACGCAGGGAATTTGa                                        >  2:461918/1‑140 (MQ=255)                                                                                                      cTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCAATGGGAAAATATGTTACGCAATGAATTTGATAACTGTc                                >  1:263779/1‑140 (MQ=255)                                                                                                           tCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTcggcgg                           >  1:132607/1‑140 (MQ=255)                                                                                                           tCAGCTACCAGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTcggcgg                           >  2:1019799/1‑140 (MQ=255)                                                                                                                    aGACGGCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATcgc                  >  2:201753/1‑140 (MQ=255)                                                                                                                        ggCAACCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATg              >  1:963946/1‑140 (MQ=255)                                                                                                                          cAACCAAACTGGCGGTGCATAAGCGTAGGTTGCAGGTGATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCa            <  1:492484/140‑1 (MQ=255)                                                                                                                            aCCAAACTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAAc          >  1:431265/1‑140 (MQ=255)                                                                                                                                aaCTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTa      <  2:540658/140‑1 (MQ=255)                                                                                                                                aaCTGGCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTa      <  2:132607/140‑1 (MQ=255)                                                                                                                                    ggCGGAGCATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGa  >  1:790673/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                      GTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGA  >  NC_000913/404168‑404438

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