Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,447,556	T→C	E64G (GAA→GGA)	aroK ←	shikimate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,447,556	0	T	C	97.1%	139.8 / ‑5.4	35	E64G (GAA→GGA)	aroK	shikimate kinase I
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (17/17);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTC  >  REL606/3447486‑3447627                                                                       |                                                                        ggATTCCACAGACCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTcg                                                                    <  1:516101/76‑1 (MQ=255)  gATTCCACAGACCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTcgc                                                                   <  2:857872/76‑1 (MQ=255)         caGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATcgcg                                                            >  2:1311353/1‑76 (MQ=255)                  cgccAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTCCGCGATCGCGGAAGGcttc                                                   <  1:1262317/76‑1 (MQ=255)                    ccAGTAGCNAGCGCAATACCCTGTATCGCGGTCGCCTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCcttctt                                                 <  2:242363/76‑1 (MQ=255)                    ccAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCcttctt                                                 <  1:469273/76‑1 (MQ=255)                     cAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCttcttc                                                <  1:723152/76‑1 (MQ=255)                       gTAGCCAGCACANTANCCTGTTTCTCGGTCAACTNATTGATTACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGc                                              <  1:1203015/76‑1 (MQ=255)                       gTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATGGCGGAAGCCTTCTTCGc                                              >  1:748927/1‑76 (MQ=255)                        tAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGcc                                             >  1:1120680/1‑76 (MQ=255)                        tAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGcc                                             <  2:1224722/76‑1 (MQ=255)                         aGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCt                                            >  2:441231/1‑76 (MQ=255)                            cAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCt                                         >  2:417784/1‑76 (MQ=255)                            cAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCt                                         <  1:1089510/76‑1 (MQ=255)                             aGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTa                                        <  2:860598/76‑1 (MQ=255)                              gCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTaa                                       <  1:695236/76‑1 (MQ=255)                                 caATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATc                                    >  2:430283/1‑76 (MQ=255)                                  aaTACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCg                                   >  2:1156813/1‑76 (MQ=255)                                  aaTACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCg                                   >  1:255344/1‑76 (MQ=255)                                  aaTACCCTGTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCg                                   >  1:1152296/1‑76 (MQ=255)                                          gTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCa                           >  1:456960/1‑76 (MQ=255)                                          gTTTCTCGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCa                           >  1:274474/1‑76 (MQ=255)                                             tctcGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGcc                        >  2:676944/1‑76 (MQ=255)                                              ctcGGTCAACTCATTGATGACCTTTCCTTCGCGATCGCAGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGccc                       <  1:1246868/76‑1 (MQ=255)                                                   tCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATc                  <  1:1143374/76‑1 (MQ=255)                                                     aaCTCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAg                <  1:973987/76‑1 (MQ=255)                                                        tCATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTc             <  1:57027/76‑1 (MQ=255)                                                         cATTGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTcc            >  2:836311/1‑76 (MQ=255)                                                           ttGATGACCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACGCAGCCCACATCAGCTCCgg          <  1:1257076/76‑1 (MQ=255)                                                            tgatgaCCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGt         >  1:1145832/1‑76 (MQ=255)                                                             gatgaCCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGtt        >  1:1092643/1‑76 (MQ=255)                                                                gaCCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGt     >  1:837257/1‑76 (MQ=255)                                                                 aCCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGtt    >  1:943484/1‑76 (MQ=255)                                                                 aCCTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGtt    <  1:457600/76‑1 (MQ=255)                                                                   cTTTCCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTtc  <  1:845581/76‑1 (MQ=255)                                                                       |                                                                        GGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTC  >  REL606/3447486‑3447627

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