Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	188,837	C→T	R316C (CGC→TGC)	PP_0165 →	GGDEF domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	188,837	0	C	T	100.0%	167.6 / NA	51	R316C (CGC→TGC)	PP_0165	GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/22);  total (29/22)

GCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCCGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACGGGGGCCAGTGATATCGACCCGGTGGCTTGTATAGGCCT  >  NC_002947/188698‑188982                                                                                                                                            |                                                                                                                                                  gctgTTGCTGCGGGTTCAGGCCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCgcgg                                                                                                                                         >  1:808486/1‑150 (MQ=255)  ctgTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCgcggc                                                                                                                                        <  1:547248/150‑1 (MQ=255)         tgcGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGt                                                                                                                                 >  2:946158/1‑150 (MQ=255)                   gggCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCt                                                                                                                       <  1:1205360/150‑1 (MQ=255)                            tggcCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCggg                                                                                                              >  1:1008591/1‑150 (MQ=255)                                    aCGCCCGCCTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGc                                                                                                      >  2:996919/1‑150 (MQ=255)                                      gcccgccTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCa                                                                                                                            <  2:1056008/126‑1 (MQ=255)                                      gcccgccTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCa                                                                                                                            >  1:1056008/1‑126 (MQ=255)                                       cccgccTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATg                                                                                                   >  2:399643/1‑150 (MQ=255)                                                   gcCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTac                                                                                       <  1:920980/150‑1 (MQ=255)                                                     cAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTacac                                                                                     >  1:400426/1‑150 (MQ=255)                                                      aGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACAcc                                                                                    <  2:1259211/150‑1 (MQ=255)                                                              gatctgttgttACAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTCCACGGGGGACGTGCATAGGTTGGCCTGCAGCCGCGGCGGCAACTCTACAGTAATGCCCCTCCACAGGGTGTATGAGGAGGCAGTACACCTTGCCCAg                                                                            <  1:1264885/145‑1 (MQ=255)                                                                  aGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCa                                                                             >  1:337422/1‑145 (MQ=255)                                                                  aGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCa                                                                             <  2:337422/145‑1 (MQ=255)                                                                   gttgttGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGctct                                                                       >  2:425225/1‑150 (MQ=255)                                                                    ttgttgCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAggagatc                                                                      >  2:1264885/1‑145 (MQ=255)                                                                      gttgCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACAcc                                                                                    <  1:1161616/134‑1 (MQ=255)                                                                      gttgCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACAcc                                                                                    >  2:1161616/1‑134 (MQ=255)                                                                                gTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTg                                                          <  1:399643/150‑1 (MQ=255)                                                                                  aGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCa                                                        <  1:996919/150‑1 (MQ=255)                                                                                      gAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGc                                                    >  2:292700/1‑150 (MQ=255)                                                                                        aCAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAg                                                                   >  1:403338/1‑133 (MQ=255)                                                                                        aCAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAg                                                                   <  2:403338/133‑1 (MQ=255)                                                                                            cTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCAc                                              >  2:785808/1‑150 (MQ=255)                                                                                                 tcgtACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAAc                                         >  2:244510/1‑150 (MQ=255)                                                                                                   gtACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACgg                                       >  1:43504/1‑150 (MQ=255)                                                                                                     aCTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACgggg                                     <  2:808486/150‑1 (MQ=255)                                                                                                      cTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACggggg                                    <  2:1008591/150‑1 (MQ=255)                                                                                                         gtgCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTg                                                 >  1:751522/1‑134 (MQ=255)                                                                                                         gtgCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTg                                                 <  2:751522/134‑1 (MQ=255)                                                                                                          tgCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACGGGGGCCAg                                >  2:539276/1‑150 (MQ=255)                                                                                                           gtcAGCTACCTGGAAACCAACGACCTCATTTCCTTCAGCCACGGCGTTGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGATGCAACGTTGCAGAGCCTGCACGAAACGGGGGCCAGt                               <  1:540388/148‑1 (MQ=255)                                                                                                             cAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATgag                                                                                                 <  2:845288/82‑1 (MQ=255)                                                                                                             cAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATgag                                                                                                 >  1:845288/1‑82 (MQ=255)                                                                                                               gCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACg                                                             <  2:88158/116‑1 (MQ=255)                                                                                                               gCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACg                                                             >  1:88158/1‑116 (MQ=255)                                                                                                                  aCCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAgg                                                                  <  2:1081248/108‑1 (MQ=255)                                                                                                                  aCCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAgg                                                                  >  1:1081248/1‑108 (MQ=255)                                                                                                                   cccGGCAACCAACGACCTCATTTCCTGCAGCCGCAGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCACACCCTGCACGAAACGGGCCCCAATGATATCgc                       >  2:567086/1‑149 (MQ=255)                                                                                                                   cccGGAAACCCACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGTGGAGTTAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACTTTGCAGAGCCTGCACGAAACGGGGGCCAGTGATATCGa                       <  1:292700/150‑1 (MQ=255)                                                                                                                    ccGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACGGGGGCCAGTGATATCGAc                      >  2:388094/1‑150 (MQ=255)                                                                                                                         aaCCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACGGGGGCCAGTGATATCGACCCGGt                 >  2:417379/1‑150 (MQ=255)                                                                                                                            cAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACGGGGGCCAGTGATATCGACCCGGTGGc              >  2:1262713/1‑150 (MQ=255)                                                                                                                                    tCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGc                                                    <  1:864766/104‑1 (MQ=255)                                                                                                                                    tCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGc                                                    >  2:864766/1‑104 (MQ=255)                                                                                                                                    tCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACgggg                                     >  2:1153698/1‑119 (MQ=255)                                                                                                                                    tCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACgggg                                     <  1:1153698/119‑1 (MQ=255)                                                                                                                                    tCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACGGGGGCCAGTGATATCGACCCGGTGGCTTGTATAg      >  2:152378/1‑150 (MQ=255)                                                                                                                                    tCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACGGGGGCCAGTGATATCGACCCGGTGGCTTGTATAg      >  1:295598/1‑150 (MQ=255)                                                                                                                                        ttCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACGGGGGCCAGTGATATCGACCCGGTGGCTTGTATAGGCCt  <  1:539276/150‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                  GCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCCGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACGGGGGCCAGTGATATCGACCCGGTGGCTTGTATAGGCCT  >  NC_002947/188698‑188982

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