Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,855,148 G→A 27.1% G194G (GGC→GGT yjeS ← predicted Fe‑S electron transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,855,1480GA27.1% 98.9 / 34.0 88G194G (GGC→GGTyjeSpredicted Fe‑S electron transport protein
Reads supporting (aligned to +/- strand):  ref base G (44/20);  new base A (7/17);  total (51/37)
Fisher's exact test for biased strand distribution p-value = 1.36e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

GAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAA  >  minE/2855081‑2855212
                                                                   |                                                                
gAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTg                                                                  >  1:1105771/1‑68 (MQ=255)
gAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTg                                                                  >  1:163987/1‑68 (MQ=255)
gAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGc                                                                 >  1:1983218/1‑69 (MQ=255)
gAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGc                                                                 >  1:2754713/1‑69 (MQ=255)
gAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGc                                                                 >  1:2371354/1‑69 (MQ=255)
gAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGc                                                                 >  1:2367294/1‑69 (MQ=255)
gAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGc                                                                 >  1:3361684/1‑69 (MQ=255)
gAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGc                                                                 >  1:416524/1‑69 (MQ=255)
     tCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCa                                                             <  1:222639/68‑1 (MQ=255)
       gACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCAtcc                                                          <  1:1196772/69‑1 (MQ=255)
       gACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCAtcc                                                          <  1:232846/69‑1 (MQ=255)
       gACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCAtcc                                                          <  1:1032617/69‑1 (MQ=255)
       gACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCAtcc                                                          <  1:3278133/69‑1 (MQ=255)
       gACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCAtcc                                                          <  1:3239159/69‑1 (MQ=255)
       gACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCAtcc                                                          <  1:1939939/69‑1 (MQ=255)
       gACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCAtcc                                                          <  1:1860301/69‑1 (MQ=255)
                tGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGAc                                                 >  1:2597160/1‑69 (MQ=255)
                tGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGAc                                                 >  1:2931522/1‑69 (MQ=255)
                tGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGAc                                                 >  1:1835567/1‑69 (MQ=255)
                tGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGAc                                                 >  1:2894155/1‑69 (MQ=255)
                tGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGAc                                                 >  1:1002773/1‑69 (MQ=255)
                          aTGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGAtac                                       <  1:1249320/69‑3 (MQ=255)
                          aTGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATcc                                       <  1:3112487/69‑1 (MQ=255)
                          aTGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATcc                                       <  1:2532289/69‑1 (MQ=255)
                          aTGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATcc                                       <  1:3227867/69‑1 (MQ=255)
                          aTGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATcc                                       <  1:2871364/69‑1 (MQ=255)
                          aTGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATcc                                       <  1:2865319/69‑1 (MQ=255)
                          aTGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATcc                                       <  1:2809574/69‑1 (MQ=255)
                          aTGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATcc                                       <  1:1047922/69‑1 (MQ=255)
                          aTGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATcc                                       <  1:1217954/69‑1 (MQ=255)
                          aTGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATcc                                       <  1:419036/69‑1 (MQ=255)
                          aTGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATcc                                       <  1:1493617/69‑1 (MQ=255)
                          aTGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATcc                                       <  1:364900/69‑1 (MQ=255)
                          aTGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATcc                                       <  1:1806411/69‑1 (MQ=255)
                          aTGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCCACTGGTTGATcc                                       <  1:291670/69‑1 (MQ=255)
                           tGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATcc                                       <  1:1944172/68‑1 (MQ=255)
                           tGGCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATcc                                       <  1:3221820/68‑1 (MQ=255)
                            ggCTCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATcc                                       <  1:647996/67‑1 (MQ=255)
                                     ctggCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCgg                            <  1:2519760/67‑1 (MQ=255)
                                     cGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCgg                            <  1:526118/69‑1 (MQ=255)
                                     cGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCgg                            <  1:3407542/69‑1 (MQ=255)
                                      gggCAAATTGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCgg                            <  1:273273/68‑1 (MQ=255)
                                               gTCATACAGGCCACGCGTTTACCGCATCCTTCCTc                                                    >  1:3084907/1‑35 (MQ=38)
                                               gTCATACAGGCCACGCATTTACCGCATCCTTCCTc                                                    >  1:394722/1‑35 (MQ=255)
                                               gTCATACAGGCCACGCATTTACCGCATCCTTCCTc                                                    >  1:641912/1‑35 (MQ=255)
                                               gTCATACAGGCCACGCATTTACCGCATCCTTCCTc                                                    >  1:245751/1‑35 (MQ=255)
                                               gTCATACAGGCCACGCATTTACCGCATCCTTCCTc                                                    >  1:1903048/1‑35 (MQ=255)
                                               gTCATACAGGCCACGCATTTACCGCATCCTTCCTc                                                    >  1:666498/1‑35 (MQ=255)
                                               gTCATACAGGCCACGCATTTAACGCATCCTTCCTc                                                    >  1:1745706/1‑35 (MQ=37)
                                                  aTACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGc               <  1:295013/69‑1 (MQ=255)
                                                  aTACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGc               <  1:3307469/69‑1 (MQ=255)
                                                  aTACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGc               <  1:934561/69‑1 (MQ=255)
                                                  aTACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGc               <  1:2681610/69‑1 (MQ=255)
                                                  aTACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGc               <  1:1008218/69‑1 (MQ=255)
                                                  aTACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGc               <  1:1680457/69‑1 (MQ=255)
                                                  aTACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGc               <  1:2403224/69‑1 (MQ=255)
                                                  aTACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGc               <  1:2438020/69‑1 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTgg                                              >  1:1883194/1‑35 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGAt                                         >  1:2620344/1‑40 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATc                      >  1:802175/1‑59 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCg                     >  1:2924184/1‑60 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCg                     >  1:594704/1‑60 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCg                     >  1:2652259/1‑60 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCg                     >  1:38543/1‑60 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCg                     >  1:383207/1‑60 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCg                     >  1:1978585/1‑60 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCg                     >  1:61778/1‑60 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCg                     >  1:627590/1‑60 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCg                     >  1:321768/1‑60 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCg                     >  1:2541640/1‑60 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCaa             >  1:645061/1‑68 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCaa             >  1:1702321/1‑68 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:1179080/1‑69 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:1118290/1‑69 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:762446/1‑69 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:810224/1‑69 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:2719063/1‑69 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:583501/1‑69 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:1474338/1‑69 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:1485228/1‑69 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:2119757/1‑69 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:2343396/1‑69 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:2409073/1‑69 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:2573547/1‑69 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:2979006/1‑69 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:2563941/1‑69 (MQ=255)
                                                     cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:2570562/1‑69 (MQ=255)
                                                     cAGGCCACGCATTTGCCGAATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAAt            >  1:1115920/1‑69 (MQ=255)
                                                     cAGGACACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGc               >  1:2546673/1‑66 (MQ=255)
                                                                  tGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTaaaa  >  1:2466993/1‑66 (MQ=255)
                                                                   |                                                                
GAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAA  >  minE/2855081‑2855212

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