Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,664,067	A→G	8.1%	R385R (CGT→CGC)	guaA ←	GMP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,664,067	0	A	G	8.1%	252.4 / 10.2	86	R385R (CGT→CGC)	guaA	GMP synthetase
Reads supporting (aligned to +/- strand):  ref base A (43/36);  new base G (2/5);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

AGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCATCTCTTTCGGC  >  Exported/1663999‑1664130                                                                     |                                                                aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTAc                                                                >  1:1754985/1‑70 (MQ=255) aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTAc                                                                >  1:1024488/1‑70 (MQ=255) aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                                               >  1:4037499/1‑71 (MQ=255) aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                                               >  1:1214456/1‑71 (MQ=255) aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                                               >  1:4000773/1‑71 (MQ=255) aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                                               >  1:2832152/1‑71 (MQ=255) aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                                               >  1:3997295/1‑71 (MQ=255) aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                                               >  1:366093/1‑71 (MQ=255)     cTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACt                                                           <  1:3449006/71‑1 (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:1895575/71‑1 (MQ=255)              ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                  <  1:202627/71‑1 (MQ=255)              ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                  <  1:3679304/71‑1 (MQ=255)              ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                  <  1:3781473/71‑1 (MQ=255)              ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                  <  1:3829704/71‑1 (MQ=255)              ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                  <  1:131728/71‑1 (MQ=255)              ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                  <  1:1249020/71‑1 (MQ=255)              ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                  <  1:4388070/71‑1 (MQ=255)              ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                  <  1:4500411/71‑1 (MQ=255)               gAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                  <  1:2531309/70‑1 (MQ=255)               gAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                  >  1:304776/1‑70 (MQ=255)                aaCGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGaa                                                >  1:1449930/1‑71 (MQ=255)                aaCGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGaa                                                >  1:1648502/1‑71 (MQ=255)                     gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                                           >  1:1057193/1‑69 (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:3863788/1‑69 (MQ=255)                     gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                                           <  1:3250474/71‑3 (MQ=255)                     gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                                           <  1:2730231/71‑3 (MQ=255)                     gTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAAAAGct                                           >  1:3510697/1‑71 (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: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)                             aCAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGc                                   >  1:1380379/1‑71 (MQ=255)                             aCAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGc                                   >  1:3385336/1‑71 (MQ=255)                             aCAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGc                                   >  1:2513688/1‑71 (MQ=255)                              cAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                                               <  1:1338403/42‑1 (MQ=255)                              cAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                                               <  1:1772229/42‑1 (MQ=255)                              cAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                                               <  1:485438/42‑1 (MQ=255)                              cAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                                               <  1:3955767/42‑1 (MQ=255)                              cAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                                               <  1:5126886/42‑1 (MQ=255)                              cAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                                                               <  1:1541885/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:2596671/71‑1 (MQ=255)                                   tGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCa                             <  1:1792390/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:2396171/1‑55 (MQ=255)                                              aGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCg                                  >  1:4174366/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:3089686/1‑71 (MQ=255)                                                   cAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCa             >  1:2172272/1‑71 (MQ=255)                                                   cAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCa             >  1:118300/1‑71 (MQ=255)                                                   cAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCa             >  1:2849757/1‑71 (MQ=255)                                                     gCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCtt                                        >  1:122214/1‑42 (MQ=255)                                                     gCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCAtc           >  1:4273777/1‑71 (MQ=255)                                                     gCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCAtc           >  1:3775184/1‑71 (MQ=255)                                                     gCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCAtc           >  1:2425185/1‑71 (MQ=255)                                                     gCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCAt            >  1:1360736/1‑70 (MQ=255)                                                       tCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCAtct          >  1:839441/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)                                                                 tCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGcc                     <  1:4463593/49‑1 (MQ=255)                                                                 tCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGcc                     <  1:435089/49‑1 (MQ=255)                                                                 tCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGcc                     <  1:3794836/49‑1 (MQ=255)                                                                 tCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGcc                     <  1:375453/49‑1 (MQ=255)                                                                 tCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGcc                     <  1:3168740/49‑1 (MQ=255)                                                                     |                                                                AGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCATCTCTTTCGGC  >  Exported/1663999‑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: ‑