Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,366,910	G→A	51.5%	H181H (CAC→CAT)	pepQ ←	proline dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,366,910	0	G	A	51.5%	‑1.8 / 88.9	97	H181H (CAC→CAT)	pepQ	proline dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (19/28);  new base A (27/23);  total (46/51)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCAGATAGGCAATGTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCAGTTTTGAAGG  >  minE/2366844‑2366974                                                                   |                                                                 ggTCAGATAGGCAATGTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTg                                                                <  1:1623736/69‑1 (MQ=255) ggTCAGATAGGCAATGTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTg                                                                <  1:3427145/69‑1 (MQ=255)     aGATAGGCAATGTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCg                                                            <  1:1765794/69‑1 (MQ=255)     aGATAGGCAATGTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCg                                                            <  1:2064063/69‑1 (MQ=255)     aGATAGGCAATGTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCg                                                            <  1:815420/69‑1 (MQ=255)       aTAGGCAATGTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGtt                                                          <  1:2869942/69‑1 (MQ=255)       aTAGGCAATGTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGtt                                                          <  1:130054/69‑1 (MQ=255)       aTAGGCAATGTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGtt                                                          <  1:2027909/69‑1 (MQ=255)       aTAGGCAATGTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGtt                                                          <  1:2419413/69‑1 (MQ=255)       aTAGGCAATGTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGGGACCGtt                                                          <  1:1631549/69‑1 (MQ=255)        tAGGCAATGTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGtt                                                          <  1:2563640/68‑1 (MQ=255)               tGTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCa                                                  <  1:672633/69‑1 (MQ=255)               tGTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCa                                                  <  1:448723/69‑1 (MQ=255)                gTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCa                                                  >  1:1179259/1‑68 (MQ=255)                gTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCa                                                  >  1:1524312/1‑68 (MQ=255)                gTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCAt                                                 <  1:2428815/69‑1 (MQ=255)                 ttGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCAtt                                                >  1:2930871/1‑69 (MQ=255)                 ttGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCAtt                                                >  1:2805443/1‑69 (MQ=255)                 ttGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCAtt                                                >  1:2285818/1‑69 (MQ=255)                 ttGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCAtt                                                >  1:291966/1‑69 (MQ=255)                 ttGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCAtt                                                >  1:2694739/1‑69 (MQ=255)                 ttGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCAtt                                                >  1:92114/1‑69 (MQ=255)                 ttGATATCG‑ACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCAttt                                               >  1:667608/1‑69 (MQ=255)                   gATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCAtttt                                              >  1:1519941/1‑69 (MQ=255)                         gAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTgcgc                                        >  1:722079/1‑69 (MQ=255)                             tCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCa                                    <  1:2532268/69‑1 (MQ=255)                             tCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCa                                    <  1:1690436/69‑1 (MQ=255)                             tCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCa                                    <  1:2213421/69‑1 (MQ=255)                             tCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCa                                    <  1:129868/69‑1 (MQ=255)                                      ccAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCt                                       >  1:2020809/1‑57 (MQ=255)                                      ccAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAg                           >  1:3087618/1‑69 (MQ=255)                                      ccAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAg                           >  1:1572289/1‑69 (MQ=255)                                      ccAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAg                           >  1:2188650/1‑69 (MQ=255)                                      ccAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAg                           >  1:3271252/1‑69 (MQ=255)                                      ccAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAg                           >  1:1135637/1‑69 (MQ=255)                                      ccAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAg                           >  1:3335491/1‑69 (MQ=255)                                      ccAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAg                           >  1:3399418/1‑69 (MQ=255)                                         gaacggaacgCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGcc                        <  1:2148899/69‑1 (MQ=255)                                         gaacggaacgCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGcc                        <  1:455731/69‑1 (MQ=255)                                         gaacggaacgCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGcc                        <  1:1133416/69‑1 (MQ=255)                                         gaacggaacgCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGcc                        <  1:156617/69‑1 (MQ=255)                                               aacgCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCg                  <  1:2713988/69‑1 (MQ=255)                                               aacgCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCg                  <  1:1524208/69‑1 (MQ=255)                                                 cgCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCa                                                  <  1:2562681/35‑1 (MQ=255)                                                     tcttctGCCGGGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:3101744/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:2887236/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:1112644/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:279989/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:956562/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:91848/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:2845335/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:809845/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:759655/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:3077651/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:2889802/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:473714/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:330012/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:3297549/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:2948199/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:1056488/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:1923652/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:2709107/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:198673/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:1664552/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:2132236/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:1496493/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:2286849/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:2364270/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:2455845/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:248721/1‑68 (MQ=255)                                                     tcttctGCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTcg            >  1:1305082/1‑68 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:755124/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:1310760/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:723905/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:1328287/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:1490656/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:1285275/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:887139/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:1171457/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:1696124/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:1937165/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:1851070/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:2044591/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:2180841/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:317471/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:3108996/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:2309961/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:2965090/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:2403718/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:2512767/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:2552048/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:2650775/68‑1 (MQ=255)                                                           gCCGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:2700930/68‑1 (MQ=255)                                                            ccGCGCGATGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCGGtttt       <  1:2587741/67‑1 (MQ=255)                                                               cgcgGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCAGTTTTGAAgg  <  1:774443/69‑1 (MQ=255)                                                               cgcgGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCAGTTTTGAAgg  <  1:1140852/69‑1 (MQ=255)                                                                gcgGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCAGTTTTGAAgg  <  1:3412115/68‑1 (MQ=255)                                                                   |                                                                 GGTCAGATAGGCAATGTTGATATCGAACTCGCTCATGCCAGAACGGAACGCTTCTTCTGCCGCGCGGTGACCGTTGACCGCCATTTTCTGCGCTTCACGCATACAGGCCAGTTCGTACTCAGTTTTGAAGG  >  minE/2366844‑2366974

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