Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	455,013	G→A	29.2%	S24S (AGC→AGT)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	455,013	0	G	A	29.2%	130.2 / 64.8	89	S24S (AGC→AGT)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (27/36);  new base A (9/17);  total (36/53)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.27e-01

GTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATT  >  minE/454947‑455081                                                                   |                                                                     gTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa                                                                  <  1:3734603/71‑1 (MQ=255) gTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTa                                                                   >  1:521910/1‑70 (MQ=255)  tGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa                                                                  >  1:872901/1‑70 (MQ=255)    aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAat                                                               <  1:1637887/71‑1 (MQ=255)     gTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                                              <  1:3885647/71‑1 (MQ=255)     gTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                                              >  1:1974156/1‑71 (MQ=255)      tATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAatat                                                             <  1:1047257/71‑1 (MQ=255)      tATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAatat                                                             >  1:808394/1‑71 (MQ=255)      tATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAatat                                                             <  1:2942021/71‑1 (MQ=255)         cGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTc                                                           >  1:3907370/1‑70 (MQ=255)         cGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTc                                                           >  1:3902917/1‑70 (MQ=255)                ttttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGAttt                                                   >  1:3281715/1‑71 (MQ=255)                ttttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGAttt                                                   <  1:3918520/71‑1 (MQ=255)                  ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                                                  <  1:1455450/70‑1 (MQ=255)                  ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                                                  <  1:1724126/70‑1 (MQ=255)                  ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                                                  <  1:3659346/70‑1 (MQ=255)                     gTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                                                >  1:1556433/1‑53 (MQ=255)                         atgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGAtcgc                                            <  1:594407/66‑1 (MQ=255)                         atgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGAtcgc                                            <  1:2110966/66‑1 (MQ=255)                          ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCaa                                          >  1:3571166/1‑70 (MQ=255)                           tGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCaa                                          >  1:2071687/1‑69 (MQ=255)                           gggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTg                                                  >  1:3637716/2‑61 (MQ=255)                            ggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACg                                        <  1:3413146/70‑1 (MQ=255)                            ggTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACg                                        <  1:1588269/70‑1 (MQ=255)                              tGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGata                                     <  1:479848/71‑3 (MQ=255)                               gagaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                                     >  1:2183757/1‑70 (MQ=255)                                agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                                   <  1:1031011/71‑1 (MQ=255)                                  aTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGc                                                 >  1:2157446/1‑55 (MQ=255)                                  aTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt                                 <  1:2289002/71‑1 (MQ=255)                                   tAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTa                                >  1:1461850/1‑71 (MQ=255)                                   tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                                              <  1:47862/41‑1 (MQ=255)                                   tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                                              <  1:2343091/41‑1 (MQ=255)                                   tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                                              <  1:43173/41‑1 (MQ=255)                                   tAGAGAATCAGCAGGACACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                <  1:220749/71‑1 (MQ=255)                                    agagAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcgc                                            <  1:2383786/58‑1 (MQ=255)                                     gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTata                              <  1:2013045/71‑3 (MQ=255)                                     gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                              <  1:1765787/71‑1 (MQ=255)                                     gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                              <  1:1366534/71‑1 (MQ=255)                                     gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                              <  1:3595261/71‑1 (MQ=255)                                      agaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                                   <  1:3048127/65‑1 (MQ=255)                                      agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATAtt                                                            >  1:1788054/1‑40 (MQ=255)                                      agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATAtt                                                            >  1:2882151/1‑40 (MQ=255)                                      agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATAtt                                                            >  1:3915264/1‑40 (MQ=255)                                       gaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatatc                            >  1:665482/1‑70 (MQ=255)                                       gaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatatc                            <  1:3408650/71‑2 (MQ=255)                                       gaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAAt                             <  1:1390922/70‑1 (MQ=255)                                            agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                                               >  1:2611027/1‑47 (MQ=255)                                            agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                       <  1:954568/71‑1 (MQ=255)                                               agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                    >  1:574366/1‑71 (MQ=255)                                                gggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGcgcgcg                     >  1:3215920/1‑69 (MQ=255)                                                gggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGcgcgcg                     >  1:1412661/1‑69 (MQ=255)                                                aggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGc                   <  1:2253888/70‑1 (MQ=255)                                                aggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGc                   <  1:742297/70‑1 (MQ=255)                                                 ggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCt                  <  1:1765940/71‑1 (MQ=255)                                                 ggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCt                  <  1:1058675/71‑1 (MQ=255)                                                  gCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACg                                        <  1:3558903/48‑1 (MQ=255)                                                  gCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtg                 <  1:523346/71‑1 (MQ=255)                                                  gCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCt                  >  1:925881/1‑70 (MQ=255)                                                   cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgc                                              >  1:2350882/1‑41 (MQ=255)                                                   cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                                               >  1:198889/1‑40 (MQ=255)                                                   cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                <  1:3820351/71‑1 (MQ=255)                                                    aCGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGAc                                      >  1:2373244/1‑48 (MQ=255)                                                    aCGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtgt                >  1:2978068/1‑70 (MQ=255)                                                     cGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgc                          >  1:1591018/1‑59 (MQ=255)                                                     cGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtgt                >  1:2146594/1‑69 (MQ=255)                                                     cGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGa              <  1:2396610/71‑1 (MQ=255)                                                     cGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGa              <  1:1320841/71‑1 (MQ=255)                                                      gCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGAt             >  1:1978969/1‑71 (MQ=255)                                                      gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTa                                <  1:1194229/52‑1 (MQ=255)                                                      gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTa                                <  1:2916924/52‑1 (MQ=255)                                                      gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag             >  1:3258357/1‑70 (MQ=255)                                                      gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag             >  1:1627916/1‑70 (MQ=255)                                                      gCATGCAGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag             >  1:513700/1‑70 (MQ=255)                                                        aTGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGc           >  1:954552/1‑71 (MQ=255)                                                         tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag             <  1:3130043/68‑2 (MQ=255)                                                           ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtg        >  1:3201955/1‑71 (MQ=255)                                                           ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt         <  1:3627767/70‑1 (MQ=255)                                                           ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt         <  1:844890/70‑1 (MQ=255)                                                           ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt         <  1:970811/70‑1 (MQ=255)                                                           ccGTAAAAGCTAAAATATTCCCAGATTTGAAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt         <  1:2950957/70‑1 (MQ=255)                                                            cGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAAt                             <  1:428715/49‑1 (MQ=255)                                                             gTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGatat                                    <  1:1236708/41‑1 (MQ=255)                                                              tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                                     <  1:3207931/39‑1 (MQ=255)                                                              tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                                     <  1:3026772/39‑1 (MQ=255)                                                              tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGttt     <  1:1840787/71‑1 (MQ=255)                                                              tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg            <  1:2582739/64‑3 (MQ=255)                                                                aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttca  <  1:3945306/71‑4 (MQ=255)                                                                aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttca  >  1:975180/1‑68 (MQ=255)                                                                 aaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAtt  <  1:1835985/71‑1 (MQ=255)                                                                  aGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa                                   >  1:1545148/1‑37 (MQ=255)                                                                  aGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                <  1:3140237/56‑1 (MQ=255)                                                                   |                                                                     GTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATT  >  minE/454947‑455081

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