Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,611,483 G→A 47.3% Y85Y (TAC→TAT eutB ← ethanolamine ammonia‑lyase, large subunit, heavy chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,611,4830GA47.3% 8.9 / 79.5 93Y85Y (TAC→TATeutBethanolamine ammonia‑lyase, large subunit, heavy chain
Reads supporting (aligned to +/- strand):  ref base G (33/16);  new base A (16/28);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 3.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 106 total aligned reads displayed.

CCACGCTGGTTTCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGTCATCTTCATAGGCAAT  >  minE/1611416‑1611545
                                                                   |                                                              
ccACGCTGGTTTCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTg                                                                <  1:1847624/68‑1 (MQ=255)
ccACGCTGGTTTCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTg                                                                <  1:1755019/68‑1 (MQ=255)
ccACGCTGGTTTCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTg                                                                <  1:2505592/68‑1 (MQ=255)
ccACGCTGGTTTCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTg                                                                <  1:2407905/68‑1 (MQ=255)
ccACGCTGGTTTCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTg                                                                <  1:3390996/68‑1 (MQ=255)
ccACGCTGGTTTCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTg                                                                <  1:562776/68‑1 (MQ=255)
ccACGCTGGTTTCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTg                                                                <  1:2074611/68‑1 (MQ=255)
          ttCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTg                                                                >  1:1562858/1‑58 (MQ=255)
          ttCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTg                                                                >  1:2591843/1‑58 (MQ=255)
          ttCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTg                                                                >  1:1951084/1‑58 (MQ=255)
          ttCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTg                                                                >  1:1285326/1‑58 (MQ=255)
          ttCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTg                                                                >  1:674211/1‑58 (MQ=255)
          ttCATCGCTCAGCACATACTCACGAAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTg                                                                >  1:170382/1‑58 (MQ=255)
              tCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGtt                                                  <  1:1472781/68‑1 (MQ=255)
              tCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGtt                                                  <  1:2861422/68‑1 (MQ=255)
              tCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGtt                                                  <  1:3423543/68‑1 (MQ=255)
                            cTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGt                                                   >  1:2669403/1‑53 (MQ=255)
                            cTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGa                                     >  1:1516486/1‑67 (MQ=255)
                            cTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGa                                     >  1:2858829/1‑67 (MQ=255)
                             tCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATc                                  >  1:637778/1‑69 (MQ=255)
                             tCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATc                                  >  1:626305/1‑69 (MQ=255)
                             tCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATc                                  >  1:2856604/1‑69 (MQ=255)
                             tCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATc                                  >  1:792129/1‑69 (MQ=255)
                                  cAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGttt                                                      >  1:3367409/1‑44 (MQ=255)
                                  cAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGcc                              >  1:713560/1‑68 (MQ=255)
                                  cAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCcg                             >  1:1001210/1‑69 (MQ=255)
                                  cAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCcg                             >  1:1888368/1‑69 (MQ=255)
                                  cAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCcg                             >  1:3433322/1‑69 (MQ=255)
                                  cAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCcg                             >  1:1862650/1‑69 (MQ=255)
                                  cAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCcg                             >  1:2031906/1‑69 (MQ=255)
                                  cAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCcg                             >  1:2281395/1‑69 (MQ=255)
                                      tCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGt                          >  1:2073365/1‑68 (MQ=255)
                                      tCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGt                          >  1:2722794/1‑68 (MQ=255)
                                      tCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGt                          >  1:3009566/1‑68 (MQ=255)
                                      tCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGt                          >  1:165508/1‑68 (MQ=255)
                                      tCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGt                          >  1:3241459/1‑68 (MQ=255)
                                      tCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGt                          >  1:409482/1‑68 (MQ=255)
                                      tCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTc                         >  1:3164064/1‑69 (MQ=255)
                                      tCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTc                         >  1:1135682/1‑69 (MQ=255)
                                      tCGCTGATGCTCCAGTTTTTAATCTGCTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTc                         >  1:1442588/1‑69 (MQ=255)
                                      tCGCTGATGCTCCAGTATTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGt                          >  1:927425/1‑68 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTg                                      >  1:625858/1‑55 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTg                                      >  1:1425119/1‑55 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCa                        >  1:2428507/1‑69 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTATAGTCCGTTTCGTTAACGTCGTCCTGAATAAGCCGCGTCa                        >  1:1877841/1‑69 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTg                                      >  1:2047793/1‑55 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGt                          >  1:1803758/1‑67 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCa                        >  1:2021616/1‑69 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCa                        >  1:3120558/1‑69 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCa                        >  1:1324786/1‑69 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCa                        >  1:711088/1‑69 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCa                        >  1:1491843/1‑69 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCa                        >  1:1461228/1‑69 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCa                        >  1:1795550/1‑69 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCa                        >  1:461474/1‑69 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCa                        >  1:2299421/1‑69 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCa                        >  1:2369253/1‑69 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCa                        >  1:532370/1‑69 (MQ=255)
                                       cGCTGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTAGTCCTGa                                     >  1:7407/1‑56 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACg                      <  1:2705206/68‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACg                      <  1:2342528/68‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACg                      <  1:2624018/68‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACg                      <  1:1916946/68‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACg                      <  1:823746/68‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACg                      <  1:2287074/68‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACg                      <  1:3252365/68‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACg                      <  1:165254/68‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:1330568/69‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:381286/69‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:3430491/69‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:1302893/69‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:2073072/69‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:2406767/69‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:2902023/69‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:2900188/69‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:2896091/69‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:2889354/69‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:178549/69‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:895294/69‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:2714784/69‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:905863/69‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:2447060/69‑1 (MQ=255)
                                          tGATGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:2402589/69‑1 (MQ=255)
                                           gATGCACCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:1674573/68‑1 (MQ=255)
                                            aTGCTCCAGTTTTTAATCTGGTTATAGGCCGTTTCGTTAACGTCGTCCTGAATCAGCCGCGTCACGc                     <  1:1502039/67‑1 (MQ=255)
                                             tGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGt                  >  1:3138853/1‑69 (MQ=255)
                                             tGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGt                  >  1:1965560/1‑69 (MQ=255)
                                              gCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGTc                 <  1:511708/69‑1 (MQ=255)
                                              gCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGTc                 <  1:1352252/69‑1 (MQ=255)
                                              gCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGTc                 <  1:1455330/69‑1 (MQ=255)
                                              gCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGTc                 <  1:1500511/69‑1 (MQ=255)
                                              gCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGTc                 <  1:3122159/69‑1 (MQ=255)
                                              gCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGTc                 <  1:1723146/69‑1 (MQ=255)
                                              gCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGTc                 <  1:2725426/69‑1 (MQ=255)
                                              gCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGTc                 <  1:908752/69‑1 (MQ=255)
                                              gCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGTc                 <  1:948435/69‑1 (MQ=255)
                                               cTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGTc                 <  1:3288447/68‑1 (MQ=255)
                                                             cTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGTCATCTTCATAGGCAAt  <  1:774959/69‑1 (MQ=255)
                                                             cTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGTCATCTTCATAGGCAAt  <  1:2907700/69‑1 (MQ=255)
                                                             cTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGTCATCTTCATAGGCAAt  <  1:1905975/69‑1 (MQ=255)
                                                                   |                                                              
CCACGCTGGTTTCATCGCTCAGCACATACTCACGCAGTTCGCTGATGCTCCAGTTTTTAATCTGGTTGTAGGCCGTTTCGTTAACATCGTCCTGAATCAGCCGCGTCACGCAGTCATCTTCATAGGCAAT  >  minE/1611416‑1611545

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: