Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,936,917 T→C 18.2% K226R (AAA→AGA)  yjgB ← predicted alcohol dehydrogenase, Zn‑dependent and NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,936,9170TC18.2% 92.1 / 15.5 44K226R (AAA→AGA) yjgBpredicted alcohol dehydrogenase, Zn‑dependent and NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (19/17);  new base C (6/2);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCTGCTCTT  >  minE/2936847‑2936983
                                                                      |                                                                  
aaaaTAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCtt                                                                    <  1:880136/71‑1 (MQ=255)
aaaaTAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCtt                                                                    >  1:3371578/1‑71 (MQ=255)
aaaaTAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCtt                                                                    >  1:3878720/1‑71 (MQ=255)
aaaaTAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCtt                                                                    >  1:3615124/1‑71 (MQ=255)
 aaaTAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCttt                                                                   <  1:691414/71‑1 (MQ=255)
 aaaTAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCttt                                                                   <  1:1284408/71‑1 (MQ=255)
  aaTAGGGCTGCCGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTCTCa                                                                 <  1:3456229/71‑1 (MQ=255)
    tAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAg                                                                <  1:1516699/71‑1 (MQ=255)
     aGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGt                                                               >  1:2123782/1‑71 (MQ=255)
            ccAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTg                                                          <  1:3738308/69‑1 (MQ=255)
            ccAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCg                                                        >  1:2858572/1‑71 (MQ=255)
              aGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGa                                                      >  1:1969579/1‑71 (MQ=255)
              aGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTCTCAGTGCCTGCGGa                                                      >  1:3788600/1‑71 (MQ=255)
                tCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTCTCAGTGCCTGCGGATc                                                    >  1:2739859/1‑71 (MQ=255)
                  gAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATcgc                                                  <  1:1639651/71‑1 (MQ=255)
                  gAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATcgc                                                  <  1:449966/71‑1 (MQ=255)
                  gAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATcgc                                                  >  1:2350762/1‑71 (MQ=255)
                  gAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTCTCAGTGCCTGCGGATcgc                                                  >  1:2694281/1‑71 (MQ=255)
                      cTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGc                                               <  1:700935/70‑1 (MQ=255)
                      cTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGc                                               <  1:1498162/70‑1 (MQ=255)
                         aCGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGc                                                         >  1:515154/1‑57 (MQ=255)
                         aCGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTAtt                                           <  1:553209/71‑1 (MQ=255)
                         aCGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTAt                                            <  1:103053/70‑1 (MQ=255)
                             tGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGc                                               >  1:110972/1‑63 (MQ=255)
                             tGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTcgc                                        >  1:117323/1‑68 (MQ=255)
                             tGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTcac                                        >  1:2670610/1‑70 (MQ=255)
                             tGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTCTCAGTGCCTGCGGATCGCGGCt                                              <  1:2525533/64‑1 (MQ=255)
                               aCGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTcacca                                      <  1:2150580/70‑1 (MQ=255)
                                cGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACt                                    >  1:1655222/1‑71 (MQ=255)
                                 ggTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACtt                                   <  1:2873037/71‑1 (MQ=255)
                                     ttGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATc                               <  1:384529/71‑1 (MQ=255)
                                       gATAATGAGATCAAACTGCCCCGCCAGTGCTCTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCgg                             >  1:347786/1‑71 (MQ=255)
                                         tAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATcgcg                                                 >  1:2242429/1‑49 (MQ=255)
                                          aaTGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCAc                          >  1:1297952/1‑71 (MQ=255)
                                               gaTCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATc                     >  1:1287691/1‑71 (MQ=255)
                                                  cAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGc                   >  1:3354982/1‑70 (MQ=255)
                                                  cAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGc                   >  1:3404720/1‑70 (MQ=255)
                                                  cAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGc                   >  1:1137046/1‑70 (MQ=255)
                                                   aaaCTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACtt                                   <  1:3434547/53‑1 (MQ=255)
                                                      cTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATc                               >  1:3230742/1‑54 (MQ=255)
                                                      cTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATc                     <  1:3708980/64‑1 (MQ=255)
                                                      cTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCa              >  1:2579728/1‑71 (MQ=255)
                                                         ccccGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACtt           >  1:2442110/1‑71 (MQ=255)
                                                         ccccGCCAGTGCTCTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCAc                          >  1:3823149/1‑56 (MQ=255)
                                                         ccccGCCAGTGCTCTCAGTGCCTGAGGATCGCGGCTATTCACCACTTTATCGGCAc                          >  1:2216787/1‑56 (MQ=39)
                                                            cGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCAc                          >  1:1695106/1‑53 (MQ=255)
                                                             gCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCTg       <  1:3686820/71‑1 (MQ=255)
                                                              ccAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCTGc      <  1:1098934/71‑1 (MQ=255)
                                                                  tGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCTGCTCtt  <  1:2464732/71‑1 (MQ=255)
                                                                  tGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCTGCTCtt  <  1:1066186/71‑1 (MQ=255)
                                                                      |                                                                  
AAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCTGCTCTT  >  minE/2936847‑2936983

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: