Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,241,827	T→C	5.9%	L88P (CTC→CCC)	pykA →	pyruvate kinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,241,827	0	T	C	5.9%	202.4 / 3.4	68	L88P (CTC→CCC)	pykA	pyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base T (39/25);  new base C (2/2);  total (41/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

GGCTATTCTGGGTGACCTCCAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAG  >  minE/1241759‑1241892                                                                     |                                                                  ggCTATTCTGGGTGACCTCCAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTc                                                                  <  1:1010332/70‑1 (MQ=255)      ttCTGGGTGACCTCCAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATAtt                                                            <  1:1884159/71‑1 (MQ=255)       tCTGGGTGACCTCCAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTg                                                           >  1:564208/1‑71 (MQ=255)         tGGGTGACCTCCAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTggg                                                         >  1:3179062/1‑71 (MQ=255)             tGACCTCCAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGa                                                       >  1:425488/1‑69 (MQ=255)             tGACCTCCAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGAt                                                      >  1:726904/1‑70 (MQ=255)               aCCTCCAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATaa                                                    >  1:3378332/1‑70 (MQ=255)                   ccAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTcc                                               >  1:2020053/1‑71 (MQ=255)                   ccAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTcc                                               >  1:2323928/1‑71 (MQ=255)                    cAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCt                                              >  1:3642796/1‑71 (MQ=255)                     aGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTg                                             >  1:1586023/1‑71 (MQ=255)                      ggggCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTgg                                                          >  1:1897942/1‑57 (MQ=255)                         gCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCt                                           >  1:3456474/1‑69 (MQ=255)                         gCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTc                                          >  1:121321/1‑70 (MQ=255)                            cAAAATCCGTGTATCCACCTTTTAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACg                                      <  1:3375660/71‑1 (MQ=255)                            cAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTc                                          <  1:1879804/67‑1 (MQ=255)                            cAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACg                                      <  1:3502816/71‑1 (MQ=255)                            cAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACg                                      >  1:2761224/1‑71 (MQ=255)                              aaaTCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGa                                        >  1:2794393/1‑67 (MQ=255)                              aaaTCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGc                                     >  1:1924452/1‑70 (MQ=255)                               aaTCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCa                                   <  1:1476661/71‑1 (MQ=255)                                aTCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCt                                                                   >  1:1475845/1‑38 (MQ=255)                                aTCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCt                                           >  1:1514238/1‑62 (MQ=255)                                      gtATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTggg                            >  1:1097911/1‑71 (MQ=255)                                      gtATCCACCTTTAAAGAAGGCAAAGTTTTCCCCAATATTGAGGaaa                                                     >  1:128858/1‑44 (MQ=37)                                         tCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTggg                                                         <  1:2368285/39‑1 (MQ=255)                                         tCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATaaa                                                   >  1:862289/1‑45 (MQ=255)                                         tCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATaaa                                                   >  1:64619/1‑45 (MQ=255)                                           cACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTc                                                >  1:425584/1‑46 (MQ=255)                                           cACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAg                       <  1:1760775/71‑1 (MQ=255)                                           cACCTTTAAAGAAGGCAAAGTTTTCCCCAATATTgagg                                                        <  1:2548414/38‑1 (MQ=37)                                           cACCTTTAAAGAAGGCAAAGTTTTCCCCAATATTgagg                                                        <  1:596877/38‑1 (MQ=37)                                            aCCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAgg                      >  1:983048/1‑71 (MQ=255)                                             ccTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTc                                                >  1:1549800/1‑44 (MQ=255)                                             ccTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTc                                                >  1:1722789/1‑44 (MQ=255)                                             ccTTTAAAGAAGGCAAAGTTTTCCCCAATATTGAGGaa                                                      >  1:1788757/1‑37 (MQ=25)                                              cTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGt                           >  1:2398537/1‑64 (MQ=255)                                              cTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTg                    <  1:3542862/71‑1 (MQ=255)                                              cTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTg                    <  1:3524818/71‑1 (MQ=255)                                              cTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTg                    <  1:3029474/71‑1 (MQ=255)                                              cTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTg                    <  1:1783126/71‑1 (MQ=255)                                                ttAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTg                    >  1:2526261/1‑69 (MQ=255)                                                 tAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAg                 >  1:2128253/1‑71 (MQ=255)                                                 tAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTAGaa                                       >  1:1376988/1‑48 (MQ=255)                                                   aagaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAgg                <  1:1684062/70‑1 (MQ=255)                                                   aagaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGc               <  1:839859/71‑1 (MQ=255)                                                     gaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGAc                                       >  1:194780/1‑45 (MQ=255)                                                       agGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCt                               <  1:2621277/51‑1 (MQ=255)                                                       agGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCt                               <  1:1403862/51‑1 (MQ=255)                                                       agGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCt                               <  1:1706894/51‑1 (MQ=255)                                                       agGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCt                               <  1:3205049/51‑1 (MQ=255)                                                       agGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCt                               <  1:1911749/51‑1 (MQ=255)                                                         gCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTgg                             >  1:1607992/1‑51 (MQ=255)                                                         gCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACaaa         <  1:2629915/71‑1 (MQ=255)                                                         gCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACaaa         <  1:1544399/71‑1 (MQ=255)                                                          cAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAg        <  1:64863/71‑1 (MQ=255)                                                          cAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAg        <  1:1695617/71‑1 (MQ=255)                                                           aaaGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAg        >  1:2514120/1‑70 (MQ=255)                                                             aGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAAc                                 >  1:449006/1‑43 (MQ=255)                                                             aGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAAc                                 >  1:286472/1‑43 (MQ=255)                                                             aGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGaaa     <  1:1629803/71‑1 (MQ=255)                                                             aGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGaa      >  1:1292498/1‑70 (MQ=255)                                                             aGTTTTCCTCAATATTGGGGATAAATTACTGCTCGACGCCAACCTg                              >  1:143794/1‑46 (MQ=255)                                                               ttttCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGaaa     >  1:2451382/1‑69 (MQ=255)                                                               ttttCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGaaa     >  1:1036136/1‑69 (MQ=255)                                                                tttCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAg  >  1:3023068/1‑71 (MQ=255)                                                                tttCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAg  >  1:2518863/1‑71 (MQ=255)                                                                tttCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAg  >  1:1272958/1‑71 (MQ=255)                                                                 ttCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAg  <  1:1289818/70‑1 (MQ=255)                                                                     |                                                                  GGCTATTCTGGGTGACCTCCAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAG  >  minE/1241759‑1241892

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