Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,194,880 C→T 14.6% A102T (GCG→ACG)  yeaZ ← predicted peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,194,8800CT14.6% 129.0 / 13.3 48A102T (GCG→ACG) yeaZpredicted peptidase
Reads supporting (aligned to +/- strand):  ref base C (22/19);  new base T (5/2);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTT  >  minE/1194813‑1194948
                                                                   |                                                                    
cAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGcccc                                                                   <  1:4005155/71‑1 (MQ=255)
 aGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCt                                                                  >  1:303020/1‑71 (MQ=255)
    aaaCTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGTCCCTTgc                                                               >  1:1278235/1‑71 (MQ=255)
       cTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGcc                                                            <  1:2092190/71‑1 (MQ=255)
          cGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCa                                                           <  1:478016/69‑1 (MQ=255)
            cccATTCGCGCGTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGt                                                       <  1:395702/71‑1 (MQ=255)
            cccATTCGCGCGTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGt                                                       <  1:1680599/71‑1 (MQ=255)
             ccATTCGCGCGTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGt                                                       <  1:2684760/70‑1 (MQ=255)
              cATTCGCGCGTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCAt                                                          >  1:3676150/1‑66 (MQ=255)
                 tCGCGCGTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGTCCCTTGCGCCATCGTCAt                                                    >  1:1241144/1‑69 (MQ=255)
                    cgcgTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCg                                               >  1:1223588/1‑71 (MQ=255)
                    cgcgTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCg                                               >  1:3359520/1‑71 (MQ=255)
                     gcgTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGTCCCTTGCGCCATCGTCATTAGCGt                                              <  1:2251377/71‑1 (MQ=255)
                     gcgTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGt                                              >  1:1173109/1‑71 (MQ=255)
                      cgTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTg                                             <  1:3904161/71‑1 (MQ=255)
                      cgTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTg                                             >  1:469159/1‑71 (MQ=255)
                        tCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGga                                           <  1:1985630/71‑1 (MQ=255)
                        tCAATGGCTGCCAGCACGCGGGTTCCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGga                                           <  1:3568465/71‑1 (MQ=255)
                           aTGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGAc                                        >  1:800181/1‑71 (MQ=255)
                            tGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACg                                       >  1:442762/1‑71 (MQ=255)
                            tGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACg                                       >  1:1996267/1‑71 (MQ=255)
                              gCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGTCCCTTgcgc                                                             >  1:3609889/1‑47 (MQ=255)
                              gCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGTCCCTTgc                                                               >  1:3360053/1‑45 (MQ=255)
                               cTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCa                                    <  1:2601683/71‑1 (MQ=255)
                                tGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCaa                                   >  1:1705402/1‑71 (MQ=255)
                                    aGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAatcat                               <  1:3335520/71‑1 (MQ=255)
                                      cACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCAt                                                    <  1:1408902/48‑1 (MQ=255)
                                      cACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCg                             >  1:2247775/1‑71 (MQ=255)
                                        cgcgGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGt                           <  1:1298982/71‑1 (MQ=255)
                                              ttGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTc                                                      <  1:3713716/38‑1 (MQ=255)
                                              ttGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTc                                                      <  1:3043389/38‑1 (MQ=255)
                                                 cgcCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAAt                        <  1:3315980/65‑1 (MQ=255)
                                                  tccGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGt                           <  1:3359994/60‑1 (MQ=255)
                                                  cccGTTTTTGCGCCACGTCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGt                           <  1:1404750/60‑1 (MQ=255)
                                                   ccGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGcc                <  1:2143103/71‑1 (MQ=255)
                                                    cGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGgaga                                         >  1:3589201/1‑45 (MQ=255)
                                                    cGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGgaga                                         >  1:1677981/1‑45 (MQ=255)
                                                    cGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGgaga                                         >  1:3668976/1‑45 (MQ=255)
                                                    cGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGgaca                                         >  1:1407346/1‑43 (MQ=255)
                                                     gTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTg                                             <  1:516537/40‑1 (MQ=255)
                                                     gTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCaa              <  1:112012/71‑1 (MQ=255)
                                                      tttttGCGCCACGTCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAg             >  1:3359334/1‑71 (MQ=255)
                                                      tttttGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAg             >  1:157239/1‑71 (MQ=255)
                                                      tttttGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATAATCGGTAATTCCGCGCCAAg             >  1:2814656/1‑71 (MQ=255)
                                                        tttGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTg                                             >  1:2901085/1‑37 (MQ=255)
                                                        tttGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTg                                             >  1:2434413/1‑37 (MQ=255)
                                                        tttGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTg                                             >  1:112006/1‑37 (MQ=255)
                                                        tttGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAgcg           >  1:2582088/1‑71 (MQ=255)
                                                                 cGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCtt  >  1:1287402/1‑71 (MQ=255)
                                                                   |                                                                    
CAGTAAACTTCGCCCATTCGCGCGTCAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTT  >  minE/1194813‑1194948

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: