BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,664,028	C→G	6.9%	L398L (CTG→CTC) ‡	guaA ←	GMP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,664,028	0	C	G	6.9%	321.3 / 10.8	101	L398L (CTG→CTC) ‡	guaA	GMP synthetase
Reads supporting (aligned to +/- strand): ref base C (45/49); new base G (2/5); total (47/54)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 109 total aligned reads displayed.

AGTACTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGC  >  Exported/1663961‑1664097

aGTACTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGc                                                                    <  1:467392/71‑1 (MQ=255)

aGTACTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGc                                                                    <  1:1457215/71‑1 (MQ=255)

    ctctTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGt                                                                <  1:3421603/71‑1 (MQ=255)

      ctttcttCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGtag                                                              >  1:4486477/1‑69 (MQ=255)

      ctttcttCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCg                                                              >  1:1043004/1‑71 (MQ=255)

      ctttcttCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCg                                                              >  1:4744374/1‑71 (MQ=255)

      ctttcttCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCg                                                              >  1:4616777/1‑71 (MQ=255)

         tcttcACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTa                                                            <  1:123529/70‑1 (MQ=255)

          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)

              acttcacCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACgg                                                         <  1:3845360/68‑1 (MQ=255)

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

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

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

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

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

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

                       aGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTc                                             <  1:2266476/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:4475641/44‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                       ggCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTa                                                            <  1:2332097/40‑1 (MQ=255)

                                       ggCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTa                                                            <  1:3573669/40‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                  gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                                    >  1:2825778/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:202627/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                    gAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                 <  1:2531309/70‑1 (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:3250474/71‑3 (MQ=255)

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

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

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

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

                                                          gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt          <  1:2730231/71‑3 (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)

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

                                                              cGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCttt      <  1:1593999/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:2513688/1‑71 (MQ=255)

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

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

                                                                   cAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                              <  1:306198/42‑1 (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:1541885/42‑1 (MQ=255)

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

AGTACTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGC  >  Exported/1663961‑1664097

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