Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,492,284	C→T	12.9%	L361L (CTG→TTG)	gyrB →	DNA gyrase, subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,492,284	0	C	T	12.9%	139.0 / 10.1	54	L361L (CTG→TTG)	gyrB	DNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base C (23/24);  new base T (2/5);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACAAACTGGTTTCTTCTGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCTGCCCGTGCCCGTGAAGC  >  minE/2492214‑2492354                                                                       |                                                                       aGACAAACTGGTTTCTTCTGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATAcc                                                                        <  1:3557769/71‑1 (MQ=255)   aCAAACTGGTTTCTTCTGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACttg                                                                      <  1:3187928/71‑1 (MQ=255)   aCAAACTGGTTTCTTCTGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACctg                                                                      >  1:1745378/1‑71 (MQ=255)      aaCTGGTTTCTTCTGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACctgctg                                                                   <  1:2419147/71‑1 (MQ=255)           gTTTCTTCTGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGaaa                                                               <  1:1672666/70‑1 (MQ=255)           gTTTCTTCTGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGaaa                                                               >  1:3096646/1‑70 (MQ=255)           gTTTCTTCTGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGaaa                                                               <  1:2477936/70‑1 (MQ=255)           gTTTCTTCTGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGaaa                                                               <  1:1402405/70‑1 (MQ=255)            tttCTTCTGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACTTGCTGGaaa                                                               <  1:2833873/69‑1 (MQ=255)            tttCTTCTGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGaa                                                                <  1:3952775/68‑1 (MQ=255)              tcttctGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGa                                                                 <  1:2204872/65‑1 (MQ=255)              tcttctGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAAcc                                                            >  1:3186111/1‑70 (MQ=255)                ttctGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACttg                                                                      >  1:1486559/1‑58 (MQ=255)                ttctGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACTTGCTGGaaa                                                               >  1:2398808/1‑65 (MQ=255)                ttctGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCa                                                          >  1:797846/1‑70 (MQ=255)                  ctGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAAcc                                                       >  1:3236190/1‑71 (MQ=255)                    gAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGa                                                     >  1:2457879/1‑71 (MQ=255)                          aaaTCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGa                                                >  1:2541923/1‑70 (MQ=255)                            aTCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGAcgc                                                  >  1:2087815/1‑66 (MQ=255)                               ggcggTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATc                                           <  1:1323279/70‑1 (MQ=255)                               ggcggTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATc                                           <  1:1642020/70‑1 (MQ=255)                                gcggTTGAACAGCAGATGAACGAACTGCTGGCAGAATACctgct                                                                    >  1:1504975/1‑44 (MQ=255)                                 cggTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAAcc                                                            <  1:611982/51‑1 (MQ=255)                                 cggTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTg                                        >  1:1145362/1‑71 (MQ=255)                                   gTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTgtt                                      <  1:1857374/71‑3 (MQ=255)                                   gTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGt                                      <  1:1530659/71‑1 (MQ=255)                                   gTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGt                                      <  1:2896565/71‑1 (MQ=255)                                    ttGAACAGCAGATGAACGAACTGCTGGCAGAATACTTGCTGGaa                                                                <  1:2922899/44‑1 (MQ=255)                                    ttGAACAGCAGATGAACGAACTGCTGGCAGAATACTTGCTGGAAAACCCAAcc                                                       <  1:3485222/53‑1 (MQ=255)                                    ttGAACAGCAGATGAACGAACTGCTGGCAGAATACTTGCTGGAAAACCCAACCGACGCGAAAATCGTGGtt                                     <  1:2431837/71‑1 (MQ=255)                                     tGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTg                                    >  1:1830541/1‑71 (MQ=255)                                     tGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTg                                    >  1:63978/1‑71 (MQ=255)                                      gAACAGCAGATGAACGAACTGCTGGCAGAATACctgc                                                                     <  1:1803688/37‑1 (MQ=255)                                       aaCAGAAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATc                                           >  1:1746166/1‑62 (MQ=255)                                        aCAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCa                                 <  1:3753644/71‑1 (MQ=255)                                              gATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTa                           >  1:1765349/1‑71 (MQ=255)                                                tGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTg                                        >  1:2006349/1‑56 (MQ=255)                                                 gaacgaacTGCTGGCAGAATACTTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATc                         >  1:808862/1‑70 (MQ=255)                                                 gaacgaacTGCTGGCAGAATACCTGCTGGAAAAccc                                                           >  1:3867955/1‑36 (MQ=255)                                                   acgaacTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGAt                      >  1:556985/1‑71 (MQ=255)                                                     gaacTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGAtgc                    >  1:1390489/1‑71 (MQ=255)                                                     gaacTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGAtgc                    >  1:2875381/1‑71 (MQ=255)                                                      aacTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATc                                           <  1:3922532/47‑1 (MQ=255)                                                      aacTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGAtgc                    >  1:626433/1‑70 (MQ=255)                                                         tgctgGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTg                                        >  1:1240832/1‑47 (MQ=255)                                                         tgctgGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGAtgct                   <  1:1852235/68‑1 (MQ=255)                                                          gctgGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAAtt                            <  1:1570278/58‑1 (MQ=255)                                                           ctgGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGtt                                     >  1:3243750/1‑48 (MQ=255)                                                           ctgGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCaaaa                              <  1:2928235/55‑1 (MQ=255)                                                              gCAGAATACTTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCtgcccgtgc           <  1:1190193/71‑1 (MQ=255)                                                                aGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCtgcccgtgccc         >  1:1786109/1‑71 (MQ=255)                                                                 gAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCtgcccgtgcccg        <  1:2351638/71‑1 (MQ=255)                                                                 gAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCtgcccgtgcccg        <  1:2876249/71‑1 (MQ=255)                                                                  aaTACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCTgcccgtgcccgt       >  1:3118229/1‑71 (MQ=255)                                                                    tACCTGCTTGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCTGCCCGTGCCCGTGa     <  1:3828870/71‑1 (MQ=255)                                                                    tACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCTGcccgtgcccgtg      >  1:2316501/1‑70 (MQ=255)                                                                    tACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCTGCCCGTGCCCGTGa     <  1:1964731/71‑1 (MQ=255)                                                                    tACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCTGCCCGTGCCCGTGa     <  1:2012192/71‑1 (MQ=255)                                                                    tACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCTGCCCGTGCCCGTGa     <  1:128667/71‑1 (MQ=255)                                                                    tACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCTGCCCGTGCCCGTGa     <  1:680698/71‑1 (MQ=255)                                                                    tACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCTGCCCGTGCCCGTGa     <  1:2807360/71‑1 (MQ=255)                                                                      ccTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCTGCCCGTGCCCGTGAAg   <  1:1973655/71‑1 (MQ=255)                                                                       ctgctgGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCTGCCCGTGCCCGTGAAgc  <  1:176732/71‑1 (MQ=255)                                                                       |                                                                       AGACAAACTGGTTTCTTCTGAGGTGAAATCGGCGGTTGAACAGCAGATGAACGAACTGCTGGCAGAATACCTGCTGGAAAACCCAACCGACGCGAAAATCGTGGTTGGCAAAATTATCGATGCTGCCCGTGCCCGTGAAGC  >  minE/2492214‑2492354

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