Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,601,306	T→C	26.1%	F319F (TTT→TTC)	prlC →	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,601,306	0	T	C	26.1%	90.5 / 29.4	80	F319F (TTT→TTC)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base T (33/26);  new base C (13/8);  total (46/34)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCGCGCGTCCACAAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAA  >  minE/2601257‑2601373                                                  |                                                                    aaCGCGCGCGTCCACAAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGcc                                                            >  1:2930834/1‑59 (MQ=255) aaCGCGCGCGTCCACAAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGc                                                   >  1:1741162/1‑68 (MQ=255) aaCGCGCGCGTCCACAAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCg                                                  >  1:3038803/1‑69 (MQ=255) aaCGCGCGCGTCCACAAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCg                                                  >  1:1923010/1‑69 (MQ=255) aaCGCGCGCGTCCACAAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCg                                                  >  1:1737348/1‑69 (MQ=255) aaCGCGCGCGTCCACAAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCg                                                  >  1:1495333/1‑69 (MQ=255) aaCGCGCGCGTCCACAAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCg                                                  >  1:1576752/1‑69 (MQ=255)             cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATg                                            <  1:2532667/63‑1 (MQ=255)             cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATg                                            <  1:1902116/63‑1 (MQ=255)             cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATg                                            <  1:1860600/63‑1 (MQ=255)             cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATg                                            <  1:1936675/63‑1 (MQ=255)             cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATg                                            <  1:1594426/63‑1 (MQ=255)             cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATg                                            <  1:1498728/63‑1 (MQ=255)             cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATg                                            <  1:2271193/63‑1 (MQ=255)             cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATg                                            <  1:1465289/63‑1 (MQ=255)             cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATg                                            <  1:2290197/63‑1 (MQ=255)             cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATg                                            <  1:3380735/63‑1 (MQ=255)             cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATg                                            <  1:2424474/63‑1 (MQ=255)             cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATg                                            <  1:1231073/63‑1 (MQ=255)             cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATg                                            <  1:781293/63‑1 (MQ=255)             cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGc                                      >  1:1653543/1‑69 (MQ=255)             cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGc                                      >  1:454502/1‑69 (MQ=255)              acaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATg                                            <  1:1568385/62‑1 (MQ=255)               caAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCa                                     <  1:2661254/68‑1 (MQ=255)                        aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAAt                                                        >  1:3063654/1‑40 (MQ=25)                        aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                           >  1:1662047/1‑69 (MQ=255)                        aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                           >  1:274868/1‑69 (MQ=255)                        aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                           >  1:2035587/1‑69 (MQ=255)                        aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                           >  1:3214155/1‑69 (MQ=255)                        aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                           >  1:3300673/1‑69 (MQ=255)                        aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                           >  1:232975/1‑69 (MQ=255)                        aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                           >  1:2672402/1‑69 (MQ=255)                        aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                           >  1:250201/1‑69 (MQ=255)                        aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                           >  1:838419/1‑69 (MQ=255)                        aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGa                            >  1:2821841/1‑68 (MQ=255)                        aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGa                            >  1:2216230/1‑68 (MQ=255)                        aaaGAGCTGGCGCAACTGCGCGCATTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTTGGa                            >  1:3059478/1‑68 (MQ=255)                              cTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCt                     <  1:489821/69‑1 (MQ=255)                              cTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCt                     <  1:932074/69‑1 (MQ=255)                               tGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCt                     <  1:1681689/68‑1 (MQ=255)                               tGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCt                     <  1:1462500/68‑1 (MQ=255)                               tGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCt                     <  1:2809880/68‑1 (MQ=255)                                 gcgcAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCt                     <  1:1391701/66‑1 (MQ=255)                                   gcAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTtac                  <  1:2881638/67‑1 (MQ=255)                                    cAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTtac                  <  1:1491606/66‑1 (MQ=255)                                      aTTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTtagtacg              >  1:2726749/1‑67 (MQ=255)                                      aTTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACa              >  1:571097/1‑68 (MQ=255)                                      aTTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACa              >  1:2604637/1‑68 (MQ=255)                                      aTTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACa              >  1:2934895/1‑68 (MQ=255)                                      aTTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACa              >  1:292138/1‑68 (MQ=255)                                      aTTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACa              >  1:1996788/1‑68 (MQ=255)                                        tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAg                                          >  1:3092107/1‑38 (MQ=255)                                        tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGt                                >  1:1058318/1‑48 (MQ=255)                                        tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                              >  1:3135559/1‑50 (MQ=255)                                        tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                              >  1:596238/1‑50 (MQ=255)                                        tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                              >  1:1366911/1‑50 (MQ=255)                                        tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                              >  1:630278/1‑50 (MQ=255)                                        tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                              >  1:1371854/1‑50 (MQ=255)                                        tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                              >  1:633967/1‑50 (MQ=255)                                        tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                              >  1:704475/1‑50 (MQ=255)                                        tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                              >  1:860827/1‑50 (MQ=255)                                        tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTg                               >  1:2949637/1‑49 (MQ=255)                                        tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTg                               >  1:2893971/1‑49 (MQ=255)                                         gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa          <  1:324720/69‑1 (MQ=255)                                         gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa          <  1:2533173/69‑1 (MQ=255)                                         gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa          <  1:1125375/69‑1 (MQ=255)                                         gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa          <  1:775103/69‑1 (MQ=255)                                         gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa          <  1:727847/69‑1 (MQ=255)                                         gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa          <  1:198917/69‑1 (MQ=255)                                         gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa          <  1:2859678/69‑1 (MQ=255)                                         gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa          <  1:618874/69‑1 (MQ=255)                                         gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAg                                          >  1:266079/1‑37 (MQ=255)                                         gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGt                                <  1:2393132/47‑1 (MQ=255)                                         gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGt                                <  1:1878936/47‑1 (MQ=255)                                         gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                              <  1:370390/49‑1 (MQ=255)                                         gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCg           >  1:606346/1‑68 (MQ=255)                                         gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCg           >  1:2710773/1‑68 (MQ=255)                                         gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCg           >  1:178467/1‑68 (MQ=255)                                         gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCg           >  1:3287138/1‑68 (MQ=255)                                         gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCg           >  1:3125840/1‑68 (MQ=255)                                                 ttGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGaaaaacaaa  <  1:1986448/69‑1 (MQ=255)                                                 ttGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGaaaaacaaa  <  1:1219168/69‑1 (MQ=255)                                                  |                                                                    AACGCGCGCGTCCACAAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAA  >  minE/2601257‑2601373

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