Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,195,448 C→T 30.3% V568V (GTG→GTA yoaA ← conserved hypothetical protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,195,4480CT30.3% 88.4 / 42.1 92V568V (GTG→GTAyoaAconserved hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base T (13/15);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

CGTCGGCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCA  >  minE/1195384‑1195514
                                                                |                                                                  
cgtcgGCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACtt                                                                <  1:1028066/69‑1 (MQ=255)
cgtcgGCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACtt                                                                <  1:1079845/69‑1 (MQ=255)
cgtcgGCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACtt                                                                <  1:2419128/69‑1 (MQ=255)
          gcgAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCa                                                              >  1:1156197/1‑61 (MQ=255)
          gcgAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCCTCGCACgg                                                      >  1:2063533/1‑69 (MQ=255)
          gcgAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACgg                                                      >  1:971574/1‑69 (MQ=255)
          gcgAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACgg                                                      >  1:3149768/1‑69 (MQ=255)
          gcgAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACgg                                                      >  1:2945335/1‑69 (MQ=255)
          gcgAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACgg                                                      >  1:2330523/1‑69 (MQ=255)
          gcgAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACgg                                                      >  1:167126/1‑69 (MQ=255)
            gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                                     <  1:2063211/68‑2 (MQ=255)
            gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                                     <  1:2442992/68‑2 (MQ=255)
            gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                                     <  1:1954520/68‑2 (MQ=255)
            gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                                     <  1:2511617/68‑2 (MQ=255)
            gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                                     <  1:2719335/68‑2 (MQ=255)
            gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                                     <  1:273599/68‑2 (MQ=255)
            gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                                     <  1:3177260/68‑2 (MQ=255)
            gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                                     <  1:1081096/68‑2 (MQ=255)
            gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                                     <  1:1279117/68‑2 (MQ=255)
            gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                                     <  1:1521142/68‑2 (MQ=255)
            gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                                     <  1:1559398/68‑2 (MQ=255)
            gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                                     <  1:2795261/68‑2 (MQ=255)
            gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                                     <  1:2927747/68‑2 (MQ=255)
             aatcaatcGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                                     <  1:1550896/67‑2 (MQ=255)
                                  cttGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                                     <  1:2088320/45‑2 (MQ=25)
                                    tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                             >  1:1572431/1‑68 (MQ=255)
                                    tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                             >  1:1109759/1‑68 (MQ=255)
                                    tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                             >  1:2753097/1‑68 (MQ=255)
                                    tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                             >  1:3281235/1‑68 (MQ=255)
                                    tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                             >  1:1388854/1‑68 (MQ=255)
                                    tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                             >  1:2606907/1‑68 (MQ=255)
                                    tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                             >  1:3343862/1‑68 (MQ=255)
                                    tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                             >  1:369364/1‑68 (MQ=255)
                                    tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                             >  1:792803/1‑68 (MQ=255)
                                    tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                             >  1:2046670/1‑68 (MQ=255)
                                    tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                             >  1:2279381/1‑68 (MQ=255)
                                    tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                             >  1:1309215/1‑68 (MQ=255)
                                      agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGTTCGCCACCGCGCAAACGACAATCTTc                           <  1:1027920/68‑1 (MQ=255)
                                      agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                           <  1:3168337/68‑1 (MQ=255)
                                      agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                           <  1:1497078/68‑1 (MQ=255)
                                      agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                           <  1:3299977/68‑1 (MQ=255)
                                      agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                           <  1:1205203/68‑1 (MQ=255)
                                      agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                           <  1:1737149/68‑1 (MQ=255)
                                              gACGGCATCTGGTAGTTGCACTTCATCTAACGGGTCGCCACCGCGCAAACGACAATCTTCCATgcgcg                   <  1:1240699/68‑1 (MQ=255)
                                              gACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATgcgcg                   <  1:592706/68‑1 (MQ=255)
                                              gACGGCATCTGGTAGTTGCACTTCATCGAACGGGGCGCCACCGCGCAAACGACAATCTTc                           <  1:1649760/60‑1 (MQ=255)
                                                  gCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCtttt              <  1:2651409/69‑1 (MQ=255)
                                                  gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTcc                          >  1:2750770/1‑57 (MQ=255)
                                                  gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt              >  1:1766677/1‑69 (MQ=255)
                                                  gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt              >  1:3059343/1‑69 (MQ=255)
                                                  gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt              >  1:2940174/1‑69 (MQ=255)
                                                  gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt              >  1:1385147/1‑69 (MQ=255)
                                                  gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt              >  1:3382154/1‑69 (MQ=255)
                                                  gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt              >  1:369199/1‑69 (MQ=255)
                                                  gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt              >  1:1372538/1‑69 (MQ=255)
                                                  gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt              >  1:2619957/1‑69 (MQ=255)
                                                  gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt              >  1:371337/1‑69 (MQ=255)
                                                  gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt              >  1:420987/1‑69 (MQ=255)
                                                  gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt              >  1:765249/1‑69 (MQ=255)
                                                  gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt              >  1:646772/1‑69 (MQ=255)
                                                     tCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAAc           <  1:3023428/69‑1 (MQ=255)
                                                       tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACaa                                >  1:297784/1‑46 (MQ=255)
                                                       tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                           <  1:3388037/51‑1 (MQ=255)
                                                       tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                           <  1:221187/51‑1 (MQ=255)
                                                       tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                           <  1:2266175/51‑1 (MQ=255)
                                                       tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                           <  1:2777270/51‑1 (MQ=255)
                                                       tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         >  1:1580254/1‑69 (MQ=255)
                                                       tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         >  1:1848657/1‑69 (MQ=255)
                                                       tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         >  1:150283/1‑69 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         <  1:426174/62‑1 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         <  1:330389/62‑1 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         <  1:1321120/62‑1 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         <  1:1516221/62‑1 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         <  1:1289505/62‑1 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         <  1:3010453/62‑1 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         <  1:2938959/62‑1 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         <  1:2898657/62‑1 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         <  1:1795558/62‑1 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         <  1:2716879/62‑1 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         <  1:2667696/62‑1 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         <  1:3148908/62‑1 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         <  1:258581/62‑1 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         <  1:2017928/62‑1 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg         <  1:2300593/62‑1 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGatca  >  1:1840823/1‑69 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGatca  >  1:2156353/1‑69 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGatca  >  1:2058449/1‑69 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGatca  >  1:2441812/1‑69 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGatca  >  1:429907/1‑69 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGatca  >  1:295028/1‑69 (MQ=255)
                                                              tGCACTTCATCGAACGGGTCGCCACCGCGCAAAAGACAATCTTCCATGCGCGCTTTTAACAGTGGatca  >  1:1423318/1‑69 (MQ=255)
                                                               gCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATgc                      >  1:1195563/1‑48 (MQ=255)
                                                                |                                                                  
CGTCGGCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCA  >  minE/1195384‑1195514

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: