Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,381,717 A→G 60.7% L94L (TTG→CTG)  yeiB ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,381,7170AG60.7% 31.8 / 41.8 66L94L (TTG→CTG) yeiBconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (16/10);  new base G (33/7);  total (49/17)
Fisher's exact test for biased strand distribution p-value = 8.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCACCAGCCCGTAAGCCAGCAAAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAACAAGACTAACAGCGTTAACCGCGACTGGATCCAGCGTCTGCCACGGGGCAGCAACATTTGC  >  minE/1381653‑1381778
                                                                |                                                             
cccACCAGCCCGTAAGCCAGCAAAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCcaacaa                                                            <  1:1238482/68‑1 (MQ=255)
cccACCAGCCCGTAAGCCAGCAAAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCcaacaa                                                            <  1:2854311/68‑1 (MQ=255)
cccACCAGCCCGTAAGCCAGCAAAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAACAAg                                                           <  1:3406825/69‑1 (MQ=255)
cccACCAGCCCGTAAGCCAGCAAAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAACAAg                                                           <  1:2131435/69‑1 (MQ=255)
cccACCAGCCCGTAAGCCAGCAAAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAACAAg                                                           <  1:2771183/69‑1 (MQ=255)
cccACCAGCCCGTAAGCCAGCAAAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAACAAg                                                           <  1:3372201/69‑1 (MQ=255)
cccACCAGCCCGTAAGCCAGCAAAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAACAAg                                                           <  1:680495/69‑1 (MQ=255)
   accaGCCCGTAAGCCAGCAAAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAACAAGACt                                                        <  1:2354007/69‑1 (MQ=255)
   accaGCCCGTAAGCCAGCAAAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAACAAGACt                                                        <  1:617227/69‑1 (MQ=255)
    ccaGCCCGTAAGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACt                                                        >  1:1631546/1‑68 (MQ=255)
    ccaGCCCGTAAGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACt                                                        >  1:1149027/1‑68 (MQ=255)
    ccaGCCCGTAAGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACt                                                        >  1:1828850/1‑68 (MQ=255)
    ccaGCCCGTAAGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACt                                                        >  1:2255182/1‑68 (MQ=255)
    ccaGCCCGTAAGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACt                                                        >  1:239462/1‑68 (MQ=255)
    ccaGCCCGTAAGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACt                                                        >  1:2515294/1‑68 (MQ=255)
    ccaGCCCGTAAGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACt                                                        >  1:2526983/1‑68 (MQ=255)
    ccaGCCCGTAAGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACt                                                        >  1:797832/1‑68 (MQ=255)
    ccaGCCCGTAAGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACt                                                        >  1:785878/1‑68 (MQ=255)
    ccaGCCCGTAAGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACt                                                        >  1:290880/1‑68 (MQ=255)
    ccaGCCCGTAAGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACt                                                        >  1:134007/1‑68 (MQ=255)
    ccaGCCCGTAAGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACt                                                        >  1:1481970/1‑68 (MQ=255)
    ccaGCCCGTAAGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACt                                                        >  1:1381123/1‑68 (MQ=255)
    ccaGCCCGTAAGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAAAAAAGCCCAGCAAGACt                                                        >  1:672072/1‑68 (MQ=255)
             aaGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCTTGAATAAAGCCCAGCAAGACTAACAGCGTTa                                              <  1:1428743/69‑1 (MQ=255)
             aaGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACTAACAGCGTTa                                              <  1:1175429/69‑1 (MQ=255)
             aaGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACTAACAGCGTTa                                              <  1:1857597/69‑1 (MQ=255)
             aaGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACTAACAGCGTTa                                              <  1:1572003/69‑1 (MQ=255)
             aaGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACTAACAGCGTTa                                              <  1:1840454/69‑1 (MQ=255)
             aaGCCAGCAGAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAGCAAGACTAACAGCGTTa                                              <  1:164796/69‑1 (MQ=255)
                     aaaaTATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAACAAGACTAACAGCGTTAACCGCGAc                                      <  1:2958829/69‑1 (MQ=255)
                     aaaaTATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAAAAAGACTAACAGCGTTAACCGCGAc                                      <  1:1659355/69‑1 (MQ=255)
                        atatCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCcaacaa                                                            >  1:2744114/1‑44 (MQ=255)
                        atatCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCcaacaa                                                            >  1:3196078/1‑44 (MQ=255)
                        atatCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCcaacaa                                                            >  1:2865062/1‑44 (MQ=255)
                        atatCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCcaacaa                                                            >  1:547202/1‑44 (MQ=255)
                        atatCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCcaacaa                                                            >  1:572082/1‑44 (MQ=255)
                        atatCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCcaacaa                                                            >  1:2360273/1‑44 (MQ=255)
                        atatCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCcaacaa                                                            >  1:1397010/1‑44 (MQ=255)
                        atatCGCCGTCCCAGAACAATAAACCGTGAATAAAACCCAACAAGACTAACAGCGTTAACCGCGACTgg                                   >  1:148622/1‑69 (MQ=255)
                        atatCGCCGTCCCAGAACAATAAACCGTGAATAAAACCCAACAAGACTAACAGCGTTAACCGCGACTgg                                   >  1:36505/1‑69 (MQ=255)
                        atatCGCCGTCCCAGAACAATAAACCGTGAATAAAACCCAACAAGACTAACAGCGTTAACCGCGACTgg                                   >  1:1691711/1‑69 (MQ=255)
                        atatCGCCGTCCCAGAACAATAAACCGTGAATAAAACCCAACAAGACTAACAGCGTTAACCGCGACTgg                                   >  1:1296172/1‑69 (MQ=255)
                        atatCGCCGTCCCAGAACAATAAACCGTGAATAAAACCCAACAAGACTAACAGCGTTAACCGCGACTgg                                   >  1:2950760/1‑69 (MQ=255)
                        atatCGCCGTACCAGAACAATAAACCGTGAATAAAACCCAACAAGACTAACAGCGTTAACCGCGACTgg                                   >  1:2954671/1‑69 (MQ=255)
                         tatCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCTGCAAGACTAACAGCGTTa                                              <  1:2797595/57‑1 (MQ=255)
                               cGTCCCAGAACAATAAACCGTGAATAAAGCCCAACAAGACTAACAGCGTTAACCGCGACTGGATCCAGc                            <  1:3053499/69‑1 (MQ=255)
                                        aCAATAAACCGTGAATAAAGCCCAACAAGACTAACAGCGTTAACCGCGACTGGATCCAGCGTCTGCCAc                   >  1:2384243/1‑69 (MQ=255)
                                        aCAATAAACCGTGAATAAAGCCCAACAAGACTAACAGCGTTAACCGCGACTGGATCCAGCGTCTGCCAc                   >  1:2598658/1‑69 (MQ=255)
                                         cAATAAACCGTGAATAAAGCCCAACAAGACTAACAGCGTTAACCGCGACTGGATCCAGCGTCTGCCACg                  >  1:3019563/1‑69 (MQ=255)
                                         cAATAAACCGTGAATAAAGCCCAACAAGACTAACAGCGTTAACCGCGACTGGATCCAGCGTCTGCCACg                  >  1:1785503/1‑69 (MQ=255)
                                         cAATAAACCGTGAATAAAGCCCAACAAGACTAACAGCGTTAACCGCGACTGGATCCAGCGTCTGCCACg                  >  1:1663686/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATc                                >  1:3429320/1‑46 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:1283483/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:1542241/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:350110/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:74145/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:1354791/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:426746/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:915816/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:3319865/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:3278849/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:3036316/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:3017932/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:272546/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:2532414/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:16710/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:1779160/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:2124205/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCACGGGGCAGCaa         >  1:1932986/1‑69 (MQ=255)
                                                  gTGAATAAAGCCCAGCAAGACTAACAGCGTTAACCGCGACTGGATCCATCGTCTGCCAAGGGGCAGCaa         >  1:2855085/1‑69 (MQ=255)
                                                         aaGCCCAACAAGACTAACAGCGTTTACCGCGACTGGATCCAGCGTCTGCCACGGGGCAGCAACATTTGc  >  1:2058723/1‑69 (MQ=255)
                                                         aaGCCCAACAAGACTAACAGCGTTAACCGCGACTGGATCCAGCGTCTGCCACGGGGCAGCAACATTTg   >  1:156589/1‑68 (MQ=255)
                                                         aaGCCCAACAAGACTAACAGCGTTAACCGCGACTGGATCCAGCGTCTGCCACGGGGCAGCAACATTTGc  >  1:2244887/1‑69 (MQ=255)
                                                         aaGCCCAACAAGACTAACAGCGTTAACCGCGACTGGATCCAGCGTCTGCCACGGGGCAGCAACATTTGc  >  1:718143/1‑69 (MQ=255)
                                                         aaGCCCAACAAGACTAACAGCGTTAACCGCGACTGGATCCAGCGTCTGCCACGGGGCAGCAACATTTGc  >  1:2120134/1‑69 (MQ=255)
                                                                |                                                             
CCCACCAGCCCGTAAGCCAGCAAAATATCGCCGTCCCAGAACAATAAACCGTGAATAAAGCCCAACAAGACTAACAGCGTTAACCGCGACTGGATCCAGCGTCTGCCACGGGGCAGCAACATTTGC  >  minE/1381653‑1381778

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: