Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,601,339	G→A	41.4%	Q330Q (CAG→CAA)	prlC →	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,601,339	0	G	A	41.4%	32.5 / 57.0	82	Q330Q (CAG→CAA)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base G (27/21);  new base A (12/22);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 7.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATC  >  minE/2601280‑2601393                                                            |                                                       aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:838419/1‑69 (MQ=255) aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:3300673/1‑69 (MQ=255) aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:3214155/1‑69 (MQ=255) aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:274868/1‑69 (MQ=255) aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:2672402/1‑69 (MQ=255) aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:250201/1‑69 (MQ=255) aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:232975/1‑69 (MQ=255) aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:2035587/1‑69 (MQ=255) aaaGAGCTGGCGCAACTGCGCGCCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGat                                               >  1:1662047/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:932074/69‑1 (MQ=255)       cTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCt                                         <  1:489821/69‑1 (MQ=255)        tGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCt                                         <  1:1462500/68‑1 (MQ=255)        tGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCt                                         <  1:1681689/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:292138/1‑68 (MQ=255)               aTTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACa                                  >  1:571097/1‑68 (MQ=255)               aTTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACa                                  >  1:2934895/1‑68 (MQ=255)               aTTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACa                                  >  1:2604637/1‑68 (MQ=255)               aTTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACa                                  >  1:1996788/1‑68 (MQ=255)                 tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGt                                                    >  1:1058318/1‑48 (MQ=255)                 tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:1366911/1‑50 (MQ=255)                 tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:860827/1‑50 (MQ=255)                 tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:704475/1‑50 (MQ=255)                 tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:633967/1‑50 (MQ=255)                 tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:630278/1‑50 (MQ=255)                 tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:596238/1‑50 (MQ=255)                 tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:3135559/1‑50 (MQ=255)                 tGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTgg                                                  >  1:1371854/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:727847/69‑1 (MQ=255)                  gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                              <  1:618874/69‑1 (MQ=255)                  gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                              <  1:198917/69‑1 (MQ=255)                  gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                              <  1:2533173/69‑1 (MQ=255)                  gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                              <  1:2859678/69‑1 (MQ=255)                  gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                              <  1:1125375/69‑1 (MQ=255)                  gcgcgcCTTCGCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                              <  1:324720/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:606346/1‑68 (MQ=255)                  gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCg                               >  1:178467/1‑68 (MQ=255)                  gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCg                               >  1:3125840/1‑68 (MQ=255)                  gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCg                               >  1:3287138/1‑68 (MQ=255)                  gCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCg                               >  1:2710773/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:1205804/1‑67 (MQ=255)                            gCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGaaaaacaaa                      >  1:270951/1‑67 (MQ=255)                            gCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGaaaaacaaa                      >  1:3420056/1‑67 (MQ=255)                            gCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGaaaaacaa                       >  1:1289367/1‑66 (MQ=255)                               aaaGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAg                 <  1:2330993/69‑1 (MQ=255)                               aaaGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAg                 <  1:2535038/69‑1 (MQ=255)                                              gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc  <  1:2650800/69‑1 (MQ=255)                                              gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc  <  1:656240/69‑1 (MQ=255)                                              gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc  <  1:2830903/69‑1 (MQ=255)                                              gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc  <  1:418675/69‑1 (MQ=255)                                              gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc  <  1:3356145/69‑1 (MQ=255)                                              gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc  <  1:2926715/69‑1 (MQ=255)                                              gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:3098895/68‑1 (MQ=255)                                              gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:1131083/68‑1 (MQ=255)                                              gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:857767/68‑1 (MQ=255)                                              gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:1342423/68‑1 (MQ=255)                                              gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:1165911/68‑1 (MQ=255)                                              gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:728344/68‑1 (MQ=255)                                              gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:1965755/68‑1 (MQ=255)                                              gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:2720271/68‑1 (MQ=255)                                              gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:264382/68‑1 (MQ=255)                                              gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:239376/68‑1 (MQ=255)                                              gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:1366775/68‑1 (MQ=255)                                              gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:2011993/68‑1 (MQ=255)                                              gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt   <  1:610645/68‑1 (MQ=255)                                              gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTACACAGCAt   <  1:1137733/68‑1 (MQ=255)                                                            |                                                       AAAGAGCTGGCGCAATTGCGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATC  >  minE/2601280‑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: ‑