Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	455,034	C→G	31.2%	L17L (CTG→CTC) ‡	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	455,034	0	C	G	31.2%	106.1 / 69.8	87	L17L (CTG→CTC) ‡	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (23/36);  new base G (9/18);  total (32/55)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

ATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAT  >  minE/454971‑455098                                                                |                                                                 atgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGAtcgc                                                             <  1:2110966/66‑1 (MQ=255) atgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGAtcgc                                                             <  1:594407/66‑1 (MQ=255)  ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCaa                                                           >  1:3571166/1‑70 (MQ=255)   tGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCaa                                                           >  1:2071687/1‑69 (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:1366534/71‑1 (MQ=255)             gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                               <  1:3595261/71‑1 (MQ=255)              agaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                                                    <  1:3048127/65‑1 (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:1058675/71‑1 (MQ=255)                         ggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCt                                   <  1:1765940/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: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:970811/70‑1 (MQ=255)                                   ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt                          <  1:3627767/70‑1 (MQ=255)                                   ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt                          <  1:844890/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)                                             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: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)                                                                |                                                                 ATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAT  >  minE/454971‑455098

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