Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 523,603 C→T 35.5% G35G (GGC→GGT moaA → molybdopterin biosynthesis protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE523,6030CT35.5% 75.2 / 60.3 104G35G (GGC→GGTmoaAmolybdopterin biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base C (30/37);  new base T (19/18);  total (49/55)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 111 total aligned reads displayed.

TTTTACTACTTGCGCCTGTCGATTACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCG  >  minE/523535‑523668
                                                                    |                                                                 
ttttACTACTTGCGCCTGTCGATTACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGc                                                                   <  1:2501686/69‑1 (MQ=255)
ttttACTACTTGCGCCTGTCGATTACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGc                                                                   <  1:3328241/69‑1 (MQ=255)
ttttACTACTTGCGCCTGTCGATTACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGc                                                                   <  1:283932/69‑1 (MQ=255)
ttttACTACTTGCGCCTGTCGATTACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGc                                                                   <  1:641197/69‑1 (MQ=255)
ttttACTACTTGCGCCTGTCGATTACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGc                                                                   <  1:1578733/69‑1 (MQ=255)
        ctTGCGCCTGTCGATTACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACaaa                                                             >  1:841572/1‑67 (MQ=255)
        ctTGCGCCTGTCGATTACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACaaa                                                             >  1:338430/1‑67 (MQ=255)
                  tCGATTACTGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCTTc                                                 >  1:2043564/1‑69 (MQ=255)
                  tCGATTACTGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTc                                                 >  1:636461/1‑69 (MQ=255)
                  tCGATTACTGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTc                                                 >  1:3248473/1‑69 (MQ=255)
                  tCGATTACTGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTc                                                 >  1:2879815/1‑69 (MQ=255)
                  tCGATTACTGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTc                                                 >  1:1451989/1‑69 (MQ=255)
                    gATTACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCa                                                >  1:1861858/1‑68 (MQ=255)
                       tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:2562805/63‑1 (MQ=255)
                       tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:1346037/63‑1 (MQ=255)
                       tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:853067/63‑1 (MQ=255)
                       tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:795582/63‑1 (MQ=255)
                       tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:2663842/63‑1 (MQ=255)
                       tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:2503464/63‑1 (MQ=255)
                       tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:1872618/63‑1 (MQ=255)
                       tACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGt                                                  <  1:2242183/63‑1 (MQ=255)
                           gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCATTaaa                                        <  1:1313295/69‑1 (MQ=255)
                           gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:1993288/69‑1 (MQ=255)
                           gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:2505961/69‑1 (MQ=255)
                           gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:276969/69‑1 (MQ=255)
                           gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:2094253/69‑1 (MQ=255)
                           gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:2581469/69‑1 (MQ=255)
                           gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:2133298/69‑1 (MQ=255)
                           gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:2324128/69‑1 (MQ=255)
                           gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:322542/69‑1 (MQ=255)
                           gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:484209/69‑1 (MQ=255)
                           gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:60510/69‑1 (MQ=255)
                           gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:1308428/69‑1 (MQ=255)
                           gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:1299777/69‑1 (MQ=255)
                           gATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATaaa                                        <  1:752604/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)
                               tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTa                                                                 >  1:524093/1‑40 (MQ=255)
                               tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:3367257/1‑69 (MQ=255)
                               tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:3379604/1‑69 (MQ=255)
                               tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:643912/1‑69 (MQ=255)
                               tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:2529525/1‑69 (MQ=255)
                               tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:878097/1‑69 (MQ=255)
                               tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    <  1:1069618/69‑1 (MQ=255)
                               tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:1428249/1‑69 (MQ=255)
                               tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:1854174/1‑69 (MQ=255)
                               tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:199023/1‑69 (MQ=255)
                               tgtgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCt                                    >  1:1081209/1‑69 (MQ=255)
                                 tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGc                                                                   <  1:2100578/36‑1 (MQ=255)
                                 tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACa                                                               >  1:2838590/1‑40 (MQ=255)
                                 tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAAcc                                                           >  1:313445/1‑44 (MQ=255)
                                 tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCg                                                          >  1:3361475/1‑45 (MQ=255)
                                 tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCg                                                          >  1:2216754/1‑45 (MQ=255)
                                 tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCg                                                          >  1:2020901/1‑45 (MQ=255)
                                 tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCttt                                  <  1:2311496/69‑1 (MQ=255)
                                 tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCttt                                  <  1:3434051/69‑1 (MQ=255)
                                 tgtAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCttt                                  <  1:3323157/69‑1 (MQ=255)
                                   tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTc                                 >  1:732674/1‑68 (MQ=255)
                                   tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:2720558/1‑69 (MQ=255)
                                   tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:1943368/1‑69 (MQ=255)
                                   tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:2164530/1‑69 (MQ=255)
                                   tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:178419/1‑69 (MQ=255)
                                   tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:3262387/1‑69 (MQ=255)
                                   tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:559583/1‑69 (MQ=255)
                                   tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:593194/1‑69 (MQ=255)
                                   tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:2595162/1‑69 (MQ=255)
                                   tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:2888113/1‑69 (MQ=255)
                                   tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:161076/1‑69 (MQ=255)
                                   tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:2342221/1‑69 (MQ=255)
                                   tAACTTTCGTTGCACCTACTGCCTGCCGGATGGTTACAAGCCGAGCGGCGTCACCAATAAAGGCTTTCt                                >  1:2785236/1‑69 (MQ=255)
                                     aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCTAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:1973863/68‑1 (MQ=255)
                                     aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:1693156/68‑1 (MQ=255)
                                     aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:1734674/68‑1 (MQ=255)
                                     aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:66522/68‑1 (MQ=255)
                                     aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:2508044/68‑1 (MQ=255)
                                     aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:2279934/68‑1 (MQ=255)
                                     aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:2741645/68‑1 (MQ=255)
                                     aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:3304073/68‑1 (MQ=255)
                                     aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:2110516/68‑1 (MQ=255)
                                     aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:3172592/68‑1 (MQ=255)
                                     aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCtt                               <  1:871775/68‑1 (MQ=255)
                                     aCTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAACAAAGGCTTTCtt                               <  1:1647970/68‑1 (MQ=255)
                                                       gcctgccGGATGGCTACAAACCTAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:3320237/61‑1 (MQ=255)
                                                       gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:1936672/61‑1 (MQ=255)
                                                       gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:1123397/61‑1 (MQ=255)
                                                       gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:191690/61‑1 (MQ=255)
                                                       gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:2270900/61‑1 (MQ=255)
                                                       gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:1274132/61‑1 (MQ=255)
                                                       gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:3130544/61‑1 (MQ=255)
                                                       gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:2633675/61‑1 (MQ=255)
                                                       gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGa                    <  1:2697785/61‑1 (MQ=255)
                                                       gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGcc            >  1:10264/1‑69 (MQ=255)
                                                       gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGcc            >  1:679959/1‑69 (MQ=255)
                                                       gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGcc            >  1:192561/1‑69 (MQ=255)
                                                       gcctgccGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGcc            >  1:2000804/1‑69 (MQ=255)
                                                       gcctgcaGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGcc            >  1:436070/1‑69 (MQ=255)
                                                                 tGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGt                          >  1:2005041/1‑45 (MQ=255)
                                                                 tGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTAcgcg  >  1:1123674/1‑69 (MQ=255)
                                                                 tGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTAcgcg  >  1:2411911/1‑69 (MQ=255)
                                                                    |                                                                 
TTTTACTACTTGCGCCTGTCGATTACCGATGTGTGTAACTTTCGTTGCACCTACTGCCTGCCGGATGGCTACAAACCGAGCGGCGTCACCAATAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCG  >  minE/523535‑523668

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: