Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	269,535	C→T	59.7%	K12K (AAG→AAA)	ispA ←	geranyltranstransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	269,535	0	C	T	59.7%	32.3 / 73.1	62	K12K (AAG→AAA)	ispA	geranyltranstransferase
Reads supporting (aligned to +/- strand):  ref base C (12/13);  new base T (25/12);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

GTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAG  >  minE/269467‑269603                                                                     |                                                                     gTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCt                                                                     >  1:491100/1‑70 (MQ=255)    tCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAAc                                                                 >  1:3179616/1‑71 (MQ=255)    tCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAAc                                                                 >  1:1414388/1‑71 (MQ=255)    tCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAAc                                                                 >  1:519696/1‑71 (MQ=255)           cGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCt                                                          <  1:4054490/71‑1 (MQ=255)             ggAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAACg                                                                <  1:1467136/63‑1 (MQ=255)              gAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCTGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCg                                                       <  1:1139820/71‑1 (MQ=255)                   ttCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAACGCa                                                              <  1:748063/59‑1 (MQ=255)                   ttCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtg                                                  >  1:2062845/1‑71 (MQ=255)                    tCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTtgc                                                 >  1:3689478/1‑71 (MQ=255)                     cTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTtgct                                                >  1:1940855/1‑71 (MQ=255)                     cTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTtgct                                                >  1:2162681/1‑71 (MQ=255)                     cTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTtgct                                                >  1:327581/1‑71 (MQ=255)                     cTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCa                                                              >  1:1462888/1‑57 (MQ=255)                          aGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCgg                                            >  1:2812796/1‑70 (MQ=255)                             gCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAg                                        >  1:2816523/1‑71 (MQ=255)                               aGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTc                                      <  1:1134130/71‑1 (MQ=255)                                 tGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGt                                                    <  1:3734071/55‑1 (MQ=255)                                 tGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTcc                                     <  1:3214205/70‑1 (MQ=255)                                 tGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCGTCGAGt                                                    <  1:538767/55‑1 (MQ=255)                                 tGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACg                                                                <  1:397209/43‑1 (MQ=255)                                 tGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCa                                    >  1:1228071/1‑71 (MQ=255)                                      cGATAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACt                               <  1:27923/71‑1 (MQ=255)                                        aTAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCg                                                       <  1:3011934/45‑1 (MQ=255)                                         tAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCa                             >  1:449865/1‑70 (MQ=255)                                         tAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCa                             >  1:69317/1‑70 (MQ=255)                                         tAAAACGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCa                             >  1:2812140/1‑70 (MQ=255)                                          aaaaCGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTtggtgc                                              >  1:741493/1‑52 (MQ=255)                                          aaaaCGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTa                                 >  1:3087025/1‑65 (MQ=255)                                          aaaaCGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTa                                 >  1:3082417/1‑65 (MQ=255)                                          aaaaCGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTa                                 >  1:1626682/1‑65 (MQ=255)                                          aaaaCGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTa                                 >  1:220531/1‑65 (MQ=255)                                          aaaaCGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCa                             >  1:2047638/1‑69 (MQ=255)                                          aaaaCGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGa                                                      >  1:2464040/1‑44 (MQ=255)                                           aaaCGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCa                             <  1:460510/68‑1 (MQ=255)                                           aaaCGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTg                          >  1:787641/1‑71 (MQ=255)                                           aaaCGGCTCAGCGCCTGGTTGGCCTGCTTAACGCa                                                              <  1:2558091/35‑1 (MQ=255)                                           aaaCGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAg                                        <  1:882320/57‑1 (MQ=255)                                           aaaCGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTcc                                     <  1:3974142/60‑1 (MQ=255)                                            aaCGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCAtt                           >  1:3070040/1‑69 (MQ=255)                                             aCGGCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTc                        <  1:2309044/71‑1 (MQ=255)                                              cGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCa                                    >  1:3677051/1‑58 (MQ=255)                                               ggCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCAt                                   <  1:1848740/58‑1 (MQ=255)                                               ggCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCAt                                   <  1:1853611/58‑1 (MQ=255)                                                gCTCAGCGCCTGGTTGTCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCg                      >  1:892110/1‑70 (MQ=255)                                                gCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCgg                     >  1:1879803/1‑71 (MQ=255)                                                gCTCAGCGCCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCg                      >  1:3539659/1‑70 (MQ=255)                                                 cTCAGCGCCTGGTTGGCCTTCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCAtt                                  >  1:3565329/1‑57 (MQ=255)                                                    aGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGgtgt                  >  1:3667820/1‑70 (MQ=255)                                                     gcgcCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCg                      >  1:1620346/1‑65 (MQ=255)                                                       gcCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGaaa                                         >  1:1974084/1‑44 (MQ=255)                                                       gcCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTaa                >  1:3616094/1‑69 (MQ=255)                                                       gcCTGGTTGGCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTaa                >  1:3464066/1‑69 (MQ=255)                                                         cTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAgg            <  1:3977354/71‑1 (MQ=255)                                                         cTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAgg            <  1:3089944/71‑1 (MQ=255)                                                            gTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGTAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGt         <  1:2572286/71‑1 (MQ=255)                                                            gTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCa                                    >  1:3230701/1‑44 (MQ=255)                                                             ttGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGa                                           <  1:173820/36‑1 (MQ=255)                                                               ggCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGgt                    >  1:1920067/1‑57 (MQ=255)                                                                gCCTGTTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCa                             <  1:1189650/47‑1 (MQ=255)                                                                gCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTaga     <  1:2905783/71‑1 (MQ=255)                                                                 ccTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGgt                    <  1:3360037/55‑1 (MQ=255)                                                                     cTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTagagag  <  1:532441/69‑1 (MQ=255)                                                                     cTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTGagagag  <  1:537087/69‑1 (MQ=255)                                                                     |                                                                     GTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAG  >  minE/269467‑269603

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