Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 523,609 A→G 37.0% K37K (AAA→AAG moaA → molybdopterin biosynthesis protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE523,6090AG37.0% 55.0 / 63.5 100K37K (AAA→AAGmoaAmolybdopterin biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base A (26/37);  new base G (19/18);  total (45/55)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

Only 100 of 103 total aligned reads displayed.

CTTGCGCCTGTCGATTACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCG  >  minE/523543‑523668
                                                                  |                                                           
ctTGCGCCTGTCGATTACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACaaa                                                             >  1:841572/1‑67 (MQ=255)
ctTGCGCCTGTCGATTACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACaaa                                                             >  1:338430/1‑67 (MQ=255)
          tCGATTACTGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCTTc                                                 >  1:2043564/1‑69 (MQ=255)
          tCGATTACTGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTc                                                 >  1:2879815/1‑69 (MQ=255)
          tCGATTACTGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTc                                                 >  1:3248473/1‑69 (MQ=255)
          tCGATTACTGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTc                                                 >  1:636461/1‑69 (MQ=255)
          tCGATTACTGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTc                                                 >  1:1451989/1‑69 (MQ=255)
            gATTACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCa                                                >  1:1861858/1‑68 (MQ=255)
               tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:853067/63‑1 (MQ=255)
               tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:2242183/63‑1 (MQ=255)
               tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:2503464/63‑1 (MQ=255)
               tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:2562805/63‑1 (MQ=255)
               tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:2663842/63‑1 (MQ=255)
               tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:1872618/63‑1 (MQ=255)
               tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:795582/63‑1 (MQ=255)
               tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:1353801/63‑1 (MQ=255)
               tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:1346037/63‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCATTaaa                                        <  1:1313295/69‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:484209/69‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:2324128/69‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:2133298/69‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:2094253/69‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:752604/69‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:2505961/69‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:1993288/69‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:2581469/69‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:1299777/69‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:276969/69‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:1308428/69‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:322542/69‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:17863/69‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:60510/69‑1 (MQ=255)
                   gATGTGTGTAACTTTCGTTGCACATACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATata                                        <  1:1719190/69‑3 (MQ=255)
                   gATGTGTGTAACTTTCGTTGAACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:3392887/69‑1 (MQ=255)
                    aTGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:1788939/68‑1 (MQ=255)
                       tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGc                                     >  1:659130/1‑68 (MQ=255)
                       tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:3367257/1‑69 (MQ=255)
                       tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:2529525/1‑69 (MQ=255)
                       tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:3379604/1‑69 (MQ=255)
                       tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:643912/1‑69 (MQ=255)
                       tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:878097/1‑69 (MQ=255)
                       tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:1854174/1‑69 (MQ=255)
                       tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:1428249/1‑69 (MQ=255)
                       tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:1081209/1‑69 (MQ=255)
                       tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:199023/1‑69 (MQ=255)
                       tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    <  1:1069618/69‑1 (MQ=255)
                         tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAAcc                                                           >  1:313445/1‑44 (MQ=255)
                         tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCg                                                          >  1:2216754/1‑45 (MQ=255)
                         tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCg                                                          >  1:3361475/1‑45 (MQ=255)
                         tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCg                                                          >  1:2875379/1‑45 (MQ=255)
                         tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCg                                                          >  1:2382243/1‑45 (MQ=255)
                         tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCg                                                          >  1:2020901/1‑45 (MQ=255)
                         tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCttt                                  <  1:3323157/69‑1 (MQ=255)
                         tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCttt                                  <  1:3434051/69‑1 (MQ=255)
                         tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCttt                                  <  1:2311496/69‑1 (MQ=255)
                           tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTc                                 >  1:732674/1‑68 (MQ=255)
                           tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:1943368/1‑69 (MQ=255)
                           tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:3315451/1‑69 (MQ=255)
                           tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:3262387/1‑69 (MQ=255)
                           tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:2164530/1‑69 (MQ=255)
                           tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:2888113/1‑69 (MQ=255)
                           tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:178419/1‑69 (MQ=255)
                           tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:2342221/1‑69 (MQ=255)
                           tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:2785236/1‑69 (MQ=255)
                           tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:2720558/1‑69 (MQ=255)
                           tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:2595162/1‑69 (MQ=255)
                           tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:161076/1‑69 (MQ=255)
                           tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:559583/1‑69 (MQ=255)
                           tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:593194/1‑69 (MQ=255)
                             aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCTAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:1973863/68‑1 (MQ=255)
                             aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:3172592/68‑1 (MQ=255)
                             aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:1734674/68‑1 (MQ=255)
                             aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:1693156/68‑1 (MQ=255)
                             aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:2279934/68‑1 (MQ=255)
                             aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:3304073/68‑1 (MQ=255)
                             aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:66522/68‑1 (MQ=255)
                             aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:2741645/68‑1 (MQ=255)
                             aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:1990048/68‑1 (MQ=255)
                             aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:871775/68‑1 (MQ=255)
                             aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:2110516/68‑1 (MQ=255)
                             aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:2508044/68‑1 (MQ=255)
                             aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAACAAAGGCTTTCtt                               <  1:1647970/68‑1 (MQ=255)
                                              tgcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:423600/51‑1 (MQ=255)
                                               gcctgccGGATGGCTACAAACCTAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:3320237/61‑1 (MQ=255)
                                               gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:1274132/61‑1 (MQ=255)
                                               gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:2270900/61‑1 (MQ=255)
                                               gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:1123397/61‑1 (MQ=255)
                                               gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:191690/61‑1 (MQ=255)
                                               gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:1936672/61‑1 (MQ=255)
                                               gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:2697785/61‑1 (MQ=255)
                                               gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:2633675/61‑1 (MQ=255)
                                               gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:260272/61‑1 (MQ=255)
                                               gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGcc            >  1:192561/1‑69 (MQ=255)
                                               gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGcc            >  1:679959/1‑69 (MQ=255)
                                               gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGcc            >  1:2000804/1‑69 (MQ=255)
                                               gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGcc            >  1:10264/1‑69 (MQ=255)
                                               gcctgcaGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGcc            >  1:436070/1‑69 (MQ=255)
                                                         tGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTAcgcg  >  1:1123674/1‑69 (MQ=255)
                                                         tGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTAcgcg  >  1:2411911/1‑69 (MQ=255)
                                                                aaaCCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:1402292/44‑1 (MQ=255)
                                                                  |                                                           
CTTGCGCCTGTCGATTACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCG  >  minE/523543‑523668

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: