Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,695,639 T→C 46.4% E124G (GAG→GGG)  yfhQ ← predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,695,6390TC46.4% 9.2 / 97.3 69E124G (GAG→GGG) yfhQpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (22/15);  new base C (20/12);  total (42/27)
Fisher's exact test for biased strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CATCGCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGC  >  minE/1695569‑1695706
                                                                      |                                                                   
cATCGCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAAct                                                                     <  1:810517/71‑1 (MQ=255)
cATCGCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAAct                                                                     <  1:2042335/71‑1 (MQ=255)
     ccAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCa                                                                <  1:2094863/71‑1 (MQ=255)
     ccAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCa                                                                <  1:743121/71‑1 (MQ=255)
     ccAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCa                                                                <  1:3080958/71‑1 (MQ=255)
      cAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCa                                                                <  1:531429/70‑1 (MQ=255)
         gTTCAGCGCGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTg                                                             >  1:2480208/1‑70 (MQ=255)
          ttCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGt                                                           <  1:227422/71‑1 (MQ=255)
            cAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTc                                                          <  1:2110550/70‑1 (MQ=255)
             aGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCa                                                         <  1:1471628/70‑1 (MQ=255)
               cGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGc                                                      >  1:1045485/1‑71 (MQ=255)
                gAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGcc                                                     <  1:2748520/71‑1 (MQ=255)
                 aGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGcc                                                     >  1:3625939/1‑70 (MQ=255)
                 aGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCg                                                    >  1:2989475/1‑71 (MQ=255)
                         tCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTc                                            >  1:275283/1‑71 (MQ=255)
                            gggTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTcgcg                                         <  1:3148538/71‑1 (MQ=255)
                               ttAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTcgcg                                         >  1:3493373/1‑68 (MQ=255)
                               ttAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGa                                        >  1:182416/1‑69 (MQ=255)
                               ttAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGa                                        >  1:851410/1‑69 (MQ=255)
                               ttAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGAcc                                      <  1:317713/71‑1 (MQ=255)
                                 aGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTc                                            >  1:1647038/1‑63 (MQ=255)
                                     gCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAacac                                >  1:1968475/1‑71 (MQ=255)
                                      cAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGAGCTCGCGACCAAacac                                >  1:2831237/1‑70 (MQ=255)
                                       aaTCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGAc                                                  <  1:3409430/51‑1 (MQ=255)
                                       aaTCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGAc                                                  <  1:2226633/51‑1 (MQ=255)
                                       aaTCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGAc                                                  <  1:3436616/51‑1 (MQ=255)
                                       aaTCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGAc                                                  <  1:2959151/51‑1 (MQ=255)
                                       aaTCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAg                                                       >  1:2140501/1‑46 (MQ=255)
                                       aaTCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAg                                                       >  1:1258708/1‑46 (MQ=255)
                                       aaTCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAg                                                       >  1:1024722/1‑46 (MQ=255)
                                       aaTCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAg                                                       >  1:1971758/1‑46 (MQ=255)
                                       aaTCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAg                                                       >  1:925438/1‑46 (MQ=255)
                                       aaTCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAg                                                       >  1:627162/1‑46 (MQ=255)
                                       aaTCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAg                                                       >  1:1427290/1‑46 (MQ=255)
                                          cgcgACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTcgcg                                         <  1:709696/57‑1 (MQ=255)
                                           gcgACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAgcgc                          >  1:705913/1‑71 (MQ=255)
                                           gcgACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAgcgc                          >  1:1886864/1‑71 (MQ=255)
                                           gcgACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAgcgc                          >  1:3436664/1‑71 (MQ=255)
                                           gcgACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAgcgc                          >  1:176674/1‑71 (MQ=255)
                                           gcgACATGAAAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAgcg                           <  1:26004/70‑1 (MQ=255)
                                                 tGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGa                                                   >  1:399337/1‑40 (MQ=255)
                                                 tGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGAc                                                  >  1:623328/1‑41 (MQ=255)
                                                 tGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGAc                                                  >  1:1210703/1‑41 (MQ=255)
                                                 tGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGAc                                                  >  1:2907832/1‑41 (MQ=255)
                                                 tGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGAc                                                  >  1:1987026/1‑41 (MQ=255)
                                                 tGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGa                                        >  1:905891/1‑51 (MQ=255)
                                                   aTAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAaca                                 >  1:1774251/1‑56 (MQ=255)
                                                   aTAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAaca                                 >  1:1756193/1‑56 (MQ=255)
                                                   aTAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAaca                                 >  1:2240129/1‑56 (MQ=255)
                                                   aTAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAaca                                 >  1:3013912/1‑56 (MQ=255)
                                                   aTAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAaca                                 >  1:2036309/1‑56 (MQ=255)
                                                   aTAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAac                                  >  1:3518973/1‑55 (MQ=255)
                                                   aTAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAgc                            >  1:313746/1‑61 (MQ=255)
                                                   aTAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACAcc                               <  1:1540999/58‑1 (MQ=255)
                                                    tAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGAc                                                  <  1:33131/38‑1 (MQ=255)
                                                    tAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACAcc                               <  1:2339373/57‑1 (MQ=255)
                                                      aTGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTAt               <  1:2746465/71‑1 (MQ=255)
                                                      aTGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTAt               <  1:1799695/71‑1 (MQ=255)
                                                      aTGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTAt               <  1:3078485/71‑1 (MQ=255)
                                                       tGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTAtt              <  1:362530/71‑1 (MQ=255)
                                                          cATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGa                                        >  1:2179767/1‑42 (MQ=255)
                                                          cATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGa                                        >  1:542446/1‑42 (MQ=255)
                                                          cATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGa                                        >  1:432334/1‑42 (MQ=255)
                                                            tttCTGCAACCCTTCATTGGTCAAGCCGACGCGCTcgcg                                         <  1:1276013/39‑1 (MQ=255)
                                                            tttCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCa                                     <  1:681890/43‑1 (MQ=255)
                                                            tttCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGc           >  1:2414561/1‑69 (MQ=255)
                                                            tttCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCg         >  1:158328/1‑71 (MQ=255)
                                                            tttCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCg         <  1:2902171/71‑1 (MQ=255)
                                                             ttCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGc           >  1:3483928/1‑68 (MQ=255)
                                                                 gCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAacac                                >  1:3452407/1‑43 (MQ=255)
                                                                   aaCCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTCCCGCTTCAGc  <  1:3316197/71‑1 (MQ=255)
                                                                      |                                                                   
CATCGCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGC  >  minE/1695569‑1695706

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: