Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,601,336	G→A	39.9%	L329L (TTG→TTA)	prlC →	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,601,336	0	G	A	39.9%	40.2 / 56.8	85	L329L (TTG→TTA)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base G (29/22);  new base A (12/22);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 7.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

CACAAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATC  >  minE/2601269‑2601393                                                                    |                                                          cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGc                                                          >  1:454502/1‑69 (MQ=255) cacaAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGc                                                          >  1:1653543/1‑69 (MQ=255)   caAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCa                                                         <  1:2661254/68‑1 (MQ=255)            aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:2035587/1‑69 (MQ=255)            aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:838419/1‑69 (MQ=255)            aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:3300673/1‑69 (MQ=255)            aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:3214155/1‑69 (MQ=255)            aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:232975/1‑69 (MQ=255)            aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:250201/1‑69 (MQ=255)            aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:274868/1‑69 (MQ=255)            aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:2672402/1‑69 (MQ=255)            aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:1662047/1‑69 (MQ=255)            aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGa                                                >  1:2216230/1‑68 (MQ=255)            aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGa                                                >  1:2821841/1‑68 (MQ=255)            aaaGAGCTGGCGCAACTGCGCGCATTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTTGGa                                                >  1:3059478/1‑68 (MQ=255)                  cTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCt                                         <  1:932074/69‑1 (MQ=255)                  cTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCt                                         <  1:489821/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:1996788/1‑68 (MQ=255)                          aTTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACa                                  >  1:292138/1‑68 (MQ=255)                          aTTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACa                                  >  1:2934895/1‑68 (MQ=255)                          aTTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACa                                  >  1:571097/1‑68 (MQ=255)                          aTTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACa                                  >  1:2604637/1‑68 (MQ=255)                            tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGt                                                    >  1:1058318/1‑48 (MQ=255)                            tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:633967/1‑50 (MQ=255)                            tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:1366911/1‑50 (MQ=255)                            tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:1371854/1‑50 (MQ=255)                            tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:860827/1‑50 (MQ=255)                            tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:3135559/1‑50 (MQ=255)                            tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:596238/1‑50 (MQ=255)                            tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:630278/1‑50 (MQ=255)                            tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:704475/1‑50 (MQ=255)                            tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTg                                                   >  1:2949637/1‑49 (MQ=255)                            tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTg                                                   >  1:2893971/1‑49 (MQ=255)                             gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                              <  1:775103/69‑1 (MQ=255)                             gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                              <  1:324720/69‑1 (MQ=255)                             gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                              <  1:727847/69‑1 (MQ=255)                             gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                              <  1:1125375/69‑1 (MQ=255)                             gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                              <  1:618874/69‑1 (MQ=255)                             gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                              <  1:2859678/69‑1 (MQ=255)                             gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                              <  1:198917/69‑1 (MQ=255)                             gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                              <  1:2533173/69‑1 (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:2710773/1‑68 (MQ=255)                             gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCg                               >  1:3287138/1‑68 (MQ=255)                             gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCg                               >  1:178467/1‑68 (MQ=255)                             gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCg                               >  1:606346/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)                                       gCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGata                                              >  1:676326/1‑43 (MQ=255)                                       gCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGaaaaacaaa                      >  1:270951/1‑67 (MQ=255)                                       gCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGaaaaacaaa                      >  1:1205804/1‑67 (MQ=255)                                       gCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGaaaaacaaa                      >  1:3420056/1‑67 (MQ=255)                                       gCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGaaaaacaa                       >  1:1289367/1‑66 (MQ=255)                                          aaaGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAg                 <  1:2535038/69‑1 (MQ=255)                                          aaaGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAg                 <  1:2330993/69‑1 (MQ=255)                                                         gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc  <  1:656240/69‑1 (MQ=255)                                                         gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc  <  1:418675/69‑1 (MQ=255)                                                         gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc  <  1:3356145/69‑1 (MQ=255)                                                         gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc  <  1:2650800/69‑1 (MQ=255)                                                         gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc  <  1:2926715/69‑1 (MQ=255)                                                         gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc  <  1:2830903/69‑1 (MQ=255)                                                         gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:728344/68‑1 (MQ=255)                                                         gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:1131083/68‑1 (MQ=255)                                                         gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:857767/68‑1 (MQ=255)                                                         gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:2011993/68‑1 (MQ=255)                                                         gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:239376/68‑1 (MQ=255)                                                         gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:1965755/68‑1 (MQ=255)                                                         gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:610645/68‑1 (MQ=255)                                                         gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:264382/68‑1 (MQ=255)                                                         gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:1165911/68‑1 (MQ=255)                                                         gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:3098895/68‑1 (MQ=255)                                                         gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:2720271/68‑1 (MQ=255)                                                         gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:1342423/68‑1 (MQ=255)                                                         gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:1366775/68‑1 (MQ=255)                                                         gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTACACAGCAt   <  1:1137733/68‑1 (MQ=255)                                                                    |                                                          CACAAGGCGAAAAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATC  >  minE/2601269‑2601393

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