BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,664,040	C→G	7.1%	P394P (CCG→CCC)	guaA ←	GMP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,664,040	0	C	G	7.1%	309.0 / 9.0	97	P394P (CCG→CCC)	guaA	GMP synthetase
Reads supporting (aligned to +/- strand): ref base C (45/45); new base G (2/5); total (47/50)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

CTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAAC  >  minE/1663971‑1664105

cttcacttcaCCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACTGTACAGCATGTCGTAc                                                                   <  1:3064787/70‑1 (MQ=255)

cttcacttcaCCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                                   <  1:1048961/70‑1 (MQ=255)

cttcacttcaCCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                                   <  1:528721/70‑1 (MQ=255)

   cacttcacCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACg                                                                  <  1:88332/68‑1 (MQ=255)

   cacttcacCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACg                                                                  <  1:1661854/68‑1 (MQ=255)

       tcacCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGcc                                                           >  1:3459166/1‑71 (MQ=255)

         acCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCa                                                         <  1:2388972/71‑1 (MQ=255)

            cAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCt                                                      >  1:3896723/1‑71 (MQ=255)

            cAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCt                                                      >  1:4179993/1‑71 (MQ=255)

            cAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCt                                                      >  1:1198864/1‑71 (MQ=255)

             aGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTc                                                     <  1:2266476/71‑1 (MQ=255)

             aGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTc                                                     <  1:1804236/71‑1 (MQ=255)

             aGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTc                                                     <  1:3074858/71‑1 (MQ=255)

                acacGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGGCGTACGGCAGGCCCAGCTCCa                                                   <  1:2878320/70‑1 (MQ=255)

                     aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCaa                                             >  1:1413873/1‑71 (MQ=255)

                       cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCa                                                               <  1:1850033/51‑1 (MQ=255)

                          cAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                                   <  1:2084574/44‑1 (MQ=255)

                          cAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                                   <  1:1933782/44‑1 (MQ=255)

                          cAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                                   <  1:4475641/44‑1 (MQ=255)

                          cAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                                   <  1:5090592/44‑1 (MQ=255)

                          cAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                                   <  1:821432/44‑1 (MQ=255)

                          cAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGGCGTAc                                                                   <  1:4117009/44‑1 (MQ=255)

                          cAAGGCCTGGGCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                                   <  1:2522327/44‑1 (MQ=255)

                            aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTAc                                       >  1:1754985/1‑70 (MQ=255)

                            aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTAc                                       >  1:1024488/1‑70 (MQ=255)

                            aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                      >  1:4000773/1‑71 (MQ=255)

                            aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                      >  1:4037499/1‑71 (MQ=255)

                            aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                      >  1:3997295/1‑71 (MQ=255)

                            aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                      >  1:366093/1‑71 (MQ=255)

                            aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                      >  1:2832152/1‑71 (MQ=255)

                            aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                      >  1:1214456/1‑71 (MQ=255)

                                cTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACt                                  <  1:3449006/71‑1 (MQ=255)

                                      ccGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCa                                              <  1:2655393/53‑1 (MQ=255)

                                      ccGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCa                                              <  1:343356/53‑1 (MQ=255)

                                      ccGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCa                                              <  1:2031173/53‑1 (MQ=255)

                                      ccGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCa                                              <  1:3514854/53‑1 (MQ=255)

                                      ccGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCa                                              <  1:1735353/53‑1 (MQ=255)

                                      ccGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCa                                              <  1:4936268/53‑1 (MQ=255)

                                        gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACGAAt                                            >  1:240993/1‑53 (MQ=255)

                                        gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                                            >  1:3406977/1‑53 (MQ=255)

                                        gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                                            >  1:3103929/1‑53 (MQ=255)

                                        gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                                            >  1:57622/1‑53 (MQ=255)

                                        gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                                            >  1:806880/1‑53 (MQ=255)

                                        gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                                            >  1:2825778/1‑53 (MQ=255)

                                        gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                                            >  1:2768266/1‑53 (MQ=255)

                                        gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                                            >  1:1217841/1‑53 (MQ=255)

                                        gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                                            >  1:2622858/1‑53 (MQ=255)

                                        gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                                            >  1:1558924/1‑53 (MQ=255)

                                        gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                                            >  1:2322804/1‑53 (MQ=255)

                                        gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                                            >  1:2270316/1‑53 (MQ=255)

                                        gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCttt                          >  1:4807250/1‑71 (MQ=255)

                                        gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCtt                           >  1:1396368/1‑70 (MQ=255)

                                         ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCttt                          >  1:244365/1‑70 (MQ=255)

                                         ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                         <  1:1249020/71‑1 (MQ=255)

                                         ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                         <  1:131728/71‑1 (MQ=255)

                                         ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                         <  1:4388070/71‑1 (MQ=255)

                                         ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                         <  1:3829704/71‑1 (MQ=255)

                                         ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                         <  1:3781473/71‑1 (MQ=255)

                                         ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                         <  1:3679304/71‑1 (MQ=255)

                                         ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                         <  1:202627/71‑1 (MQ=255)

                                          gAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                         >  1:304776/1‑70 (MQ=255)

                                           aaCGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGAc                                                >  1:3666372/1‑46 (MQ=255)

                                           aaCGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGaa                       >  1:1449930/1‑71 (MQ=255)

                                           aaCGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGaa                       >  1:1648502/1‑71 (MQ=255)

                                                gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                  <  1:5014705/71‑3 (MQ=255)

                                                gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                  <  1:118674/71‑3 (MQ=255)

                                                gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                  <  1:4002955/71‑3 (MQ=255)

                                                gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                  <  1:3250474/71‑3 (MQ=255)

                                                gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                  >  1:3863788/1‑69 (MQ=255)

                                                gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                  <  1:2730231/71‑3 (MQ=255)

                                                gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                  >  1:1057193/1‑69 (MQ=255)

                                                gTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAAAAGct                  >  1:3510697/1‑71 (MQ=255)

                                                 tGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCa                                              <  1:3795137/42‑1 (MQ=255)

                                                   aCGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCa                                    >  1:1793439/1‑50 (MQ=255)

                                                   aCGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCAc                                   >  1:689712/1‑51 (MQ=255)

                                                   aCGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGctct                <  1:2207894/70‑1 (MQ=255)

                                                    cGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCttt              <  1:1593999/71‑1 (MQ=255)

                                                    cGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCttt              <  1:2413908/71‑1 (MQ=255)

                                                     ggTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGaa                       >  1:667252/1‑61 (MQ=255)

                                                     ggTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCttt              >  1:3410301/1‑70 (MQ=255)

                                                      gTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCa            <  1:1852451/71‑1 (MQ=255)

                                                       tACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCt                                          <  1:142339/40‑1 (MQ=255)

                                                        aCAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGc          >  1:1380379/1‑71 (MQ=255)

                                                        aCAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGc          >  1:2513688/1‑71 (MQ=255)

                                                         cAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                      <  1:5126886/42‑1 (MQ=255)

                                                         cAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                      <  1:1772229/42‑1 (MQ=255)

                                                         cAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                      <  1:1338403/42‑1 (MQ=255)

                                                         cAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                      <  1:485438/42‑1 (MQ=255)

                                                         cAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                      <  1:1541885/42‑1 (MQ=255)

                                                         cAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                      <  1:3955767/42‑1 (MQ=255)

                                                         cAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                      <  1:306198/42‑1 (MQ=255)

                                                         cAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAACCTTACg                                      <  1:2766882/42‑1 (MQ=255)

                                                             aTGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTc     <  1:4370710/71‑1 (MQ=255)

                                                              tGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCa    <  1:4073385/71‑1 (MQ=255)

                                                              tGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCa    <  1:1792390/71‑1 (MQ=255)

                                                              tGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCa    <  1:2596671/71‑1 (MQ=255)

                                                                tCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGa                        >  1:4037875/1‑49 (MQ=255)

                                                                tCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGa                        >  1:3474764/1‑49 (MQ=255)

                                                                tCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAAc  >  1:1757827/1‑71 (MQ=255)

                                                                tCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAAAAGCTCTTTCAGCGGTTCAAc  >  1:4832436/1‑71 (MQ=255)

CTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAAC  >  minE/1663971‑1664105

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