Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,195,481 A→G 24.4% C557C (TGT→TGC 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,4810AG24.4% 127.7 / 37.7 106C557C (TGT→TGCyoaAconserved hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base A (43/37);  new base G (13/13);  total (56/50)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

Only 100 of 106 total aligned reads displayed.

TGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGAT  >  minE/1195420‑1195543
                                                             |                                                              
tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                                                          >  1:369364/1‑68 (MQ=255)
tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                                                          >  1:1109759/1‑68 (MQ=255)
tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                                                          >  1:3343862/1‑68 (MQ=255)
tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                                                          >  1:3281235/1‑68 (MQ=255)
tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                                                          >  1:2753097/1‑68 (MQ=255)
tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                                                          >  1:2606907/1‑68 (MQ=255)
tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                                                          >  1:2279381/1‑68 (MQ=255)
tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                                                          >  1:2046670/1‑68 (MQ=255)
tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                                                          >  1:1309215/1‑68 (MQ=255)
tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                                                          >  1:1572431/1‑68 (MQ=255)
tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                                                          >  1:1388854/1‑68 (MQ=255)
  agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGTTCGCCACCGCGCAAACGACAATCTTc                                                        <  1:1027920/68‑1 (MQ=255)
  agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                                                        <  1:1205203/68‑1 (MQ=255)
  agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                                                        <  1:1737149/68‑1 (MQ=255)
  agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                                                        <  1:3168337/68‑1 (MQ=255)
  agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                                                        <  1:3299977/68‑1 (MQ=255)
  agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                                                        <  1:1497078/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:371337/1‑69 (MQ=255)
              gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt                                           >  1:369199/1‑69 (MQ=255)
              gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt                                           >  1:3382154/1‑69 (MQ=255)
              gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt                                           >  1:420987/1‑69 (MQ=255)
              gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt                                           >  1:3059343/1‑69 (MQ=255)
              gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt                                           >  1:646772/1‑69 (MQ=255)
              gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt                                           >  1:2940174/1‑69 (MQ=255)
              gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt                                           >  1:765249/1‑69 (MQ=255)
              gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt                                           >  1:2619957/1‑69 (MQ=255)
              gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt                                           >  1:1372538/1‑69 (MQ=255)
              gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt                                           >  1:1385147/1‑69 (MQ=255)
              gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt                                           >  1:1766677/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:330389/62‑1 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg                                      <  1:2716879/62‑1 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg                                      <  1:2667696/62‑1 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg                                      <  1:1795558/62‑1 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg                                      <  1:1321120/62‑1 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg                                      <  1:1516221/62‑1 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg                                      <  1:2898657/62‑1 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg                                      <  1:2300593/62‑1 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg                                      <  1:426174/62‑1 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg                                      <  1:2938959/62‑1 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg                                      <  1:3148908/62‑1 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg                                      <  1:1289505/62‑1 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGatca                               >  1:295028/1‑69 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGatca                               >  1:1840823/1‑69 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGatca                               >  1:2441812/1‑69 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGatca                               >  1:429907/1‑69 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGatca                               >  1:2156353/1‑69 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGatca                               >  1:2058449/1‑69 (MQ=255)
                          tGCACTTCATCGAACGGGTCGCCACCGCGCAAAAGACAATCTTCCATGCGCGCTTTTAACAGTGGatca                               >  1:1423318/1‑69 (MQ=255)
                           gCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATgc                                                   >  1:1195563/1‑48 (MQ=255)
                                     gAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAgg                     >  1:2986479/1‑68 (MQ=255)
                                     gAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAgg                     >  1:3427274/1‑68 (MQ=255)
                                     gAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAgg                     >  1:1480454/1‑68 (MQ=255)
                                     gAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAgg                     >  1:600611/1‑68 (MQ=255)
                                     gAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAgg                     >  1:663170/1‑68 (MQ=255)
                                     gAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAgg                     >  1:3098724/1‑68 (MQ=255)
                                     gAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAgg                     >  1:1624991/1‑68 (MQ=255)
                                     gAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAgg                     >  1:3017162/1‑68 (MQ=255)
                                     gAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAgg                     >  1:2136300/1‑68 (MQ=255)
                                     gAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAgg                     >  1:241317/1‑68 (MQ=255)
                                     gAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAgg                     >  1:275285/1‑68 (MQ=255)
                                     gAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAgg                     >  1:2844536/1‑68 (MQ=255)
                                     gAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAgg                     >  1:3043167/1‑68 (MQ=255)
                                     gAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAgg                     >  1:1871134/1‑68 (MQ=255)
                                     caaCGGGTCGCCACCGCGCAAACGACAATCTTCCa                                                      >  1:2363111/2‑35 (MQ=255)
                                             caccACCGCGCAAACGGCAATCTTCCATGCGTTCTTTTAATACTGGATCATCCGGCGAGGTAAACGGCa            <  1:3216065/67‑1 (MQ=255)
                                             caccACCGCGCAAACGGCAATCTTCCATGCGTGCTTTTAATAGTGGATCATCCGGCGAGGTAAACGGCa            <  1:2399720/67‑1 (MQ=255)
                                             caccACCGCGCAAACGGCAATCTTCCATGCGTGCTTTTAATAGTGGATCATCCGGCGAGGTAAACGGCa            <  1:2006091/67‑1 (MQ=255)
                                             caccACCGCGCAAACGGCAATCTTCCATGCGTGCTTTTAATAGTGGATCATCCGGCGAGGTAAACGGCa            <  1:1421814/67‑1 (MQ=255)
                                             caccACCGCGCAAACGGCAATCTTCCATGCGTGCTTTTAATAGTGGATCATCCGGCGAGGTAAACGGCa            <  1:2199367/67‑1 (MQ=255)
                                             caccACCGCGCAAACGGCAATCTTCCATGCGTGCTTTTAATAGTGGATCATCCGGCGAGGTAAACGGCa            <  1:428971/67‑1 (MQ=255)
                                             caccACCGCGCAAACGGCAATCTTCCATGCGTGCTTTTAATAGTGGATCATCCGGCGAGGTAAACGGCa            <  1:1006284/67‑1 (MQ=255)
                                             caccACCGCGCAAACGGCAATCTTCCATGCGTGCTTTTAATAGTGGATCATCCGGCGAGGTAAACGGCa            <  1:2535322/67‑1 (MQ=255)
                                             caccACCGCGCAAACGGCAATCTTCCATGCGTGCTTTTAATAGTGGATCATCCGGCGAGGTAAACGGCa            <  1:1781118/67‑1 (MQ=255)
                                             caccACCGCGCAAACGGCAATCTTCCATGCGTGCTTTTAATAGTGGATCATCCGGCGAGGTAAACGGCa            <  1:1172921/67‑1 (MQ=255)
                                             caccACCGCGCAAACGGCAATCTTCCATGCGTGCTTTTAATAGTGGATCATCCGGCGAGGTAAACGGCa            <  1:3163385/67‑1 (MQ=255)
                                             caccACCGCGCAAACGGCAATCTTCCATGCGTGCTTTTAATAGTGGATCATCCGGCGAGGTAAACGGCa            <  1:2817698/67‑1 (MQ=255)
                                             caccACCGCGCAAACGGCAATCTTCCATGCGTGCTTTTAATAGTGGATCATCCGGCGAGGTAAACGGCa            <  1:3248214/67‑1 (MQ=255)
                                             cGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGc             <  1:2195419/68‑1 (MQ=255)
                                                    gcgcAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTc     >  1:1274295/1‑69 (MQ=255)
                                                    gcgcAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTc     >  1:607118/1‑69 (MQ=255)
                                                       cAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGAt  <  1:1565001/69‑1 (MQ=255)
                                                       cAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGAt  <  1:1751163/69‑1 (MQ=255)
                                                       cAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGAt  <  1:3257127/69‑1 (MQ=255)
                                                       cAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGAt  <  1:2094635/69‑1 (MQ=255)
                                                       cAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTAGTCGAt  <  1:441649/69‑1 (MQ=255)
                                                        aaaCGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGAt  <  1:660227/68‑1 (MQ=255)
                                                             |                                                              
TGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGAT  >  minE/1195420‑1195543

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: