BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,196,247	A→G	L375M (TTA→CTG)	icd →	isocitrate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,196,247	0	A	G	95.0%	36.1 / NA	21	L375M (TTA→CTG)	icd	isocitrate dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (16/3); minor base T (0/1); total (17/4)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCT  >  NC_000913/1196119‑1196309

gcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCGATACGGGTTGGACTGAAGCGGCTGACCTGATTGAGAAAGGTAt                                                  >  2:611305/1‑143 (MQ=38)

gcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGGCTGAAGCGGCTGACCTGCTTGTCAAAGGTAt                                                  >  1:702711/1‑143 (MQ=38)

 cccTGTTTGAAGCCTCTCACGTTTCTGCTCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATg                                                 <  1:710364/143‑1 (MQ=25)

  ccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATgg                                                >  1:387224/1‑143 (MQ=255)

      tttGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAgg                                            >  1:92103/1‑143 (MQ=255)

      tttGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGAGCCATATGGGTTGGACTGAAGACGCTGACCTGACTGATACTGGAATGGAAgg                                            >  2:304815/1‑143 (MQ=37)

        tGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCCAAAATTGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcg                                          >  1:200897/1‑143 (MQ=37)

            gCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAACTAAACCCTGGCTCTATTATTCTCTCCGCTGAGATCATGATGCGACACATGGGTTGGACTCAAGCGGCTTACCTCAATGCTAAAGGTATCGAAGGCGCaac                                      >  2:545423/1‑142 (MQ=255)

                          tGCGCCGAAATATGCCGGTCAGGACAAAGTACAACCTGTAGCTAATATTCTCTCCGCTGAGATGATGCTGCGCAATATGGGTTGGACTGACGCGGCTGACCTGATTGTTAAAGGTATGGAAGTAGCAATCAATGCCaagatcg                        >  2:397298/1‑138 (MQ=17)

                          tGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGGTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCg                        >  1:252810/1‑143 (MQ=18)

                            cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACGGAAGCGGCTGACCTGAGTGTTAAGGGTATGGAAGGCGCAATCAGTGCCAAGACCGTa                      >  1:356933/1‑143 (MQ=18)

                             gcctgagtatGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCGCCGCTGAGATGATGCGGCGCCAGAGGGGTTGGACTGCAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTaa                     <  2:287278/136‑1 (MQ=17)

                             gcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTaa                     >  2:680619/1‑143 (MQ=18)

                               cGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTCTGGTAGGCGCAATCAATGCCAAGACCGTAAct                   >  1:353338/1‑142 (MQ=18)

                               cGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATACCAAGACCGTAAct                   >  1:614476/1‑142 (MQ=18)

                                  aaTATGCCGGTCAGGGCAAAGTAAATCCTGGCTCTATCATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTAt                <  1:119872/143‑1 (MQ=14)

                                  aaTATGCCGGTCAGGACAGAGTGAAGCCTGGCTTTATTATGCTCTCCGCTGAGAGGATGCTGAGCCATATGGGTTGTACTGAAGCGTCTGACCTTATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTAt                <  2:265552/143‑1 (MQ=14)

                                      tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATGTGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACt            >  2:432158/1‑143 (MQ=17)

                                        ccGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTc          >  1:418114/1‑143 (MQ=17)

                                          ggTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGGCGCGGCTGACCTGATTGTTAAAGGTATGGAAAGCGCAATAAATGCCCAGACCGTAACTTATGACtacaa        >  2:337384/1‑139 (MQ=17)

                                                gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCt  >  1:791664/1‑143 (MQ=11)

GCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCT  >  NC_000913/1196119‑1196309

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