Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	455,022	T→C	33.3%	E21E (GAA→GAG)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	455,022	0	T	C	33.3%	88.4 / 73.8	84	E21E (GAA→GAG)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base T (21/35);  new base C (10/18);  total (31/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01

TATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTC  >  minE/454952‑455092                                                                       |                                                                       tATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAatat                                                                        <  1:1047257/71‑1 (MQ=255) tATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAatat                                                                        <  1:2942021/71‑1 (MQ=255) tATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAatat                                                                        >  1:808394/1‑71 (MQ=255)    cGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTc                                                                      >  1:3902917/1‑70 (MQ=255)    cGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTc                                                                      >  1:3907370/1‑70 (MQ=255)           ttttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGAttt                                                              >  1:3281715/1‑71 (MQ=255)           ttttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGAttt                                                              <  1:3918520/71‑1 (MQ=255)             ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                                                             <  1:3659346/70‑1 (MQ=255)             ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                                                             <  1:1724126/70‑1 (MQ=255)             ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                                                             <  1:1455450/70‑1 (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)                              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:3595261/71‑1 (MQ=255)                                gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                         <  1:1366534/71‑1 (MQ=255)                                 agaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                                              <  1:3048127/65‑1 (MQ=255)                                 agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATAtt                                                                       >  1:1788054/1‑40 (MQ=255)                                 agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATAtt                                                                       >  1:3915264/1‑40 (MQ=255)                                 agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATAtt                                                                       >  1:2882151/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:3215920/1‑69 (MQ=255)                                           gggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGcgcgcg                                >  1:1412661/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: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:970811/70‑1 (MQ=255)                                                      ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt                    <  1:844890/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:3761966/48‑1 (MQ=255)                                                                 aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                               <  1:2976604/48‑1 (MQ=255)                                                                  aaaTATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                         <  1:3947282/37‑1 (MQ=255)                                                                  aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                         >  1:768312/1‑53 (MQ=255)                                                                  aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                         >  1:3120921/1‑53 (MQ=255)                                                                    atatTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa               <  1:2360379/61‑1 (MQ=255)                                                                    atatTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCAtaataa     <  1:1649166/71‑1 (MQ=255)                                                                     tatTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                           >  1:1749513/1‑48 (MQ=255)                                                                      atTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAAtt   >  1:1736829/1‑71 (MQ=255)                                                                       ttCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTc  >  1:3255694/1‑71 (MQ=255)                                                                       |                                                                       TATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTC  >  minE/454952‑455092

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