Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,241,927	C→T	19.0%	V121V (GTC→GTT)	pykA →	pyruvate kinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,241,927	0	C	T	19.0%	146.8 / 20.8	94	V121V (GTC→GTT)	pykA	pyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base C (36/40);  new base T (10/8);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACT  >  minE/1241861‑1241980                                                                   |                                                      ccTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGAcgtcgt                                                     <  1:1337163/69‑1 (MQ=255)  ccgggTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTg                                                    <  1:1293258/67‑1 (MQ=255)  cTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTg                                                    <  1:871450/69‑1 (MQ=255)  cTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTg                                                    <  1:2654595/69‑1 (MQ=255)  cTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTg                                                    <  1:2650654/69‑1 (MQ=255)  cTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTg                                                    <  1:1948843/69‑1 (MQ=255)  cTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTg                                                    <  1:2309327/69‑1 (MQ=255)  cTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTg                                                    <  1:734589/69‑1 (MQ=255)       tAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGGCGTTGTGCCTg                                                >  1:1756794/1‑68 (MQ=255)       tAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTg                                                >  1:1642176/1‑68 (MQ=255)       tAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTg                                                >  1:620714/1‑68 (MQ=255)       tAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTg                                                >  1:1737727/1‑68 (MQ=255)       tAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTg                                                >  1:215646/1‑68 (MQ=255)       tAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTg                                                >  1:2466265/1‑68 (MQ=255)       tAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTg                                                >  1:2555302/1‑68 (MQ=255)       tAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTg                                                >  1:1111952/1‑68 (MQ=255)       tAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTg                                                >  1:3075312/1‑68 (MQ=255)       tAAAGGGGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTg                                                >  1:3476/1‑68 (MQ=255)            ttGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACa                                          >  1:1002990/2‑69 (MQ=255)            gTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACa                                          >  1:166974/1‑69 (MQ=255)            gTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACa                                          >  1:3329575/1‑69 (MQ=255)            gTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACa                                          >  1:3425478/1‑69 (MQ=255)            gTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACa                                          >  1:771578/1‑69 (MQ=255)            gTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACa                                          >  1:1384763/1‑69 (MQ=255)              gAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTGGTGACATc                                        <  1:2383152/69‑1 (MQ=255)              gAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTGGTGACATc                                        <  1:3135342/69‑1 (MQ=255)              gAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTGGTGACATc                                        <  1:3354882/69‑1 (MQ=255)              gAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTGGTGACATc                                        <  1:769552/69‑1 (MQ=255)              gAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTGGTGACATc                                        <  1:2492977/69‑1 (MQ=255)              gAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTGGTGACATc                                        <  1:2335110/69‑1 (MQ=255)              gAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTGGTGACATc                                        <  1:2090617/69‑1 (MQ=255)              gAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTTGTGCCTGGTGACATc                                        <  1:1884999/69‑1 (MQ=255)                aGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTg                                                    <  1:571853/55‑1 (MQ=255)                  gCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCctgc                                    >  1:2660492/1‑69 (MQ=255)                  gCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCctgc                                    >  1:239372/1‑69 (MQ=255)                       aaaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCctgctgctg                               <  1:1814998/69‑1 (MQ=255)                       aaaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCctgctgctg                               <  1:2692983/69‑1 (MQ=255)                       aaaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCctgctgctg                               <  1:1699618/69‑1 (MQ=255)                       aaaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCctgctgctg                               <  1:241971/69‑1 (MQ=255)                        aaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCctgctgctg                               <  1:1500175/68‑1 (MQ=255)                        aaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCctgctgctg                               <  1:2764033/68‑1 (MQ=255)                                         gACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCa              >  1:2122511/1‑68 (MQ=255)                                            tACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           >  1:605228/1‑68 (MQ=255)                                            tACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           >  1:1147233/1‑68 (MQ=255)                                            tACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           >  1:263088/1‑68 (MQ=255)                                            tACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           >  1:2713296/1‑68 (MQ=255)                                            tACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           >  1:2492138/1‑68 (MQ=255)                                            tACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           >  1:1769686/1‑68 (MQ=255)                                            tACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           >  1:679603/1‑68 (MQ=255)                                            tACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           >  1:2114525/1‑68 (MQ=255)                                            tACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           >  1:3232546/1‑68 (MQ=255)                                            tACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           >  1:1542815/1‑68 (MQ=255)                                                aaGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAg      <  1:2806932/69‑1 (MQ=255)                                                aaGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAg      <  1:1524247/69‑1 (MQ=255)                                                aaGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAg      <  1:3204044/69‑1 (MQ=255)                                                aaGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAg      <  1:1276837/69‑1 (MQ=255)                                                aaGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAg      <  1:2345760/69‑1 (MQ=255)                                                aaGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAg      <  1:1964180/69‑1 (MQ=255)                                                aaGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAg      <  1:527748/69‑1 (MQ=255)                                                 aGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCctgctgctg                               <  1:2031299/43‑1 (MQ=255)                                                 aGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAg      <  1:1727238/68‑1 (MQ=255)                                                 aGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAg      <  1:1568325/68‑1 (MQ=255)                                                 aGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGt     <  1:472642/69‑1 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTaaaa       >  1:2941057/1‑63 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:790998/1‑68 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:586840/1‑68 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:735332/1‑68 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:3145051/1‑68 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:1077364/1‑68 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:1131164/1‑68 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:132763/1‑68 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:1371921/1‑68 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:1296848/1‑68 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:1165870/1‑68 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:1770648/1‑68 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:18237/1‑68 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:2700600/1‑68 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:2825941/1‑68 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:2988659/1‑68 (MQ=255)                                                     ctgcctgcTGACGTCGTGCCTGCTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt  >  1:2546089/1‑68 (MQ=255)                                                       gcctgcTGACGTCGTGCCTGGTGACATCctgctgctg                               <  1:878923/37‑1 (MQ=255)                                                              gACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           <  1:2866080/50‑1 (MQ=255)                                                              gACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           <  1:1144160/50‑1 (MQ=255)                                                              gACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           <  1:898667/50‑1 (MQ=255)                                                              gACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           <  1:770144/50‑1 (MQ=255)                                                              gACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           <  1:907933/50‑1 (MQ=255)                                                              gACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           <  1:1650476/50‑1 (MQ=255)                                                              gACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           <  1:1723635/50‑1 (MQ=255)                                                              gACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           <  1:560324/50‑1 (MQ=255)                                                              gACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           <  1:1824098/50‑1 (MQ=255)                                                              gACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           <  1:3401521/50‑1 (MQ=255)                                                              gACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           <  1:2825074/50‑1 (MQ=255)                                                              gACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           <  1:2862087/50‑1 (MQ=255)                                                              gACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGCCCAGtt           <  1:124/50‑1 (MQ=255)                                                              gACGTCGAGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt           <  1:171976/50‑1 (MQ=255)                                                                   |                                                      CCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACT  >  minE/1241861‑1241980

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