Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,876,560 A→G 44.0% S288S (TCT→TCC eno ← enolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,876,5600AG44.0% 13.7 / 59.3 75S288S (TCT→TCCenoenolase
Reads supporting (aligned to +/- strand):  ref base A (13/29);  new base G (15/18);  total (28/47)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCCCAGAACTTTGGTCTGGTATGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGAATTCTTCAGAGGTGAACGCTTT  >  minE/1876495‑1876628
                                                                 |                                                                    
tCGGCCAGAACTTTGGTCTGGTATGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATAGa                                                                    <  1:2939390/68‑1 (MQ=255)
tCGCCCAGAACTTTGGTCTGGTATGCGAAACCGTCCCAGTCAGATTCGTCCAGGCCGTCTTCGATAGa                                                                    <  1:2012573/68‑1 (MQ=255)
tCGCCCAGAACTTTGGTCTGGTATGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATAGa                                                                    <  1:1047383/68‑1 (MQ=255)
tCGCCCAGAACTTTGGTCTGGTATGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATAGa                                                                    <  1:2462590/68‑1 (MQ=255)
tCGCCCAGAACTTTGGTCTGGTATGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATAGa                                                                    <  1:2753824/68‑1 (MQ=255)
tCGCCCAGAACTTTGGTCTGGTATGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATAGa                                                                    <  1:528408/68‑1 (MQ=255)
       gAACTTTGGTCTGGTATGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCg                                                            <  1:2582547/69‑1 (MQ=255)
         aCTTTGGTCTGGTATGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCggg                                                          <  1:1286209/69‑1 (MQ=255)
         aCTTTGGTCTGGTATGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCggg                                                          <  1:2448634/69‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCTGGTACTGTTTgg                                                <  1:1531761/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:958992/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:2177034/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:2207255/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:220946/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:2289036/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:936573/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:2523804/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:809440/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:2795096/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:760257/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:319760/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:320179/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:456760/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:724818/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:1071322/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTgg                                                <  1:1536985/66‑1 (MQ=255)
                      aTGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGGACTGTTTgg                                                <  1:1408914/66‑1 (MQ=255)
                       tGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAg                                            <  1:3066420/69‑1 (MQ=255)
                       tGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAg                                            <  1:424853/69‑1 (MQ=255)
                        gCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAg                                            <  1:124502/68‑1 (MQ=255)
                                gTCCCAGTCAGATTCTTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGtt                                                   >  1:3037148/1‑53 (MQ=255)
                                gTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCa                                    >  1:1582961/1‑68 (MQ=255)
                                gTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCa                                    >  1:3151500/1‑68 (MQ=255)
                                gTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCa                                    >  1:3236309/1‑68 (MQ=255)
                                gTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCa                                    >  1:2952434/1‑68 (MQ=255)
                                gTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCa                                    >  1:418227/1‑68 (MQ=255)
                                gTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCa                                    >  1:438788/1‑68 (MQ=255)
                                gTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCa                                    >  1:2610507/1‑68 (MQ=255)
                                gTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCa                                    >  1:457217/1‑68 (MQ=255)
                                gTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCa                                    >  1:2466464/1‑68 (MQ=255)
                                gTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCa                                    >  1:2462612/1‑68 (MQ=255)
                                gTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCa                                    >  1:559312/1‑68 (MQ=255)
                                gTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCa                                    >  1:1079869/1‑68 (MQ=255)
                                gTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCa                                    >  1:2143369/1‑68 (MQ=255)
                                gTCCCAGTCAGATTCGTCCAGACCGTCTTCGATGCAAACGCTCGGGTACTGtt                                                   >  1:1085958/1‑53 (MQ=255)
                                   ccAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGaa                                >  1:2766231/1‑69 (MQ=255)
                                   ccAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGaa                                >  1:1296091/1‑69 (MQ=255)
                                   ccAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGaa                                >  1:2823640/1‑69 (MQ=255)
                                   ccAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGaa                                >  1:1786621/1‑69 (MQ=255)
                                   ccAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGaa                                >  1:1790579/1‑69 (MQ=255)
                                   ccAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGaa                                >  1:2151256/1‑69 (MQ=255)
                                   ccAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGaa                                >  1:3075424/1‑69 (MQ=255)
                                   ccAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGa                                 >  1:3079833/1‑68 (MQ=255)
                                                  cAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAGt                                           >  1:2036137/1‑43 (MQ=255)
                                                  cAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAagt                              >  1:213033/1‑56 (MQ=255)
                                                  cAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAagt                              >  1:2371561/1‑56 (MQ=255)
                                                  cAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAagt                              >  1:2913831/1‑56 (MQ=255)
                                                     aCCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGAATTCTTCaga              >  1:143900/1‑69 (MQ=255)
                                                           ttCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGaa                       <  1:2372558/54‑1 (MQ=255)
                                                           ttCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGaa                       <  1:1786813/54‑1 (MQ=255)
                                                           ttCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGaa                       <  1:780710/54‑1 (MQ=255)
                                                           ttCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGaa                       <  1:641490/54‑1 (MQ=255)
                                                           ttCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGaa                       <  1:3059005/54‑1 (MQ=255)
                                                           ttCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGaa                       <  1:1540903/54‑1 (MQ=255)
                                                           ttCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGaa                       <  1:1507413/54‑1 (MQ=255)
                                                           ttCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGaa                       <  1:41629/54‑1 (MQ=255)
                                                           ttCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGaa                       <  1:3277289/54‑1 (MQ=255)
                                                           ttCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGaa                       <  1:2974210/54‑1 (MQ=255)
                                                            tCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGaa                       <  1:1488912/53‑1 (MQ=255)
                                                               aTAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGAATTCTTCAGAGGTGAACGc     <  1:1684456/68‑1 (MQ=255)
                                                               aTAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGAATTCTTCAGAGGTGAACGc     <  1:1257857/68‑1 (MQ=255)
                                                               aTAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGAATTCTTCAGAGGTGAACGc     <  1:736590/68‑1 (MQ=255)
                                                               aTAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGAATTCTTCAGAGGTGAACGc     <  1:914088/68‑1 (MQ=255)
                                                               aTAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGAATTCTTCAGAGGTGAACGc     <  1:1939060/68‑1 (MQ=255)
                                                                tAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGAATTCTTCAGAGGTGAACGc     <  1:1382426/67‑1 (MQ=255)
                                                                 aGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGAATTCTTCAGAGGTGAACGCttt  <  1:2953443/69‑1 (MQ=255)
                                                                 aGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGAATTCTTCAGAGGTGAACGCttt  <  1:834628/69‑1 (MQ=255)
                                                                 aGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGAATTCTTCAGAGGTGAACGCttt  <  1:864323/69‑1 (MQ=255)
                                                                 aGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGAATTCTTCAGAGGTGAACGCttt  <  1:1947499/69‑1 (MQ=255)
                                                                 |                                                                    
TCGCCCAGAACTTTGGTCTGGTATGCGAAACCGTCCCAGTCAGATTCGTCCAGACCGTCTTCGATAGAAACGATCGGGTACTGTTTGGTCAGTTCTTCCAGGAAGTGAGTGAATTCTTCAGAGGTGAACGCTTT  >  minE/1876495‑1876628

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: