BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,196,220	C→T	H366H (CAC→CAT)	icd →	isocitrate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,196,220	0	C	T	97.0%	72.4 / NA	33	H366H (CAC→CAT)	icd	isocitrate dehydrogenase
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (17/15); minor base G (0/1); total (17/16)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGAT  >  NC_000913/1196080‑1196312

atcggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCACGGCTCTATGACTCTCTTCGGTTCTATGGTGGTGCGCCatat                                                                                            >  1:1165230/1‑143 (MQ=255)

     tatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGtt                                                                                       <  1:1131143/143‑1 (MQ=255)

     tatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGtt                                                                                       <  1:556856/143‑1 (MQ=255)

       tcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTgg                                                                                     >  2:898247/1‑143 (MQ=255)

                 tgtaaaCATGGGTGACGAATGCGCCCTGGTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCCGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCgg                                                                           <  1:894134/140‑1 (MQ=255)

                 tGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCgg                                                                           <  2:374812/143‑1 (MQ=255)

                    aaaCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTGTTCTCTCCGCTGGGGTGATGCTGCGACATATGGGGTGGACTGGAGCGGCTg                                                                        >  1:9734/1‑143 (MQ=255)

                      aCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAc                                                                      <  2:95730/143‑1 (MQ=255)

                        aTCGGTGACGAATGCGTACGGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAcct                                                                    <  1:680672/143‑3 (MQ=38)

                              gACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGt                                                              <  2:242443/143‑1 (MQ=38)

                                       gtccggTTTGAAGCCACCCGCGGGACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTAt                                                     <  2:990364/138‑1 (MQ=25)

                                       gcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCGGGACAAAGTAAATCATCGCTCTATTATTCTCTCCGCTGAGATGAGGCTGCGCCATATGGGTTGGACTGAAGCGACTGACCTGATTGGGAAAGGCAt                                                     >  1:204971/1‑143 (MQ=37)

                                         ccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGAATGTTAAAGATATgg                                                   >  2:320292/1‑143 (MQ=37)

                                         ccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTAATATTCACTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATgg                                                   >  2:327413/1‑143 (MQ=37)

                                             tttgtttccacccacGGTGCTGCGGCGAAATAGGCCGGTCAGGACAAAGGAAAGCCTGGCTCTATGATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAgg                                               <  1:99086/136‑1 (MQ=21)

                                                      acccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGGTGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATcca                                      >  2:1025803/1‑141 (MQ=25)

                                                       cccccGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCAGATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaat                                     <  1:274424/143‑2 (MQ=25)

                                                       cccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGCGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaat                                     <  1:182569/143‑2 (MQ=25)

                                                        ccacGGTACTGCGCCGAAATATGTCGGTCAGGACAGAGTAAATCCTGGCTCTATTATTCGCTCCGCTGAGATGCTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaatg                                    <  2:477533/143‑3 (MQ=21)

                                                        ccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaatg                                    >  1:38774/1‑141 (MQ=25)

                                                         cacGGTACTGCGCCGAAATATGCCGGTCAGGACACAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAAtgc                                   <  1:641/143‑1 (MQ=21)

                                                          acGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAACGGTATGGAAGGCGCACTCAAtgcc                                  >  1:1141906/1‑142 (MQ=25)

                                                            ggTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGACGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTCAACGGTATGGAAGGCGCAATCAATGCCaa                                >  2:903520/1‑143 (MQ=21)

                                                                   cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGATATGGAAGGCGCAATCAATGCCAAGACCGTa                         >  2:52179/1‑143 (MQ=18)

                                                                     ccGAAATATGCCGGTCAGGACAAAGTCAATCCTGGCTCTATGATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAAc                       <  2:6212/143‑1 (MQ=17)

                                                                      cGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCCATGCCAAGACCGTAAct                      >  1:421933/1‑142 (MQ=18)

                                                                      cGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAAct                      >  1:250556/1‑142 (MQ=18)

                                                                        aaaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTa                    <  1:517158/143‑1 (MQ=17)

                                                                              gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCCATCAATGCCAAGACCGTAACTTATGACtt              >  2:288567/1‑143 (MQ=17)

                                                                                  gTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGaa          >  1:1243610/1‑142 (MQ=14)

                                                                                   tCAGGGCAAAGTAAATCCTGGCTCTATTATGCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGaac         <  2:519055/143‑3 (MQ=11)

                                                                                    cAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACg        >  1:522953/1‑143 (MQ=14)

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

                                                                                          aaaGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATCTGGGTTGGACTGAAGCGGCTGACCTGACTGTTAAAGGTATGGAAGGCGCAATCAATGCCACGACCGTGACTTATCACTTCGAACGTCTgat  >  2:56187/1‑143 (MQ=14)

ATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGAT  >  NC_000913/1196080‑1196312

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