BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,267,217	G→T	M98I (ATG→ATT)	ychQ →	SirB family protein YchQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,267,217	0	G	T	92.6%	77.1 / 1.0	27	M98I (ATG→ATT)	ychQ	SirB family protein YchQ
Reads supporting (aligned to +/- strand): ref base G (1/1); new base T (11/14); total (12/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GCGGCATTGCGTTGATGGCTAAAGCGCACATCCTGCCATTTTCCGGGCAGGCACAGTGGCTGACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATGCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTACATCATCATTAAACTCGCCACCACAAAAGTACCGTTACTGGGGTAAG  >  NC_000913/1267087‑1267317

gCGGCATTGCATTGATGGCTAAAGCGCACATCCTGCCATTTTCCGGGCAGGCACAGTGGCTGACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTca                                                                                                    >  2:124786/1‑133 (MQ=255)

  ggCATTGCGTTGATGGCTAAAGCGCACATCCTGCCATTTTCCGGGCAGGCACAGTGGCTGACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTcaca                                                                                                  <  1:333522/133‑1 (MQ=255)

      ttGCGTTGATGGCTAAAGCGCACATCCTGCCATTTTCCGGGCAGGCACAGTGGCTGACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTca                                                                                              >  2:284062/1‑133 (MQ=255)

             gATGGCTAAAGCGCACATCCTGCCATTTTCCGGGCAGGCACAGTGGCTGACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGgcg                                                                                       <  1:420635/133‑1 (MQ=255)

              aTGGCTAAAGCGCACATCCTGCCATTTTCCGGGCAGGCACAGTGGCTGACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGgcgc                                                                                      <  2:323731/133‑1 (MQ=255)

                          cacaTCCTGCCATTTTCCGGGCAGGCACAGTGGCTGACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCATATTCACAGTCAGCAGACGCGAATTATTGCCt                                                                          >  2:219873/1‑133 (MQ=255)

                           acaTCCTGCCATTTTCCGGGCAGGCACAGTGGCTGACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCtt                                                                         >  2:26026/1‑133 (MQ=255)

                               ccTGCCATATTCCGGTCATGCACAGTGGCTTACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCg                                                                     <  1:244076/133‑1 (MQ=255)

                                       tttCCGGGCAGGCACAGTGGCTGACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGtt                                                             >  1:110088/1‑133 (MQ=255)

                                         tCCGGGCAGGCACAGTGGCTGACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCGCAGTCAGCAGGCACGCATTATTGCCATCCCTCTGGCGTtgg                                                           >  2:220575/1‑133 (MQ=255)

                                           cGGGCAGGCACAGTGGCTGACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTtggtg                                                         >  2:83209/1‑133 (MQ=255)

                                              gcaggcaCAGTGGCTGACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGgtggtg                                                      <  2:344684/133‑1 (MQ=255)

                                                 ggcaCAGTGGCTGACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTg                                                   <  1:284062/133‑1 (MQ=255)

                                                      aGTGGCTGACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTAcat                                              <  2:110088/133‑1 (MQ=255)

                                                            tgactgaAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTAcatcatcat                                        >  1:274776/1‑133 (MQ=255)

                                                            tgactgaAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTAcatcatcat                                        >  2:504071/1‑133 (MQ=255)

                                                            tgactgaAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATGCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTAcat                                              >  2:261811/1‑127 (MQ=255)

                                                            tgactgaAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATGCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTAcat                                              <  1:261811/127‑1 (MQ=255)

                                                                        tGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTACATCATCATTAAACTCGccac                            <  2:55752/133‑1 (MQ=255)

                                                                        tGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTACATCATCATTAAACTCGccac                            <  2:283035/133‑1 (MQ=255)

                                                                                     cATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTACATCATCATTAAACTCGCCACCACAAAAGTACCg               <  1:504071/133‑1 (MQ=255)

                                                                                         tataTCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTACATCATCATTAAACTCGCCACCACAAAAGTACCGTTAc           <  1:388510/133‑1 (MQ=255)

                                                                                         tataTCGTATTGGGTGTTATTGCACTCGATGATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTACATCATCATTAAACTCGCCACCACAAAAGTACCGTTAc           <  1:219873/133‑1 (MQ=255)

                                                                                               gTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTACATCATCATTAAACTCGCCACCACAAAAGTACCGTTACTGGGGt     <  2:423442/133‑1 (MQ=255)

                                                                                               gTTTTGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTACATCATCATTAAACTCGCCACCACAAAAGTACCGTTACTGGGGt     <  2:382009/133‑1 (MQ=255)

                                                                                                ttttGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTACATCATCATTAAACTCGCCACCACAAAAGTACCGTTACTGGGGTa    >  1:424460/1‑133 (MQ=255)

                                                                                                  ttGGGTTTTATTGCACTCGATTATCGTCGTATTCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTACATCACCATTAAACACGCTACCACAAAAGTAACGTTGCTGGGGTAAg  >  2:162144/1‑133 (MQ=255)

GCGGCATTGCGTTGATGGCTAAAGCGCACATCCTGCCATTTTCCGGGCAGGCACAGTGGCTGACTGAAAAGCTGTTTGGAGTTATCATTTATATCGTTTTGGGTTTTATTGCACTCGATTATCGTCGTATGCACAGTCAGCAGGCGCGCATTATTGCCTTCCCGCTGGCGTTGGTGGTGCTGTACATCATCATTAAACTCGCCACCACAAAAGTACCGTTACTGGGGTAAG  >  NC_000913/1267087‑1267317

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