Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	455,055	A→C	19.5%	A10A (GCT→GCG)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	455,055	0	A	C	19.5%	115.6 / 38.4	72	A10A (GCT→GCG)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (30/28);  new base C (7/7);  total (37/35)
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.41e-01

AGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGT  >  minE/454990‑455121                                                                  |                                                                   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:1058675/71‑1 (MQ=255)      ggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCt                                                          <  1:1765940/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)         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)           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: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)                   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)                       aGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                        <  1:3140237/56‑1 (MQ=255)                          tAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                        <  1:3642785/53‑1 (MQ=255)                          tAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa                                            <  1:2425680/65‑1 (MQ=255)                           aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                                                            <  1:2976604/48‑1 (MQ=255)                           aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                                                            <  1:3761966/48‑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)                                  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:316874/1‑71 (MQ=255)                                                   aaCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGt             >  1:854809/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)                                                         aTAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGt        >  1:1375216/1‑70 (MQ=255)                                                          tAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa      <  1:2576594/71‑1 (MQ=255)                                                          tAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa      <  1:1664833/71‑1 (MQ=255)                                                             tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATa    <  1:3792915/70‑1 (MQ=255)                                                             tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATa    <  1:591364/70‑1 (MQ=255)                                                             tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAg   >  1:2831583/1‑71 (MQ=255)                                                             tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAg   <  1:855418/71‑1 (MQ=255)                                                             tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAg   >  1:1353590/1‑71 (MQ=255)                                                              aaaTAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGt  >  1:1821709/1‑71 (MQ=255)                                                               aaTAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa                      >  1:3001603/1‑50 (MQ=255)                                                               aaTAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGt  >  1:2918522/1‑70 (MQ=255)                                                                 tAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTc                               <  1:1720520/39‑1 (MQ=255)                                                                  aGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGt        >  1:35047/1‑61 (MQ=255)                                                                  |                                                                   AGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGT  >  minE/454990‑455121

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