BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	232,635	T→C	8.3%	N43D (AAC→GAC)	hemB ←	porphobilinogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	232,635	0	T	C	8.3%	130.3 / 4.0	48	N43D (AAC→GAC)	hemB	porphobilinogen synthase
Reads supporting (aligned to +/- strand): ref base T (22/22); new base C (2/2); total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

CAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACAT  >  minE/232566‑232704

cAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTg                                                                      >  1:13303/1‑71 (MQ=255)

 aaTTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGt                                                                     >  1:3327146/1‑71 (MQ=255)

  aTTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTc                                                                    <  1:1444981/71‑1 (MQ=255)

   tttCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTc                                                                    >  1:591414/1‑70 (MQ=255)

       gcgTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAAtt                                                                <  1:3206977/70‑1 (MQ=255)

       gcgTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAAtt                                                                <  1:3162463/70‑1 (MQ=255)

            ccAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttcttct                                                          <  1:3294288/71‑1 (MQ=255)

             cAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttc                                                              <  1:448631/66‑1 (MQ=255)

             cAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttc                                                              <  1:1804289/66‑1 (MQ=255)

             cAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttc                                                              <  1:2899318/66‑1 (MQ=255)

              aGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttcttcttc                                                        <  1:2808563/71‑1 (MQ=255)

              aGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATttcttctt                                                         >  1:1379381/1‑70 (MQ=255)

              aGATGTTTCTCTGGAATGCGCAACACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATttcttctt                                                         >  1:3407350/1‑70 (MQ=255)

               gATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttcttcttc                                                        >  1:1658605/1‑70 (MQ=255)

                aTGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaa                                                      >  1:1305333/1‑71 (MQ=255)

                  gTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCa                                                                   <  1:895459/56‑1 (MQ=255)

                     tctcTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaaca                                                    >  1:2645609/1‑68 (MQ=255)

                     tctcTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaaca                                                    >  1:2109814/1‑68 (MQ=255)

                      ctctGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAg                                                 <  1:2986580/70‑1 (MQ=255)

                         tGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaaca                                                    >  1:1968962/1‑64 (MQ=255)

                         tGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGAt                                               <  1:3005569/69‑1 (MQ=255)

                           gAATGCGCATCACGCCTGGCCTGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATc                                              <  1:2311889/68‑1 (MQ=255)

                           gAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGAta                                              <  1:1054542/68‑2 (MQ=255)

                             aTGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaaca                                                    >  1:2476028/1‑60 (MQ=255)

                             aTGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCaa                                         <  1:410571/71‑1 (MQ=255)

                              tGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAac                                        <  1:2714456/71‑1 (MQ=255)

                                  cATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAacac                                      >  1:72958/1‑69 (MQ=255)

                                  cATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAacac                                      >  1:77291/1‑69 (MQ=255)

                                   aTCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGc                                           >  1:2691094/1‑63 (MQ=255)

                                   aTCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGc                                           >  1:1962779/1‑63 (MQ=255)

                                   aTCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAg                                   >  1:1563925/1‑71 (MQ=255)

                                     cACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCTGCa                                          >  1:3626839/1‑62 (MQ=255)

                                     cACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGc                                           >  1:881434/1‑61 (MQ=255)

                                     cACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCa                                          >  1:455136/1‑62 (MQ=255)

                                      aCGCGTGGCATGGCTTCAACGGGTTTGTAGTTGTCAATttc                                                              <  1:2221621/41‑1 (MQ=38)

                                      aCGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttctt                                                            <  1:831238/43‑1 (MQ=255)

                                          cTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTa                            <  1:3062739/71‑1 (MQ=255)

                                            ggCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAg                          >  1:2779618/1‑71 (MQ=255)

                                             gCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAgg                         <  1:2565668/71‑1 (MQ=255)

                                               aTGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCt                       <  1:2429256/71‑1 (MQ=255)

                                                tGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGt                              >  1:3086716/1‑63 (MQ=255)

                                                 ggCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCt                       >  1:2149020/1‑69 (MQ=255)

                                                      cAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACAcc                                     >  1:3110058/1‑50 (MQ=255)

                                                      cAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAgtg                  <  1:2463504/69‑1 (MQ=255)

                                                      cAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAgtg                  <  1:491353/69‑1 (MQ=255)

                                                           gCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGtctc           <  1:2671273/71‑1 (MQ=255)

                                                           gCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGtctc           >  1:3162152/1‑71 (MQ=255)

                                                                  aGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCaaa     <  1:368078/70‑1 (MQ=255)

                                                                    ttGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACa   >  1:2633766/1‑70 (MQ=255)

                                                                    ttGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACAt  <  1:3619622/71‑1 (MQ=255)

CAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACAT  >  minE/232566‑232704

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