Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	455,019	A→G	35.0%	Y22Y (TAT→TAC)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	455,019	0	A	G	35.0%	77.2 / 74.9	80	Y22Y (TAT→TAC)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (22/30);  new base G (10/18);  total (32/48)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAA  >  minE/454950‑455089                                                                      |                                                                       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:2942021/71‑1 (MQ=255)   tATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAatat                                                                     >  1:808394/1‑71 (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:3659346/70‑1 (MQ=255)               ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                                                          <  1:1724126/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:2343091/41‑1 (MQ=255)                                tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                                                      <  1:43173/41‑1 (MQ=255)                                tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                                                      <  1:47862/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:1366534/71‑1 (MQ=255)                                  gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                      <  1:3595261/71‑1 (MQ=255)                                  gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                      <  1:1765787/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:3408650/71‑2 (MQ=255)                                    gaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatatc                                    >  1:665482/1‑70 (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:1412661/1‑69 (MQ=255)                                             gggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGcgcgcg                             >  1:3215920/1‑69 (MQ=255)                                             aggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGc                           <  1:742297/70‑1 (MQ=255)                                             aggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGc                           <  1:2253888/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:2916924/52‑1 (MQ=255)                                                   gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTa                                        <  1:1194229/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:844890/70‑1 (MQ=255)                                                        ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt                 <  1:970811/70‑1 (MQ=255)                                                        ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt                 <  1:3627767/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:975180/1‑68 (MQ=255)                                                             aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttca          <  1:3945306/71‑4 (MQ=255)                                                              aaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAtt          <  1:1835985/71‑1 (MQ=255)                                                               aGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa                                           >  1:1545148/1‑37 (MQ=255)                                                               aGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                        <  1:3140237/56‑1 (MQ=255)                                                                  tAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                        <  1:3642785/53‑1 (MQ=255)                                                                  tAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa            <  1:2425680/65‑1 (MQ=255)                                                                   aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTa                                        <  1:2826580/36‑1 (MQ=255)                                                                   aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTa                                        <  1:834973/36‑1 (MQ=255)                                                                   aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                            <  1:2976604/48‑1 (MQ=255)                                                                   aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                            <  1:3761966/48‑1 (MQ=255)                                                                    aaaTATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                      <  1:3947282/37‑1 (MQ=255)                                                                    aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                      >  1:3120921/1‑53 (MQ=255)                                                                    aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                      >  1:768312/1‑53 (MQ=255)                                                                      atatTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa            <  1:2360379/61‑1 (MQ=255)                                                                      atatTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCAtaataa  <  1:1649166/71‑1 (MQ=255)                                                                      |                                                                       AGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAA  >  minE/454950‑455089

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