Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,664,079	T→C	10.3%	K381K (AAA→AAG)	guaA ←	GMP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,664,079	0	T	C	10.3%	188.5 / 10.6	68	K381K (AAA→AAG)	guaA	GMP synthetase
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base C (2/5);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCATCTCTTTCGGCAGGCCGCCCACGTTGT  >  Exported/1664011‑1664146                                                                     |                                                                    gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCttt                                                                   >  1:4807250/1‑71 (MQ=255) gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCtt                                                                    >  1:1396368/1‑70 (MQ=255)  ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCttt                                                                   >  1:244365/1‑70 (MQ=255)  ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                                  <  1:4500411/71‑1 (MQ=255)  ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                                  <  1:4388070/71‑1 (MQ=255)  ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                                  <  1:3679304/71‑1 (MQ=255)  ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                                  <  1:1249020/71‑1 (MQ=255)  ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                                  <  1:131728/71‑1 (MQ=255)  ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                                  <  1:202627/71‑1 (MQ=255)  ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                                  <  1:3781473/71‑1 (MQ=255)  ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                                  <  1:1895575/71‑1 (MQ=255)  ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg                                                                  <  1:3829704/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:2730231/71‑3 (MQ=255)         gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                                                           <  1:4002955/71‑3 (MQ=255)         gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                                                           <  1:5014705/71‑3 (MQ=255)         gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                                                           <  1:3250474/71‑3 (MQ=255)         gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                                                           >  1:3863788/1‑69 (MQ=255)         gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                                                           <  1:118674/71‑3 (MQ=255)         gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt                                                           >  1:1057193/1‑69 (MQ=255)         gTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAAAAGct                                                           >  1:3510697/1‑71 (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:2513688/1‑71 (MQ=255)                 aCAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGc                                                   >  1:1380379/1‑71 (MQ=255)                 aCAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGc                                                   >  1:3385336/1‑71 (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: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:4562129/1‑55 (MQ=255)                                  aGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCg                                                  >  1:4174366/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:3089686/1‑71 (MQ=255)                                       cAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCa                             >  1:118300/1‑71 (MQ=255)                                       cAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCa                             >  1:2172272/1‑71 (MQ=255)                                       cAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCa                             >  1:2849757/1‑71 (MQ=255)                                         gCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCtt                                                        >  1:122214/1‑42 (MQ=255)                                         gCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCAtc                           >  1:2425185/1‑71 (MQ=255)                                         gCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCAtc                           >  1:4273777/1‑71 (MQ=255)                                         gCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCAtc                           >  1:3775184/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)                                                              cTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCATCTCTTTc                     <  1:4511355/56‑1 (MQ=255)                                                                tCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCATCTCTTTCGGCAGgcc             <  1:2070448/62‑1 (MQ=255)                                                                 cGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCATCTCtt                       >  1:3127307/1‑51 (MQ=255)                                                                  gTCTTTGAACAGCTCTTTCAGCGGTTCATCCAGGCCCATCTTCATCTCTTTCGGCAGGCCGCCCACGTtgt  <  1:469488/71‑1 (MQ=255)                                                                     |                                                                    GGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCATCTCTTTCGGCAGGCCGCCCACGTTGT  >  Exported/1664011‑1664146

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