Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,357,646 G→A 55.7% S492S (AGC→AGT yehB ← predicted outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,357,6460GA55.7% 8.4 / 47.5 61S492S (AGC→AGTyehBpredicted outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (11/16);  new base A (24/10);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 3.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

ACGTACTTAATGACACAGACCCCCAACCTTCTGGCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATA  >  minE/1357603‑1357711
                                           |                                                                 
aCGTACTTAATGACACAGACCCCCAACCTTCTGGCAATGACTgg                                                                   >  1:2530222/1‑44 (MQ=255)
aCGTACTTAATGACACAGACCCCCAACCTTCTGGCAATGACTGGCTCATATTGGCGGAAAAGCTAtttt                                          >  1:1885812/1‑69 (MQ=255)
aCGTACTTAATGACACAGACCCCCAACCTTCTGGCAATGACTGGCTCATATTGGCGGAAAAGCTAtttt                                          >  1:3389586/1‑69 (MQ=255)
  gTACTTAAGGACACAGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTAttttt                                         <  1:3196040/68‑1 (MQ=255)
  gTACTTAAGGACACAGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTAttttt                                         <  1:831685/68‑1 (MQ=255)
  gTACTTAAGGACACAGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTAttttt                                         <  1:2927100/68‑1 (MQ=255)
  gTACTTAAGGACACAGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTAttttt                                         <  1:3115287/68‑1 (MQ=255)
  gTACTTAAGGACACAGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTAttttt                                         <  1:1286219/68‑1 (MQ=255)
  gTACTTAAGGACACAGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTAttttt                                         <  1:115224/68‑1 (MQ=255)
  gTACTTAAGGACACAGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTAttttt                                         <  1:2499152/68‑1 (MQ=255)
  gTACTTAAGGACACAGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTAttttt                                         <  1:2241261/68‑1 (MQ=255)
   tACTTAAGGACACAGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTAttttt                                         <  1:3090883/67‑1 (MQ=255)
    aCTTAAGGACACAGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTAttttt                                         <  1:466182/66‑1 (MQ=255)
       tAATGACACAGACCCCCAACCTTCTGGCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGc                                    >  1:3221267/1‑68 (MQ=255)
       tAATGACACAGACCCCCAACCTTCTGGCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGc                                    >  1:2151817/1‑68 (MQ=255)
       tAATGACACAGACCCCCAACCTTCTGGCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGc                                    >  1:1118018/1‑68 (MQ=255)
             cacaGAACCCCAGGCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATc                             >  1:546988/1‑69 (MQ=255)
             cacaGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGaaa                                                   >  1:853398/1‑47 (MQ=25)
             cacaGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATc                             >  1:2695743/1‑69 (MQ=255)
             cacaGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATc                             >  1:2696574/1‑69 (MQ=255)
             cacaGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATc                             >  1:2546857/1‑69 (MQ=255)
             cacaGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATc                             >  1:2911088/1‑69 (MQ=255)
             cacaGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATc                             >  1:2997955/1‑69 (MQ=255)
             cacaGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATc                             >  1:3244417/1‑69 (MQ=255)
             cacaGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATc                             >  1:3355810/1‑69 (MQ=255)
             cacaGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATc                             >  1:2239027/1‑69 (MQ=255)
             cacaGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATc                             >  1:1505019/1‑69 (MQ=255)
             cacaGAACCCCAGCCTTCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATc                             >  1:1578874/1‑69 (MQ=255)
              acaGACCCCCAACCTTCTGGCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCg                            >  1:1517152/1‑69 (MQ=255)
                     cccAACCTTCTGGCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGt                     <  1:482071/69‑1 (MQ=255)
                     cccAACCTTCTGGCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGt                     <  1:2547312/69‑1 (MQ=255)
                     cccAACCTTCTGGCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGt                     <  1:839273/69‑1 (MQ=255)
                     cccAACCTTCTGGCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGt                     <  1:2919100/69‑1 (MQ=255)
                          ccTTCTGGCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGtaataa                >  1:945887/1‑69 (MQ=255)
                          ccTTCTGGCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGtaataa                >  1:1087366/1‑69 (MQ=255)
                          ccTTCTGGCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGtaataa                >  1:700766/1‑69 (MQ=255)
                            ttCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCaaa                                >  1:663810/1‑51 (MQ=255)
                            ttCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATCGTTCTGGTAATAATc              >  1:2306662/1‑69 (MQ=255)
                            ttCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATCGTTCTGGTAATAATc              >  1:1355142/1‑69 (MQ=255)
                            ttCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATCGTTCTGGTAATAATc              >  1:400223/1‑69 (MQ=255)
                            ttCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATCGTTCTGGTAATAATc              >  1:1518613/1‑69 (MQ=255)
                            ttCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATCGTTCTGGTAATAATc              >  1:152954/1‑69 (MQ=255)
                            ttCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATCGTTCTGGTAATAATc              >  1:2788566/1‑69 (MQ=255)
                            ttCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATCGTTCTGGTAATAATc              >  1:1720188/1‑69 (MQ=255)
                            ttCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATCGTTCTGGTAATAATc              >  1:1889626/1‑69 (MQ=255)
                            ttCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATCGTTCTGGTAATAATc              >  1:2535246/1‑69 (MQ=255)
                            ttCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATCGTTCTGGTAATAATc              >  1:2156822/1‑69 (MQ=255)
                            ttCTGGCAATGACTGACTCATATTGGCAGAAAAGCTATTTTTACGGCCAAAATCGTTCTGGTAATAATc              >  1:2303503/1‑69 (MQ=255)
                                ggCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGc           <  1:947429/68‑1 (MQ=255)
                                ggCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGc           <  1:3235645/68‑1 (MQ=255)
                                ggCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGc           <  1:3046458/68‑1 (MQ=255)
                                ggCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGc           <  1:2629943/68‑1 (MQ=255)
                                ggCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGc           <  1:2453492/68‑1 (MQ=255)
                                ggCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGc           <  1:2273310/68‑1 (MQ=255)
                                ggCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGc           <  1:1993626/68‑1 (MQ=255)
                                ggCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGc           <  1:164335/68‑1 (MQ=255)
                                ggCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGc           <  1:1122404/68‑1 (MQ=255)
                                ggCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGc           <  1:1054130/68‑1 (MQ=255)
                                 gCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGc           <  1:14537/67‑1 (MQ=255)
                                     tGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGt      >  1:1960271/1‑68 (MQ=255)
                                     tGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCAAAAATCGTTCTGGTAATAATCGGCAATGt      >  1:3340320/1‑68 (MQ=255)
                                         tGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATa  <  1:1009898/68‑1 (MQ=255)
                                           |                                                                 
ACGTACTTAATGACACAGACCCCCAACCTTCTGGCAATGACTGGCTCATATTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATA  >  minE/1357603‑1357711

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: