BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,870,791	T→C	V136A (GTT→GCT)	cdgI →	putative c‑di‑GMP binding protein CdgI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,870,791	0	T	C	100.0%	57.1 / NA	17	V136A (GTT→GCT)	cdgI	putative c‑di‑GMP binding protein CdgI
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (15/2); total (15/2)

TCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGAT  >  NC_000913/1870722‑1870858

tCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGCTGAAACCATTg                                                          >  2:330592/1‑81 (MQ=255)

 cACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGCTGAAACCATTg                                                          >  2:330595/1‑80 (MQ=255)

  aCTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGCTGAAACCATTGt                                                         >  2:330597/1‑80 (MQ=255)

  aCTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGCTGAAACCATTGt                                                         >  1:330599/1‑80 (MQ=255)

   cTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGCTGAAACCATTGTaa                                                       >  1:330603/1‑81 (MQ=255)

                     cAACACTCACTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGCTGAAACCATTGTAATTATCCATAAACCAAtt                                      >  1:330605/1‑80 (MQ=255)

                      aaCACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGCTGAAACCATTGTAATTATCCATAAACCAATTaa                                    >  1:330624/1‑81 (MQ=255)

                        cacTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGCTGAAACCATTGTAATTATCCATAAACCAATTAACg                                  >  2:330593/1‑81 (MQ=255)

                            cATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGCTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAg                              >  2:330604/1‑81 (MQ=255)

                                  aaGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGCTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTacac                         >  2:330598/1‑80 (MQ=255)

                                  aaGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGCTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTacac                         >  2:330602/1‑80 (MQ=255)

                                  aaGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGCTGAAACCATTGTAATTAGCCATAAACCAATTAACGGCAGTacac                         >  1:330606/1‑80 (MQ=255)

                                          cATTTTTAAGTGGCCTTATCTATTTAGCTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTAATCCAGa                >  1:330608/1‑81 (MQ=255)

                                          cATTTTTAAGTGGCCTTATCTATTTAGCTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTAATCCAGa                >  2:330600/1‑81 (MQ=255)

                                                      gCCTTATCTATTTAGCTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTAATCCAGACAAAGTCGAatg    <  1:330584/81‑1 (MQ=255)

                                                      gCCTTATCTATTTAGCTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTAATCCAGACAAAGTCGAatg    <  1:330583/81‑1 (MQ=255)

                                                        cTTATCTATTTAGCTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTAATCCAGACAAAGTCGAatgat  >  1:330601/1‑81 (MQ=255)

TCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGAT  >  NC_000913/1870722‑1870858

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