Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,186,871 A→G 24.2% Q263Q (CAA→CAG yeaH → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,186,8710AG24.2% 142.1 / 40.8 115Q263Q (CAA→CAGyeaHconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (45/42);  new base G (8/20);  total (53/62)
Fisher's exact test for biased strand distribution p-value = 4.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 115 total aligned reads displayed.

CCGATCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTCAGCAGAA  >  minE/1186810‑1186929
                                                             |                                                          
   aTCCCTCCAGCCAGGCAGTGATGTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:2011439/69‑1 (MQ=255)
    tCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:704243/69‑1 (MQ=255)
    tCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:3270374/69‑1 (MQ=255)
    tCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:3272623/69‑1 (MQ=255)
    tCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:3200551/69‑1 (MQ=255)
    tCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:1439268/69‑1 (MQ=255)
    tCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:2875813/69‑1 (MQ=255)
    tCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:616007/69‑1 (MQ=255)
    tCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:835989/69‑1 (MQ=255)
    tCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:1130728/69‑1 (MQ=255)
    tCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:205646/69‑1 (MQ=255)
    tCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:1038807/69‑1 (MQ=255)
    tCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:2058373/69‑1 (MQ=255)
    tCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:2069523/69‑1 (MQ=255)
    tCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:2097930/69‑1 (MQ=255)
     cccTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTTAAGa                                                 <  1:1711958/68‑1 (MQ=255)
     cccTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGat                                                >  1:2523008/1‑69 (MQ=255)
     cccTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGat                                                >  1:1213110/1‑69 (MQ=255)
     cccTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGat                                                >  1:3371741/1‑69 (MQ=255)
     cccTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:1512668/68‑1 (MQ=255)
     cccTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:2107322/68‑1 (MQ=255)
     cccTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:2659190/68‑1 (MQ=255)
     cccTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:31722/68‑1 (MQ=255)
     cccTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:2145143/68‑1 (MQ=255)
     cccTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:2564888/68‑1 (MQ=255)
     cccTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:1074939/68‑1 (MQ=255)
          cagccagGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATgg                                            >  1:2892010/1‑68 (MQ=255)
          cagccagGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATgg                                            >  1:2872774/1‑68 (MQ=255)
          cagccagGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATgg                                            >  1:912647/1‑68 (MQ=255)
          cagccagGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATgg                                            >  1:1022590/1‑68 (MQ=255)
          cagccagGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATgg                                            >  1:1012097/1‑68 (MQ=255)
          cagccagGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATgg                                            >  1:3286007/1‑68 (MQ=255)
          cagccagGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATgg                                            >  1:2114174/1‑68 (MQ=255)
          cagccagGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATgg                                            >  1:620716/1‑68 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATTGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:2652969/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:2434678/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:2605421/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:1006911/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:262730/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:31546/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:3367537/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:482580/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:554740/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:852122/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:928837/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:1159872/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:1439753/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:116394/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:118345/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:2161879/69‑1 (MQ=255)
                     tGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAGTCCACTAAAGATATGGCTAAACGttttt                                <  1:1783117/69‑1 (MQ=255)
                      gATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 >  1:1006497/1‑51 (MQ=255)
                      gATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 >  1:159564/1‑51 (MQ=255)
                      gATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 >  1:1461261/1‑51 (MQ=255)
                      gATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 >  1:3149041/1‑51 (MQ=255)
                      gATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 >  1:3104282/1‑51 (MQ=255)
                      gATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 >  1:3100237/1‑51 (MQ=255)
                      gATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 >  1:2892973/1‑51 (MQ=255)
                      gATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 >  1:53558/1‑51 (MQ=255)
                      gATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 >  1:130434/1‑51 (MQ=255)
                      gATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 >  1:2324937/1‑51 (MQ=255)
                      gATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACAAAAg                                                  >  1:2483790/1‑50 (MQ=255)
                                  gATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTa                   >  1:2699304/1‑69 (MQ=255)
                                  gATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTa                   >  1:1804603/1‑69 (MQ=255)
                                  gATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTa                   >  1:430843/1‑69 (MQ=255)
                                  gATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTa                   >  1:1351736/1‑69 (MQ=255)
                                  gATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTa                   >  1:847625/1‑69 (MQ=255)
                                  gATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTa                   >  1:35052/1‑69 (MQ=255)
                                  gATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTa                   >  1:3411337/1‑69 (MQ=255)
                                  gATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTa                   >  1:2453624/1‑69 (MQ=255)
                                      gACGTTTCCGGTTCAATGGATCAATCCACTATATATATGGCTAAGCGTTTTTATATTCTGCTGTATCTg               <  1:589540/69‑1 (MQ=255)
                                      gACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:534284/35‑1 (MQ=255)
                                      gACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:161624/35‑1 (MQ=255)
                                      gACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:2040270/35‑1 (MQ=255)
                                      gACGTTTCCGGTTCAATGGATCAATCCACTAAAGa                                                 <  1:2275380/35‑1 (MQ=255)
                                      gACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTg               <  1:1596306/69‑1 (MQ=255)
                                      gACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTg               <  1:766196/69‑1 (MQ=255)
                                      gACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTg               <  1:1969145/69‑1 (MQ=255)
                                      gACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTg               <  1:1394493/69‑1 (MQ=255)
                                      gACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTg               <  1:855632/69‑1 (MQ=255)
                                      gACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTg               <  1:928581/69‑1 (MQ=255)
                                        cGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTatat                            >  1:2514205/1‑54 (MQ=255)
                                        cGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTatat                            >  1:1951912/1‑54 (MQ=255)
                                        cGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTatat                            >  1:2454728/1‑54 (MQ=255)
                                                 ttCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTcagcag    <  1:2888686/69‑1 (MQ=255)
                                                 ttCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTcagcag    <  1:497725/69‑1 (MQ=255)
                                                 ttCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTcagcag    <  1:1221570/69‑1 (MQ=255)
                                                   cAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTCAGCAGaa  >  1:2491711/1‑69 (MQ=255)
                                                   cAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTCAGCAGaa  >  1:82292/1‑69 (MQ=255)
                                                   cAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTCAGCAGaa  >  1:2395605/1‑69 (MQ=255)
                                                   cAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTCAGCAGaa  >  1:736313/1‑69 (MQ=255)
                                                   cAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTCAGCAGaa  >  1:1234879/1‑69 (MQ=255)
                                                   cAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTCAGCAGaa  >  1:2618770/1‑69 (MQ=255)
                                                   cAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTCAGCAGaa  >  1:1908032/1‑69 (MQ=255)
                                                   cAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTCAGCAGaa  >  1:1713586/1‑69 (MQ=255)
                                                   cAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTCAGCAGaa  >  1:1713401/1‑69 (MQ=255)
                                                   cAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTCAGCAGaa  >  1:3391040/1‑69 (MQ=255)
                                                   cAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTCAGCAGaa  >  1:3183811/1‑69 (MQ=255)
                                                   cAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTCAGCAGaa  >  1:1699406/1‑69 (MQ=255)
                                                   cAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTCAGCAGaa  >  1:239583/1‑69 (MQ=255)
                                                             |                                                          
CCGATCCCTCCAGCCAGGCAGTGATGTTTTGCCTGATGGACGTTTCCGGTTCAATGGATCAATCCACTAAAGATATGGCTAAGCGTTTTTATATTCTGCTGTATCTGTTCCTCAGCAGAA  >  minE/1186810‑1186929

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: