Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,356,788 G→A 41.1% G778G (GGC→GGT yehB ← predicted outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,356,7880GA41.1% 51.8 / 84.1 119G778G (GGC→GGTyehBpredicted outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (59/11);  new base A (24/25);  total (83/36)
Fisher's exact test for biased strand distribution p-value = 8.43e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

GCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATCTGCTC  >  minE/1356723‑1356852
                                                                 |                                                                
gCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                                >  1:2305260/1‑68 (MQ=255)
gCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                                >  1:3141242/1‑68 (MQ=255)
gCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                                >  1:1901040/1‑68 (MQ=255)
gCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                                >  1:1245179/1‑68 (MQ=255)
gCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                                >  1:1221167/1‑68 (MQ=255)
gCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                                >  1:1196243/1‑68 (MQ=255)
gCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                                >  1:2194755/1‑68 (MQ=255)
gCCACATTAACCGATGGCGGTAATTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                                >  1:1080600/1‑68 (MQ=255)
           cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                                    >  1:764061/1‑69 (MQ=255)
           cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                                    <  1:3098399/69‑1 (MQ=255)
           cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                                    <  1:1840741/69‑1 (MQ=255)
           cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                                    >  1:3198238/1‑69 (MQ=255)
           cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                                    >  1:2956026/1‑69 (MQ=255)
           cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                                    >  1:3361085/1‑69 (MQ=255)
           cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                                    >  1:1506661/1‑69 (MQ=255)
           cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                                    >  1:1912687/1‑69 (MQ=255)
           cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                                    >  1:2693085/1‑69 (MQ=255)
           cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                                    >  1:2659810/1‑69 (MQ=255)
           cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                                    >  1:2638099/1‑69 (MQ=255)
           cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                                    >  1:804839/1‑69 (MQ=255)
           cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                                    >  1:2569197/1‑69 (MQ=255)
                     aCTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCa                                          >  1:1878420/1‑69 (MQ=255)
                     aCTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCa                                          >  1:1617332/1‑69 (MQ=255)
                     aCTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCa                                          >  1:315718/1‑69 (MQ=255)
                     aCTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCa                                          >  1:465554/1‑69 (MQ=255)
                     aCTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCa                                          >  1:3152147/1‑69 (MQ=255)
                      cTTCATTGGTGCTAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:539308/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGAGAACACCAATATTATGACCATGGATATCa                                          >  1:3073885/1‑68 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:2599338/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:3348389/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:3351704/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:2945617/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:423863/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:47849/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:636197/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:696276/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:740612/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:2527842/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:850853/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:874709/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:1097336/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:222785/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:1165359/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:1392046/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:2104761/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:1535949/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:2086121/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:2066708/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:1974189/47‑1 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATc                                           >  1:784994/1‑67 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATc                                           >  1:1711695/1‑67 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                                          >  1:3075841/1‑68 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                                          >  1:3227897/1‑68 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                                          >  1:1569791/1‑68 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                                          >  1:2270864/1‑68 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                                          >  1:1531107/1‑68 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                                          >  1:15306/1‑68 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                                          >  1:348941/1‑68 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                                          >  1:1411655/1‑68 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                                          >  1:623078/1‑68 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                                          >  1:635818/1‑68 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                                          >  1:1308377/1‑68 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                                          >  1:2490050/1‑68 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                                          >  1:1876992/1‑68 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                                          >  1:2674093/1‑68 (MQ=255)
                      cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                                          >  1:2902950/1‑68 (MQ=255)
                       ttCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                                               <  1:2472329/46‑1 (MQ=255)
                                 cGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCAtt                                        >  1:1160694/1‑59 (MQ=255)
                                 cGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCAtt                                        >  1:2785965/1‑59 (MQ=255)
                                 cGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCAtt                                        >  1:2469854/1‑59 (MQ=255)
                                     taaataaCTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGatat                                            >  1:3306812/1‑51 (MQ=255)
                                     taaataaCTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGatat                                            >  1:408426/1‑51 (MQ=255)
                                     taaataaCTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGatat                                            >  1:2684454/1‑51 (MQ=255)
                                     taaataaCTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGatat                                            >  1:3255215/1‑51 (MQ=255)
                                     taaataaCTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGatat                                            >  1:1895138/1‑51 (MQ=255)
                                     taaataaCTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGatat                                            >  1:1112243/1‑51 (MQ=255)
                                      aaataaCTGACTTCCCTGGCCGACAACGCCAatat                                                           <  1:2099650/35‑1 (MQ=255)
                                             tGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAAt                  <  1:2463426/69‑1 (MQ=255)
                                                     cTGGCCGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGcc           <  1:1167134/68‑1 (MQ=255)
                                                     cTGGCCGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGcc           <  1:45735/68‑1 (MQ=255)
                                                      tGGCCGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGccc          <  1:3378849/68‑1 (MQ=255)
                                                        gCCGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATc       >  1:635121/1‑69 (MQ=255)
                                                        gCCGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATc       >  1:573176/1‑69 (MQ=255)
                                                        gCCGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATc       >  1:215638/1‑69 (MQ=255)
                                                        gCCGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATc       >  1:2477098/1‑69 (MQ=255)
                                                        gCCGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATc       >  1:1833499/1‑69 (MQ=255)
                                                          cGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATCTg     <  1:801639/69‑1 (MQ=255)
                                                          cGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATCTg     <  1:2272797/69‑1 (MQ=255)
                                                          cGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATCTg     <  1:3091515/69‑1 (MQ=255)
                                                              aaCGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGt                      >  1:1765405/1‑48 (MQ=255)
                                                              aaCGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATCTGctc  >  1:1453885/1‑68 (MQ=255)
                                                              aaCGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATCTGctc  >  1:1314945/1‑68 (MQ=255)
                                                              aaCGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATCTGctc  >  1:673897/1‑68 (MQ=255)
                                                              aaCGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATCTGctc  >  1:1459354/1‑68 (MQ=255)
                                                              aaCGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATCTGctc  >  1:69794/1‑68 (MQ=255)
                                                              aaCGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATCTGctc  >  1:2629545/1‑68 (MQ=255)
                                                              aaCGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATCTGctc  >  1:3259983/1‑68 (MQ=255)
                                                              aaCGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATCTGctc  >  1:1639870/1‑68 (MQ=255)
                                                              aaCGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATCTGctc  >  1:1752191/1‑68 (MQ=255)
                                                              aaCGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATCTGctc  >  1:2958873/1‑68 (MQ=255)
                                                                 |                                                                
GCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAATGATTGCCCATCTGCTC  >  minE/1356723‑1356852

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