BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,664,061	A→G	7.4%	I387I (ATT→ATC)	guaA ←	GMP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,664,061	0	A	G	7.4%	275.7 / 9.3	94	I387I (ATT→ATC)	guaA	GMP synthetase
Reads supporting (aligned to +/- strand): ref base A (55/32); new base G (2/5); total (57/37)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

AACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCATCTCTTTCGGC  >  minE/1663992‑1664130

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

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

       aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTAc                                                                >  1:1754985/1‑70 (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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                   gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                                                                     >  1:1217841/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)

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

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

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

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

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

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

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

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

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

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

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

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

                     gAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                  <  1:2531309/70‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                           gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                                           <  1:118674/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)

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

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

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

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

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

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

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

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

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

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

                                         tGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCa                             <  1:1792390/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)

                                                  gCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGccc                    <  1:1088964/71‑1 (MQ=255)

                                                   cAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGccc                    <  1:3178043/70‑1 (MQ=255)

                                                   cAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGcc                     >  1:4670816/1‑69 (MQ=255)

                                                    aGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCg                                  >  1:4174366/1‑55 (MQ=255)

                                                    aGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCg                                  >  1:2396171/1‑55 (MQ=255)

                                                    aGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCg                                  >  1:4562129/1‑55 (MQ=255)

                                                      gCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCt                <  1:4625295/71‑1 (MQ=255)

                                                         cAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTc              >  1:306222/1‑70 (MQ=255)

                                                         cAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCa             >  1:2849757/1‑71 (MQ=255)

                                                         cAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCa             >  1:2172272/1‑71 (MQ=255)

                                                         cAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCa             >  1:118300/1‑71 (MQ=255)

                                                         cAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCa             >  1:3089686/1‑71 (MQ=255)

                                                           gCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCtt                                        >  1:122214/1‑42 (MQ=255)

                                                           gCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCAtc           >  1:2425185/1‑71 (MQ=255)

                                                           gCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCAtc           >  1:3775184/1‑71 (MQ=255)

                                                           gCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCAtc           >  1:4273777/1‑71 (MQ=255)

                                                           gCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCAt            >  1:1360736/1‑70 (MQ=255)

                                                             tCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCAtct          >  1:2933719/1‑70 (MQ=255)

                                                             tCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCACATCTTCAtct          >  1:4144983/1‑70 (MQ=255)

                                                                    cAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCATCTCTTTCGGc  <  1:644616/71‑1 (MQ=255)

AACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCATCTCTTTCGGC  >  minE/1663992‑1664130

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