Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,045,790	C→T	20.0%	L99L (TTG→TTA)	ydhA ←	predicted lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,045,790	0	C	T	20.0%	80.5 / 15.4	40	L99L (TTG→TTA)	ydhA	predicted lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (16/16);  new base T (5/3);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

TGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACAT  >  minE/1045721‑1045848                                                                      |                                                           tGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTc                                                            >  1:2322911/1‑70 (MQ=255)   ggTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTTaa                                                          <  1:414072/70‑1 (MQ=255)     tGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGAc                                                       >  1:3082827/1‑71 (MQ=255)        cGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTTAAGACGAt                                                    >  1:1977539/1‑71 (MQ=255)        cGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGAt                                                    >  1:693571/1‑71 (MQ=255)         gCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATg                                                   <  1:3043182/71‑1 (MQ=255)            aTTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGAc                                                       <  1:1678727/64‑1 (MQ=255)             ttGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCgg                                                <  1:50182/70‑1 (MQ=255)             ttGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCgg                                                <  1:2254520/70‑1 (MQ=255)              tgtgCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTTAAGACGCTGCgg                                                >  1:1756501/1‑69 (MQ=255)                   gccgccCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGAc                                                       <  1:1855729/57‑1 (MQ=255)                      gccCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTtat                                      <  1:236674/71‑1 (MQ=255)                        cccTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTcgcg                                           <  1:1061831/64‑1 (MQ=255)                        cccTGGAAAAATCTCAACGCTGAGGATTTTGTAACTGACAGTTATTCAAGAc                                                       <  1:2865088/52‑1 (MQ=255)                         ccTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTtataca                                   <  1:3320226/71‑3 (MQ=255)                          cTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGa                                                     >  1:56187/1‑52 (MQ=255)                           tGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTTAAGACGATGCGGTCGCGTTtat                                      <  1:3537114/66‑1 (MQ=255)                            ggAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTTAAGACGATGCGGTCGCGTTTATAGACAg                                >  1:2632498/1‑71 (MQ=255)                            ggAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAg                                >  1:1971908/1‑71 (MQ=255)                            ggAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAg                                >  1:2629737/1‑71 (MQ=255)                              aaaaaTCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAg                                <  1:1769959/69‑1 (MQ=255)                              aaaaaTCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTc                              <  1:7245/71‑1 (MQ=255)                                aaaTCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGtt                                         >  1:409921/1‑58 (MQ=255)                                aaaTCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGAc                                  >  1:414364/1‑65 (MQ=255)                                 aaTCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTta                                       >  1:931448/1‑59 (MQ=255)                                     tcAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCa                        >  1:2626340/1‑70 (MQ=255)                                       aaCGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGa                                   >  1:3948895/1‑57 (MQ=255)                                       aaCGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATc                      <  1:3216167/70‑1 (MQ=255)                                            tgtgGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCtt                 >  1:1467933/1‑70 (MQ=255)                                                  ttttGTAACTGACAGTTATTTAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCa                        <  1:1975451/57‑1 (MQ=255)                                                    ttGTAACTGACAGTTATTCAAGACGATGCGGTCGCGt                                          >  1:921090/1‑37 (MQ=255)                                                     tGTAACTGACAGTTATTTAAGACGATGCGGTCGCGTTTATAGACAGTCGCt                           >  1:1899011/1‑51 (MQ=255)                                                     tGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGa        >  1:1588499/1‑70 (MQ=255)                                                      gTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGaa       <  1:2233562/70‑1 (MQ=255)                                                       tAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGaaa      <  1:1297665/70‑1 (MQ=255)                                                          cTGACAGTTATTTAAGACGATGCGGTCGCGTTTATAGACAGTCGCt                           >  1:329065/1‑46 (MQ=255)                                                          cTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCt                           >  1:3273336/1‑46 (MQ=255)                                                          cTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCa                        <  1:730633/49‑1 (MQ=255)                                                          cTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGaa       >  1:578405/1‑66 (MQ=255)                                                          cTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACAt  <  1:3626540/71‑1 (MQ=255)                                                                   aTTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCtt                          >  1:2068303/1‑38 (MQ=255)                                                                      |                                                           TGGGTGACGCTATTGTGCGCCGCCCCTGGAAAAATCTCAACGCTGTGGATTTTGTAACTGACAGTTATTCAAGACGATGCGGTCGCGTTTATAGACAGTCGCTTCATCGCCTTTCGACCAGAAAACAT  >  minE/1045721‑1045848

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