Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	455,045	C→T	25.0%	V14I (GTC→ATC)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	455,045	0	C	T	25.0%	141.3 / 45.8	80	V14I (GTC→ATC)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (26/34);  new base T (9/11);  total (35/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCAC  >  minE/454977‑455113                                                                     |                                                                     gagaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                                                                     >  1:2183757/1‑70 (MQ=255)  agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                                                                   <  1:1031011/71‑1 (MQ=255)    aTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt                                                                 <  1:2289002/71‑1 (MQ=255)     tAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTa                                                                >  1:1461850/1‑71 (MQ=255)     tAGAGAATCAGCAGGACACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                                                <  1:220749/71‑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)         gaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatatc                                                            >  1:665482/1‑70 (MQ=255)         gaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatatc                                                            <  1:3408650/71‑2 (MQ=255)         gaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAAt                                                             <  1:1390922/70‑1 (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: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)                    gCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtg                                                 <  1:523346/71‑1 (MQ=255)                    gCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCt                                                  >  1:925881/1‑70 (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:1320841/71‑1 (MQ=255)                       cGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGa                                              <  1:2396610/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:1627916/1‑70 (MQ=255)                        gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag                                             >  1:3258357/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:3627767/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:3026772/39‑1 (MQ=255)                                tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                                                                     <  1:3207931/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:834973/36‑1 (MQ=255)                                        aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTa                                                                <  1:2826580/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)                                            tatTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                >  1:1749513/1‑48 (MQ=255)                                             atTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAAtt                        >  1:1736829/1‑71 (MQ=255)                                              ttCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTc                       >  1:3255694/1‑71 (MQ=255)                                               tCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTc                       >  1:2045534/1‑70 (MQ=255)                                                 ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcac                    <  1:3367121/71‑1 (MQ=255)                                                 ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCca                     <  1:2367512/70‑1 (MQ=255)                                                   aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcaca                   >  1:3313784/1‑70 (MQ=255)                                                   aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACaa                  <  1:775624/71‑1 (MQ=255)                                                    gATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGATTATTCATAATAATTCCACAAt                 >  1:2068334/1‑71 (MQ=255)                                                         gCAGCGCATCGACATAGTAAATAGCGCTCGGCTGTGATGCGTGTTTa                                    <  1:3225894/47‑1 (MQ=39)                                                          cAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcaca                   >  1:2811476/1‑63 (MQ=255)                                                           aGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACgg           >  1:2086334/1‑70 (MQ=255)                                                              gcAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAAc       >  1:1308397/1‑71 (MQ=255)                                                                aaCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGt     >  1:197753/1‑71 (MQ=255)                                                                aaCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGt     >  1:854809/1‑71 (MQ=255)                                                                aaCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGt     >  1:316874/1‑71 (MQ=255)                                                                  cGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAAtt                        >  1:2443820/1‑50 (MQ=255)                                                                  cGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAATCATAATAATTCCACAATTGACGGCa         >  1:2142190/1‑65 (MQ=255)                                                                   gACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCAc  >  1:149303/1‑71 (MQ=255)                                                                     |                                                                     GAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCAC  >  minE/454977‑455113

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