Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,117,769 C→T 38.2% K218K (AAG→AAA yraL ← predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,117,7690CT38.2% 32.6 / 36.3 60K218K (AAG→AAAyraLpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (5/32);  new base T (7/16);  total (12/48)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

TGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAG  >  minE/2117710‑2117831
                                                           |                                                              
tgtgCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTAcc                                                          <  1:2842798/66‑1 (MQ=255)
tgtgCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTAcc                                                          <  1:702984/66‑1 (MQ=255)
 gtgCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCAc                                                       <  1:2959514/68‑1 (MQ=255)
  tgCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTAcc                                                          <  1:554223/64‑1 (MQ=255)
    cTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGc                                                     <  1:2532129/67‑1 (MQ=255)
    cTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGc                                                     <  1:95854/67‑1 (MQ=255)
          gACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCagcagc                                              <  1:375520/68‑1 (MQ=255)
          gACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCagcagc                                              <  1:2419477/68‑1 (MQ=255)
          gACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCagcagc                                              <  1:2138174/68‑1 (MQ=255)
          gACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCagcagc                                              <  1:2137635/68‑1 (MQ=255)
          gACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCagcagc                                              <  1:2124668/68‑1 (MQ=255)
          gACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCagcagc                                              <  1:3184271/68‑1 (MQ=255)
          gACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCagcagc                                              <  1:3203401/68‑1 (MQ=255)
          gACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCagcagc                                              <  1:3265722/68‑1 (MQ=255)
          gACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCagcagc                                              <  1:45833/68‑1 (MQ=255)
           aCCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTc                                            >  1:1896421/1‑69 (MQ=255)
              ttCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCTTTTACCCACGCCAGCAGCTCGcc                                         >  1:1650597/1‑69 (MQ=255)
              ttCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCTTTTACCCACGCCAGCAGCTCGcc                                         >  1:1538776/1‑69 (MQ=255)
              ttCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCTTTTACCCACGCCAGCAGCTCGcc                                         >  1:1440300/1‑69 (MQ=255)
              ttCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCTTTTACCCACGCCAGCAGCTCGcc                                         >  1:1410010/1‑69 (MQ=255)
              ttCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCTTTTACCCACGCCAGCAGCTCGcc                                         >  1:2677673/1‑69 (MQ=255)
              ttCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCTTTTACCCACGCCAGCAGCTCGcc                                         >  1:2614995/1‑69 (MQ=255)
              ttCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCTTTTACCCACGCCAGCAGCTCGc                                          >  1:1616722/1‑68 (MQ=255)
                     aTCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACggg                                   >  1:3063467/1‑68 (MQ=255)
                       cAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGgcgc                                >  1:535230/1‑69 (MQ=255)
                       cAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGgcgc                                >  1:1024286/1‑69 (MQ=255)
                       cAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGgcg                                 >  1:950510/1‑68 (MQ=255)
                        aGCACCATTTCGCCTTTGCGACTGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGcc                                         <  1:123374/59‑1 (MQ=255)
                        aGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCATCTCGcc                                         <  1:988964/59‑1 (MQ=255)
                        aGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGcc                                         <  1:2322525/59‑1 (MQ=255)
                        aGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGcc                                         <  1:922399/59‑1 (MQ=255)
                        aGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGcc                                         <  1:2244720/59‑1 (MQ=255)
                        aGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGcc                                         <  1:2012483/59‑1 (MQ=255)
                        aGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGcc                                         <  1:905560/59‑1 (MQ=255)
                        aGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGcc                                         <  1:3046149/59‑1 (MQ=255)
                        aGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGcc                                         <  1:724284/59‑1 (MQ=255)
                        aGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGcc                                         <  1:1639481/59‑1 (MQ=255)
                        aGCACCATTTCGCCTTAGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGcc                                         <  1:1729891/59‑1 (MQ=255)
                         gCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGcc                                         <  1:1764554/58‑1 (MQ=255)
                           aCCATCTCGCCTTTGCTAGGGTTTTCATCTTCTTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:826162/69‑1 (MQ=255)
                           aCCATCTCGCCTTTGCTAGGGTTTTCATCTTCTTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:343799/69‑1 (MQ=255)
                           aCCATCTCGCCTTTGCTAGGGTTTTCATCTTCTTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:1677852/69‑1 (MQ=255)
                           aCCATCTCGCCTTTGCTAGGGTTTTCATCTTCTTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:1621926/69‑1 (MQ=255)
                           aCCATCTCGCCTTTGCTAGGGTTTTCATCTTCTTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:1488913/69‑1 (MQ=255)
                           aCCATCTCGCCTTTGCTAGGGTTTTCATCTTCTTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:1021509/69‑1 (MQ=255)
                           aCCATCTCGCCTTTGCTAGGGTTTTCATCTTCTTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:1298352/69‑1 (MQ=255)
                           aCCATCTCGCCTTTGCTAGGGTTTTCATCTTCTTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:2516590/69‑1 (MQ=255)
                           aCCATCTCGCCTTTGCTAGGGTTTTCATCTTCTTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:3331439/69‑1 (MQ=255)
                           aCCATCTCGCCTTTGCTAGGGTTTTCATCTTCTTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:3185971/69‑1 (MQ=255)
                           aCCATCTCGCCTTTGCTAGGGTTTTCATCTTCTTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:2821540/69‑1 (MQ=255)
                           aCCATCTCGCCTTTGCTAGGGTTTTCATCTTCTTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:2003196/69‑1 (MQ=255)
                           aCCATCTCGCCTTTGCTAGGGTTTTCATCTTATTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:397237/69‑1 (MQ=255)
                            ccATCTCGCCTTTGCTAGGGTTTTCATCTTCTTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:455454/68‑1 (MQ=255)
                            ccATCTCGCCTTTGCTAGGGTTTTCATCTTCTTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:3096337/68‑1 (MQ=255)
                             cATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTg                          <  1:343785/69‑1 (MQ=255)
                             cATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTg                          <  1:2444317/69‑1 (MQ=255)
                             cATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTg                          <  1:137833/69‑1 (MQ=255)
                             cATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTg                          <  1:2977454/69‑1 (MQ=255)
                                 tCGCCTTTGCTAGGGTTTTCATCTTCTTTGATCCACGCCAGCAGTTCGCCGACCGACGTGCCg                            <  1:60605/63‑1 (MQ=255)
                                     cTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGtt                   <  1:1969299/68‑1 (MQ=255)
                                                      tcttctTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTTAg  <  1:2473716/68‑1 (MQ=255)
                                                      tcttctTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTTAg  <  1:852710/68‑1 (MQ=255)
                                                      tcttctTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTTAg  <  1:767572/68‑1 (MQ=255)
                                                      tcttctTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTTAg  <  1:677633/68‑1 (MQ=255)
                                                      tcttctTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTTAg  <  1:1969373/68‑1 (MQ=255)
                                                      tcttctTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCAGTGAATGGTTTCCCAGGTTTTGGTTAg  <  1:3173322/68‑1 (MQ=255)
                                                      tcttctTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCAGTGAATGGTTTCCCAGGTTTTGGTTAg  <  1:650221/68‑1 (MQ=255)
                                                           |                                                              
TGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAG  >  minE/2117710‑2117831

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: