Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,833,531 C→T 41.3% R92C (CGC→TGC)  ygbL → predicted class II aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,833,5310CT41.3% 28.5 / 51.7 46R92C (CGC→TGC) ygbLpredicted class II aldolase
Reads supporting (aligned to +/- strand):  ref base C (14/13);  new base T (11/8);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGA  >  minE/1833463‑1833598
                                                                    |                                                                   
aaTGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCgtgc                                                                   <  1:2143878/71‑1 (MQ=255)
    gTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTa                                                               <  1:1469001/71‑1 (MQ=255)
               aCAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGTGCTGTAAAGCggtgg                                                      >  1:1999314/1‑69 (MQ=255)
               aCAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGTGCTGTAAAGCggtgg                                                      >  1:1544047/1‑69 (MQ=255)
               aCAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGTGCTGTAAAGCggtgg                                                      >  1:484380/1‑69 (MQ=255)
               aCAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGTGCTGTAAAGCggtgg                                                      >  1:486984/1‑69 (MQ=255)
               aCAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGTGCTGTAAAGCggtg                                                       >  1:412958/1‑68 (MQ=255)
               aCAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGTGCTGTAAAGCggtg                                                       >  1:571993/1‑68 (MQ=255)
                cAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGCTGGTGc                                                   <  1:738210/71‑1 (MQ=255)
                 aaaCCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCa                                                  >  1:926848/1‑71 (MQ=255)
                   aCCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGTGCTGTAAAGCGGTGGTGCAtt                                                >  1:524186/1‑71 (MQ=255)
                     ccTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGTGCTGTAAAGCGGTGGTGCATTTg                                              <  1:1026852/71‑1 (MQ=255)
                     ccTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGTGCTGTAAAGCGGTGGTGCATTTg                                              <  1:441724/71‑1 (MQ=255)
                       tCGAAAGATGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGc                                                   <  1:110350/64‑1 (MQ=255)
                         gAAAGAGGTGCTCTTTCATCTGGCTCTGTATCGCAACAATCCGTGCTGTAAAGCggtgg                                                      >  1:1794293/1‑59 (MQ=255)
                         gAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGTGCTGTAAAGCGGTGGTGCATTTgca                                            <  1:139737/69‑1 (MQ=255)
                         gAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgca                                            >  1:1139802/1‑69 (MQ=255)
                                 tGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGTGCTGTAAAGCGGTGGTGCATTTGCACAGCACATgg                                  <  1:1070975/71‑1 (MQ=255)
                                     ctTTCATCTGGCGCTGTATCGCAACAATCCTCGCTGTAAAGCGGTGGTGCATTTgcaca                                          <  1:1204445/59‑1 (MQ=255)
                                     ctTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGa                              >  1:243732/1‑71 (MQ=255)
                                      tttCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCggt                                                        <  1:682128/44‑1 (MQ=255)
                                      tttCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCggt                                                        <  1:1855314/44‑1 (MQ=255)
                                        tCATCTGGCGCTGTATCGCAACAATCCGTGCTGTAAAGCGGTGGTGCATTTgcacagcaca                                     >  1:94713/1‑61 (MQ=255)
                                        tCATCTGGCGCTGTATCGCAACAATCCGTGCTGTAAAGCGGTGGTGCATTTgcacagcaca                                     >  1:639439/1‑61 (MQ=255)
                                        tCATCTGGCGCTGTATCGCAACAATCCGTGCTGTAAAGCGGTGGTGCATTTgcacagcaca                                     >  1:2182581/1‑61 (MQ=255)
                                        tCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCggtg                                                       >  1:2161693/1‑43 (MQ=255)
                                        tCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACgg                           >  1:171382/1‑71 (MQ=255)
                                              ggCGCTGTATCGCAACAATCCGCGCTGTAAAGCggtggt                                                     >  1:1834796/1‑39 (MQ=255)
                                               gcgcTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCtt                      <  1:1278915/69‑1 (MQ=255)
                                                   tGTATCGCAACAATCCGTGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCttt                     <  1:2116278/66‑1 (MQ=255)
                                                     tATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCttt                     >  1:550471/1‑64 (MQ=255)
                                                       tCGCAACAATCCGTGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCctgcctgc            <  1:1474944/71‑1 (MQ=255)
                                                       tCGCAACAATCCGTGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCctgcctgc            <  1:1088157/71‑1 (MQ=255)
                                                       tCGCAACAATCCGCGCTGTAAAGCGGTTGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCctgcctgc            <  1:1206888/71‑1 (MQ=255)
                                                       tCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTcctgcctg             >  1:1481293/1‑70 (MQ=255)
                                                       tCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTcctgcctg             >  1:881752/1‑70 (MQ=255)
                                                       tCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTcctg                 <  1:1603360/66‑1 (MQ=255)
                                                       tCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCctgcctgc            <  1:302785/71‑1 (MQ=255)
                                                       tCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCctgcctgc            <  1:1731938/71‑1 (MQ=255)
                                                       tCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCctgcctgc            <  1:1248207/71‑1 (MQ=255)
                                                         gCAACAATCCGTGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCctgcctgc            <  1:503158/69‑1 (MQ=255)
                                                         gCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCctgcctgc            >  1:1030593/1‑69 (MQ=255)
                                                         gCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCctgcctgc            >  1:972241/1‑69 (MQ=255)
                                                         gCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCctgcctgc            >  1:125385/1‑69 (MQ=255)
                                                          gaacaaTCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAg         >  1:940765/2‑71 (MQ=255)
                                                            acaaTCCTCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAggg       <  1:136197/71‑1 (MQ=255)
                                                              aaTCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCt     >  1:245978/1‑71 (MQ=255)
                                                                  cgcgCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGa  >  1:978695/1‑70 (MQ=255)
                                                                    |                                                                   
AATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCTGGA  >  minE/1833463‑1833598

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: