Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,878,402	C→A	L230M (CTG→ATG)	dgcP →	diguanylate cyclase DgcP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,878,402	0	C	A	100.0%	178.2 / NA	51	L230M (CTG→ATG)	dgcP	diguanylate cyclase DgcP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (30/21);  total (30/21)

CGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAAC  >  NC_000913/1878267‑1878537                                                                                                                                        |                                                                                                                                        cggcggGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGtt                                                                                                                                     >  1:934342/1‑140 (MQ=255)  ggcgggcgATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGttt                                                                                                                                    >  1:238501/1‑140 (MQ=255)     gggcgATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTtctc                                                                                                                                 <  2:747223/140‑1 (MQ=255)       gcgATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCAt                                                                                                                               <  2:647794/140‑1 (MQ=255)       gcgATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCAt                                                                                                                               <  2:1319619/140‑1 (MQ=255)          aTTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCa                                                                                                                            >  2:1006930/1‑140 (MQ=255)            ttttGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAgg                                                                                                                          <  2:688110/140‑1 (MQ=255)                    aTCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATGACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAg                                                                                                                  <  2:165551/140‑1 (MQ=255)                    aTCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAg                                                                                                                  <  1:1019514/140‑1 (MQ=255)                    aTCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAg                                                                                                                  <  2:1263754/140‑1 (MQ=255)                    aTCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAg                                                                                                                  <  2:744314/140‑1 (MQ=255)                    aTCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAg                                                                                                                  <  1:1189532/140‑1 (MQ=255)                    aTCAGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAg                                                                                                                  <  1:421818/140‑1 (MQ=255)                              acTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATa                                                                                                        <  2:598539/140‑1 (MQ=255)                                    ttCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTcc                                                                                                  >  1:276759/1‑140 (MQ=255)                                            cccGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAAt                                                                                          <  1:312551/140‑1 (MQ=255)                                            cccGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAAt                                                                                          >  2:1140270/1‑140 (MQ=255)                                                gCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCg                                                                                      <  2:276759/140‑1 (MQ=255)                                                           aTAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGgt                                                                           >  2:1040543/1‑140 (MQ=255)                                                             aaGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGgtcg                                                                         >  2:1102245/1‑140 (MQ=255)                                                             aaGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGgtcg                                                                         >  2:1095732/1‑140 (MQ=255)                                                              aGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGgtcgt                                                                        >  2:656732/1‑140 (MQ=255)                                                               gATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGgtcgtc                                                                       >  2:137536/1‑140 (MQ=255)                                                                  aaTGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTcg                                                                    >  2:1009233/1‑140 (MQ=255)                                                                   aTGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTcgg                                                                   >  2:455860/1‑140 (MQ=255)                                                                     gATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTcggcg                                                                 <  2:1132863/140‑1 (MQ=255)                                                                        cgcgTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGtg                                                              >  2:811211/1‑140 (MQ=255)                                                                        cgcgTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGtg                                                              >  1:213776/1‑140 (MQ=255)                                                                        cgcgTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGtg                                                              >  2:406990/1‑140 (MQ=255)                                                                          cgTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGtggt                                                            >  2:1191056/1‑138 (MQ=255)                                                                          cgTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGtgat                                                            <  1:838578/140‑1 (MQ=255)                                                                          cgTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGtgat                                                            <  1:1261773/140‑1 (MQ=255)                                                                          cgTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGtgat                                                            <  2:213776/140‑1 (MQ=255)                                                                                    tCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGgtttttt                                                  >  2:1177598/1‑138 (MQ=255)                                                                                    tCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAg                                                  <  1:811211/140‑1 (MQ=255)                                                                                        gATAACTTCAAATTAGTCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGGTTTTAgttgt                                              >  2:887256/1‑140 (MQ=255)                                                                                         aTAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAtttgtt                                             >  1:389587/1‑140 (MQ=255)                                                                                             cTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTTTTTGTTTTAc                                         >  2:288325/1‑140 (MQ=255)                                                                                               tCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTa                                       <  1:583649/140‑1 (MQ=255)                                                                                               tCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTa                                       <  1:1102245/140‑1 (MQ=255)                                                                                                cAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTaa                                      >  1:211647/1‑140 (MQ=255)                                                                                                  aaTTAATCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTTGTTGTTTCACTaaaa                                    >  1:423130/1‑139 (MQ=255)                                                                                                       aTCAATGATCGTTTTGGTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGAAGAGGTTTTTGTTGTTTTAATAAACAACga                               >  2:481188/1‑140 (MQ=255)                                                                                                                      ggTCATAATAGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTtc                >  1:339257/1‑140 (MQ=255)                                                                                                                           taataGTGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGc           <  2:211647/140‑1 (MQ=255)                                                                                                                                 tGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTTTTTTACTAAACAACGAGAATGGGGATATTTCGTCGCTGCgag     >  2:196540/1‑140 (MQ=255)                                                                                                                                 tGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCgag     >  1:387680/1‑140 (MQ=255)                                                                                                                                 tGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCgag     >  1:1013054/1‑140 (MQ=255)                                                                                                                                 tGGCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGGTTTTAGTTGTTTTACTAAACAAAGAGAATGCGGATATTTCGTCGCTGCgag     >  1:17358/1‑140 (MQ=255)                                                                                                                                   gCGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGaa   >  1:604452/1‑140 (MQ=255)                                                                                                                                    cGATATGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATTTTTTGTCGCTGCGAGAAc  >  1:338016/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                        CGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAAC  >  NC_000913/1878267‑1878537

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