Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,664,022	A→G	7.8%	R400R (CGT→CGC)	guaA ←	GMP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,664,022	0	A	G	7.8%	263.1 / 9.7	89	R400R (CGT→CGC)	guaA	GMP synthetase
Reads supporting (aligned to +/- strand):  ref base A (38/44);  new base G (2/5);  total (40/49)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

AGTACTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCT  >  minE/1663961‑1664091                                                              |                                                                      aGTACTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGc                                                              <  1:1457215/71‑1 (MQ=255) aGTACTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGc                                                              <  1:467392/71‑1 (MQ=255)     ctctTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGt                                                          <  1:3421603/71‑1 (MQ=255)       ctttcttCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGtag                                                        >  1:4486477/1‑69 (MQ=255)       ctttcttCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCg                                                        >  1:1043004/1‑71 (MQ=255)       ctttcttCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCg                                                        >  1:4616777/1‑71 (MQ=255)       ctttcttCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCg                                                        >  1:4744374/1‑71 (MQ=255)          tcttcACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTa                                                      <  1:123529/70‑1 (MQ=255)           cttcacttcaCCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACTGTACAGCATGTCGTAc                                                     <  1:3064787/70‑1 (MQ=255)           cttcacttcaCCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                     <  1:528721/70‑1 (MQ=255)           cttcacttcaCCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                     <  1:1048961/70‑1 (MQ=255)              cacttcacCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACg                                                    <  1:88332/68‑1 (MQ=255)              cacttcacCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACg                                                    <  1:1661854/68‑1 (MQ=255)               acttcacCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACgg                                                   <  1:3845360/68‑1 (MQ=255)                  tcacCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGcc                                             >  1:3459166/1‑71 (MQ=255)                    acCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCa                                           <  1:2388972/71‑1 (MQ=255)                       cAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCt                                        >  1:4179993/1‑71 (MQ=255)                       cAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCt                                        >  1:3896723/1‑71 (MQ=255)                       cAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCt                                        >  1:908798/1‑71 (MQ=255)                       cAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCt                                        >  1:1198864/1‑71 (MQ=255)                        aGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTc                                       <  1:1804236/71‑1 (MQ=255)                        aGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTc                                       <  1:2266476/71‑1 (MQ=255)                        aGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTc                                       <  1:3074858/71‑1 (MQ=255)                           acacGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGGCGTACGGCAGGCCCAGCTCCa                                     <  1:2878320/70‑1 (MQ=255)                                aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCaa                               >  1:1413873/1‑71 (MQ=255)                                  cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCa                                                 <  1:1850033/51‑1 (MQ=255)                                     cAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                     <  1:821432/44‑1 (MQ=255)                                     cAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                     <  1:5090592/44‑1 (MQ=255)                                     cAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                     <  1:1933782/44‑1 (MQ=255)                                     cAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                     <  1:4475641/44‑1 (MQ=255)                                     cAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                     <  1:2084574/44‑1 (MQ=255)                                     cAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGGCGTAc                                                     <  1:4117009/44‑1 (MQ=255)                                     cAAGGCCTGGGCCCGGGAACGGGTGACGGTACAGCATGTCGTAc                                                     <  1:2522327/44‑1 (MQ=255)                                       aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTAc                         >  1:1754985/1‑70 (MQ=255)                                       aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTAc                         >  1:1024488/1‑70 (MQ=255)                                       aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                        >  1:3997295/1‑71 (MQ=255)                                       aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                        >  1:4000773/1‑71 (MQ=255)                                       aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                        >  1:1214456/1‑71 (MQ=255)                                       aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                        >  1:4037499/1‑71 (MQ=255)                                       aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                        >  1:366093/1‑71 (MQ=255)                                       aGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACg                        >  1:2832152/1‑71 (MQ=255)                                        ggCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTa                                                      <  1:2332097/40‑1 (MQ=255)                                        ggCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTa                                                      <  1:3573669/40‑1 (MQ=255)                                           cTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACt                    <  1:3449006/71‑1 (MQ=255)                                                 ccGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCa                                <  1:343356/53‑1 (MQ=255)                                                 ccGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCa                                <  1:1503209/53‑1 (MQ=255)                                                 ccGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCa                                <  1:3514854/53‑1 (MQ=255)                                                 ccGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCa                                <  1:4936268/53‑1 (MQ=255)                                                 ccGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCa                                <  1:2031173/53‑1 (MQ=255)                                                 ccGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCa                                <  1:2655393/53‑1 (MQ=255)                                                 ccGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCa                                <  1:1735353/53‑1 (MQ=255)                                                   gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACGAAt                              >  1:240993/1‑53 (MQ=255)                                                   gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                              >  1:1217841/1‑53 (MQ=255)                                                   gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                              >  1:2270316/1‑53 (MQ=255)                                                   gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                              >  1:3406977/1‑53 (MQ=255)                                                   gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                              >  1:3103929/1‑53 (MQ=255)                                                   gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                              >  1:2322804/1‑53 (MQ=255)                                                   gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                              >  1:1558924/1‑53 (MQ=255)                                                   gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                              >  1:2825778/1‑53 (MQ=255)                                                   gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                              >  1:2768266/1‑53 (MQ=255)                                                   gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                              >  1:57622/1‑53 (MQ=255)                                                   gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                              >  1:2622858/1‑53 (MQ=255)                                                   gggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAAt                              >  1:806880/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:3829704/71‑1 (MQ=255)                                                    ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg           <  1:3781473/71‑1 (MQ=255)                                                    ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg           <  1:1895575/71‑1 (MQ=255)                                                    ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg           <  1:3679304/71‑1 (MQ=255)                                                    ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg           <  1:131728/71‑1 (MQ=255)                                                    ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg           <  1:4388070/71‑1 (MQ=255)                                                    ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg           <  1:1249020/71‑1 (MQ=255)                                                    ggAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTg           <  1:4500411/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)                                                      aaCGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGAc                                  >  1:3666372/1‑46 (MQ=255)                                                      aaCGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGaa         >  1:1648502/1‑71 (MQ=255)                                                      aaCGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGaa         >  1:1449930/1‑71 (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:118674/71‑3 (MQ=255)                                                           gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt    <  1:2730231/71‑3 (MQ=255)                                                           gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt    >  1:1057193/1‑69 (MQ=255)                                                           gTGGCGGTAGACCATGTCGTAGGGCAGGCCCAGCTCCAGGCCGATCTTGCGCACTTCGTCCTTGAACAgtt    >  1:3863788/1‑69 (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)                                                              |                                                                      AGTACTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCT  >  minE/1663961‑1664091

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