BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	753,415	A→G	Y352H (TAC→CAC)	gltA ←	citrate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	753,415	0	A	G	96.2%	79.9 / ‑3.8	26	Y352H (TAC→CAC)	gltA	citrate synthase
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (13/12); total (14/12)
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.52e-01

AAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTCCAGCAGGTCATCCTTCGTGCCCAGCTCTTTCAG  >  NC_000913/753337‑753490

aaGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTGCg                                                                           >  2:160186/1‑81 (MQ=255)

 aGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTGCgg                                                                          <  2:160155/81‑1 (MQ=255)

 aGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTGCgg                                                                          <  1:160154/81‑1 (MQ=255)

     ggAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTGCGGGTCg                                                                      >  1:160189/1‑81 (MQ=255)

        aTACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTGCGGGTCGTTc                                                                   <  2:160156/81‑1 (MQ=255)

         tACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTGCGGGTCGTTCa                                                                  <  1:160157/81‑1 (MQ=255)

            ccATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTGCGGGTCGTTCAgcg                                                               >  1:160222/1‑81 (MQ=255)

               tCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTGCGGGTCGTTCAGCGCGa                                                            <  1:160158/81‑1 (MQ=255)

                 gCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTGCGGGTCGTTCAGCGCGATg                                                          >  1:160198/1‑81 (MQ=255)

                       aGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTGCGGGTCGTTCAGCGCGATGTTTTc                                                     >  2:160199/1‑80 (MQ=255)

                               aCCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCGTTCAGCGCGATGTTTTCCAGCTCCa                                             >  2:160197/1‑80 (MQ=255)

                                         aaaTCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTGCGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTcc                                  <  2:160159/81‑1 (MQ=255)

                                          aaTCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTGCGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTCca                                 <  1:160160/81‑1 (MQ=255)

                                               aCGTTCGGGTACAGTTTCTTCTCGATAAAGTGCGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTCCAGCAgg                            >  1:160206/1‑81 (MQ=255)

                                                   tCGGGTACAGTTTCTTCTCGATAAAGTGCGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTCCAGCAGGTCAt                        >  1:160190/1‑81 (MQ=255)

                                                    cGGGTACAGTTTCTTCTCGATAAAGTGCGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTCCAGCAGGTCATc                       >  2:160195/1‑81 (MQ=255)

                                                         aCAGTTTCTTCTCGATAAAGTGCGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTCCAGCAGGTCATCCTTCg                  >  1:160202/1‑81 (MQ=255)

                                                           aGTTTCTTCTCGATAAAGTGCGGGTCGTTCAGCGCGATGTTTTCCTGCTCCATAGCCACTTCCAGCAGGTCATCCTTCGTg                >  1:160208/1‑81 (MQ=255)

                                                                   ctcGATAAAGTGCGGGTCGTTCAGCGCGGTGTTTTCCAGCTCCATAGCCACTTCCAGCAGGTCATCCTTCGTGCCCAg           <  2:160161/78‑1 (MQ=255)

                                                                   ctcGATAAAGTGCGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTCCAGCAGGTCATCCTTCGTGCCCAGctc        >  1:160203/1‑81 (MQ=255)

                                                                      gATAAAGTGCGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTCCAGCAGGTCATCCTTCGTGCCCAGCTCtt      <  2:160162/80‑1 (MQ=255)

                                                                      gATAAAGTGCGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTCCAGCAGGTCATCCTTCGTGCCCAGCTCtt      <  2:160164/80‑1 (MQ=255)

                                                                      gATAAAGTGCGGGTCGTTCAGCGCGATGTATTCCAGCTCCATAGCCACTTCCAGCAGGTCATCCTTCGTGCCCAGctct       <  2:160163/79‑1 (MQ=255)

                                                                       aTAAAGTGCGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTCCAGCAGGTCATCCTTCGTGCCCAGCTCTTTc    <  1:160165/81‑1 (MQ=255)

                                                                         aaaGTGCGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTCCAGCAGGTCATCCTTCGTGCCCAGCTCTTTCAg  >  1:160205/1‑81 (MQ=255)

                                                                         aaaGTGCGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTCCAGCAGGTCATCCTTCGTGCCCAGCTCTTTCAg  >  2:160200/1‑81 (MQ=255)

AAGACGGAATACCCATCGCTTTCAGGATGATACCAGAGTAGAAATCGACGTTCGGGTACAGTTTCTTCTCGATAAAGTACGGGTCGTTCAGCGCGATGTTTTCCAGCTCCATAGCCACTTCCAGCAGGTCATCCTTCGTGCCCAGCTCTTTCAG  >  NC_000913/753337‑753490

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