Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,878,561	G→T	E283* (GAA→TAA)	dgcP →	diguanylate cyclase DgcP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,878,561	0	G	T	100.0%	162.9 / NA	48	E283* (GAA→TAA)	dgcP	diguanylate cyclase DgcP
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/32);  total (16/32)

GCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGAGAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTAT  >  NC_000913/1878427‑1878680                                                                                                                                       |                                                                                                                        gcCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAatat                                                                                                                    >  1:227570/1‑140 (MQ=255)  ccTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATc                                                                                                                   <  1:532853/140‑1 (MQ=255)     tAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTTACTAAAAAACGAGAATTCGGATATTTTGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACt                                                                                                                >  2:815910/1‑140 (MQ=255)      aaTACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACtt                                                                                                               <  2:728266/140‑1 (MQ=255)            cTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGa                                                                                                         <  2:1127325/140‑1 (MQ=255)             tCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTTACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAAAACGTGGATAATATCACTTAGGGGAt                                                                                                        >  2:423420/1‑140 (MQ=255)             tCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGAt                                                                                                        >  1:125804/1‑140 (MQ=255)                   aaaaTGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATCCGTGGATAATATCACTTAGGTGATGTTGAt                                                                                                  <  1:199316/140‑1 (MQ=255)                   aaaaTGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGAt                                                                                                  <  1:595879/140‑1 (MQ=255)                   aaaaTGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGAt                                                                                                  <  2:1282453/140‑1 (MQ=255)                   aaaaTGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGAt                                                                                                  <  2:1502967/140‑1 (MQ=255)                   aaaaTGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGAt                                                                                                  <  2:227570/140‑1 (MQ=255)                      aTGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTg                                                                                               <  1:449956/140‑1 (MQ=255)                                   tGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGAt                                                                                                  <  2:1044224/124‑1 (MQ=255)                                   tGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGAt                                                                                                  >  1:1044224/1‑124 (MQ=255)                                    ggTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGc                                                                                 <  1:677942/140‑1 (MQ=255)                                         tctcGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGcc                                                                                >  1:1034463/1‑136 (MQ=255)                                         tctcGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGcc                                                                                <  2:1034463/136‑1 (MQ=255)                                         tctcGGCGGTGATGAGTTTTTAGTTGTTTCACTAAAAAAAGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTTAGCAGCAAAAACGTGGATAATATCACTTAGGGGATGTTGATTTGTATTATCCCGGTGCCAGTc                                                                            >  1:856950/1‑140 (MQ=255)                                          ctcGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCt                                                                           <  1:485890/140‑1 (MQ=255)                                              gcggTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGc                                                                       <  2:125804/140‑1 (MQ=255)                                               cggTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCa                                                                      <  2:336504/140‑1 (MQ=255)                                                         tttttAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGt                                                            <  1:1130062/140‑1 (MQ=255)                                                             tAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGAt                                                        <  2:191803/140‑1 (MQ=255)                                                              aGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATc                                                       <  1:522994/140‑1 (MQ=255)                                                              aGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATc                                                       <  2:1278043/140‑1 (MQ=255)                                                                   tttCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGaa                                                  >  2:31440/1‑140 (MQ=255)                                                                   tttCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGGGCCAGTCTTGGCATAGTAGAAGTCGATCCTGaa                                                  >  2:1136404/1‑140 (MQ=255)                                                                       aCTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAaacaa                                              <  1:1398982/140‑1 (MQ=255)                                                                       aCTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAaacaa                                              <  1:1497766/140‑1 (MQ=255)                                                                       aCTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAaacaa                                              <  2:1026570/140‑1 (MQ=255)                                                                           aacaacGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAaacaa                                              <  2:540719/136‑1 (MQ=255)                                                                           aacaacGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAaacaa                                              >  1:540719/1‑136 (MQ=255)                                                                            acaacGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGAt                                         <  1:738231/140‑1 (MQ=255)                                                                             caacGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATg                                        >  1:1453511/1‑140 (MQ=255)                                                                              aacGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGc                                       <  2:873518/140‑1 (MQ=255)                                                                              aacGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGc                                       <  2:1292938/140‑1 (MQ=255)                                                                                cGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAg                                     >  2:1098303/1‑140 (MQ=255)                                                                                   gaATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACa                                  <  1:1490510/140‑1 (MQ=255)                                                                                              tttCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCAt                       >  1:1429779/1‑140 (MQ=255)                                                                                                 cgtcgCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCAtgct                    >  2:746629/1‑140 (MQ=255)                                                                                                   tcgCTGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCAtgctgc                  <  2:1453511/140‑1 (MQ=255)                                                                                                       tGCGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGat              >  2:5767/1‑140 (MQ=255)                                                                                                         cGAGAACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGatat            >  2:1126052/1‑140 (MQ=255)                                                                                                            gaACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGc         <  1:1098303/140‑1 (MQ=255)                                                                                                            gaACGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGc         <  2:856950/140‑1 (MQ=255)                                                                                                              aCGCATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGa       <  1:815910/140‑1 (MQ=255)                                                                                                                 cATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCAtgctgc                  <  2:234945/126‑1 (MQ=255)                                                                                                                 cATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCAtgctgc                  >  1:234945/1‑126 (MQ=255)                                                                                                                 cATTCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGt    >  2:1008870/1‑140 (MQ=255)                                                                                                                   ttCAGCAGCAAATACGTGGATAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTAt  >  2:1237677/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                        GCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGAGAATATCACTTAGGTGATGTTGATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTAT  >  NC_000913/1878427‑1878680

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