BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,664,058	A→G	6.8%	G388G (GGT→GGC)	guaA ←	GMP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,664,058	0	A	G	6.8%	296.6 / 9.6	103	G388G (GGT→GGC)	guaA	GMP synthetase
Reads supporting (aligned to +/- strand): ref base A (57/39); new base G (2/5); total (59/44)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

Only 100 of 103 total aligned reads displayed.

AACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCATCT  >  Exported/1663992‑1664122

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

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

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

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

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

                    ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                          <  1:3829704/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)

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

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

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

                           gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                                   >  1:1057193/1‑69 (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:3385336/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                           tCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGa                                         >  1:4037875/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:4562129/1‑55 (MQ=255)

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

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

                                                         cAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCa     >  1:2172272/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)

AACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCATCT  >  Exported/1663992‑1664122

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