Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 851,042 G→A 50.0% R53R (CGC→CGT yciA ← predicted hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE851,0420GA50.0% ‑4.2 / 45.2 52R53R (CGC→CGTyciApredicted hydrolase
Reads supporting (aligned to +/- strand):  ref base G (14/12);  new base A (23/3);  total (37/15)
Fisher's exact test for biased strand distribution p-value = 1.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

AGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCA  >  minE/850975‑851104
                                                                   |                                                              
agcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAcgc                                                               >  1:1151886/1‑69 (MQ=255)
agcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAcgc                                                               >  1:2731749/1‑69 (MQ=255)
       cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACGGTCACTACACGAc                                                            <  1:1120890/65‑1 (MQ=255)
        acaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACGGTCACTACACGAc                                                            <  1:1493995/64‑1 (MQ=255)
        acaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACGGTCACTACACGAc                                                            <  1:2471109/64‑1 (MQ=255)
                         ctgCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                                      <  1:3422524/67‑1 (MQ=255)
                         cGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                                      <  1:3100567/69‑1 (MQ=255)
                         cGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                                      <  1:1417919/69‑1 (MQ=255)
                         cGGCCGTAAGAAAGTCATTCCTTCAACCCCCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                                      <  1:1065536/69‑1 (MQ=255)
                         cGGCCGTAAGAAAGTCATTCCTACAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                                      <  1:1058667/69‑1 (MQ=255)
                         cGGCCGTAAGAAAGTCATTCCATCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                                      <  1:2242693/69‑1 (MQ=255)
                          ggCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                                      >  1:1065227/1‑68 (MQ=255)
                          ggCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                                      >  1:1779479/1‑68 (MQ=255)
                            ccGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAAt                                   <  1:196863/69‑1 (MQ=255)
                                       tCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCaa                        >  1:2589372/1‑69 (MQ=255)
                                       tCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCaa                        >  1:2163290/1‑69 (MQ=255)
                                        cATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAat                       >  1:3182290/1‑69 (MQ=255)
                                        cATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAat                       >  1:415248/1‑69 (MQ=255)
                                        cATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAat                       >  1:787755/1‑69 (MQ=255)
                                        cATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAat                       >  1:2445577/1‑69 (MQ=255)
                                        cATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAat                       >  1:2019077/1‑69 (MQ=255)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATgggg                                                     >  1:721540/1‑35 (MQ=25)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:1904338/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:1008285/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:778052/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:553606/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:1249013/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:1266828/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:3269759/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:136036/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:3016844/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:3006899/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:3000926/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:2922481/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:1434385/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:247747/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:1793633/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:1976684/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:1507194/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:1633796/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:2379315/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:1605331/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:253272/1‑36 (MQ=37)
                                           tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                                     >  1:1535841/1‑36 (MQ=37)
                                            ccTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATcc                   >  1:1897177/1‑69 (MQ=255)
                                            ccTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATcc                   >  1:2849274/1‑69 (MQ=255)
                                                      gCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACa         <  1:30582/69‑1 (MQ=255)
                                                      gCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACa         <  1:218347/69‑1 (MQ=255)
                                                      gCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACa         <  1:234321/69‑1 (MQ=255)
                                                      gCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACa         <  1:1490792/69‑1 (MQ=255)
                                                      gCACAGTCACAACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACa         <  1:2878724/69‑1 (MQ=255)
                                                             cACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTaacca  >  1:1423009/1‑69 (MQ=255)
                                                             cACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTaacca  >  1:2427618/1‑69 (MQ=255)
                                                             cACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTaacca  >  1:1492905/1‑69 (MQ=255)
                                                                   gcgACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACa         <  1:3037269/56‑1 (MQ=255)
                                                                   |                                                              
AGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCA  >  minE/850975‑851104

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: