Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 740,069 C→G 34.0% S246S (TCG→TCC fhuE ← ferric‑rhodotorulic acid outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE740,0690CG34.0% 97.2 / 74.0 97S246S (TCG→TCCfhuEferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (49/15);  new base G (30/3);  total (79/18)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

AATGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAGTCGTTATTCTGGTAGCCGCCGACAA  >  minE/740003‑740136
                                                                  |                                                                   
aaggcgcTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:1461141/4‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:186307/1‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:1730249/1‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:3093899/1‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:1768701/1‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:3176334/1‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:1827986/1‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:1830411/1‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:3330765/1‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:1381926/1‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:1268540/1‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:525542/1‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:2469168/1‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:2291003/1‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:2092162/1‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:897495/1‑68 (MQ=255)
aaTGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCg                                                                    >  1:2154517/1‑68 (MQ=255)
    cgcTGATATTCGTAACCGGCTGAAAGCGTGGTAAGATCGCCTAAATCAGCATCGACAATGCCCgaga                                                                 >  1:2006449/1‑67 (MQ=255)
    cgcTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCgaga                                                                 >  1:1584982/1‑67 (MQ=255)
            ttCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCtt                                                       <  1:997458/69‑1 (MQ=255)
            ttCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCtt                                                       <  1:2540824/69‑1 (MQ=255)
            ttCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCtt                                                       <  1:1469272/69‑1 (MQ=255)
            ttCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCtt                                                       <  1:668248/69‑1 (MQ=255)
             tCGTACCCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATAGACAATGCCCGAGAAGAAGGTCtt                                                       <  1:243381/68‑1 (MQ=255)
             tCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCtt                                                       <  1:248952/68‑1 (MQ=255)
              cGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCtttt                                                     >  1:2177172/1‑69 (MQ=255)
              cGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCtttt                                                     >  1:2673882/1‑69 (MQ=255)
              cGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCtttt                                                     >  1:1068031/1‑69 (MQ=255)
              cGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCtttt                                                     >  1:3406992/1‑69 (MQ=255)
                tAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTc                                                    >  1:1062725/1‑68 (MQ=255)
                tAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTc                                                    >  1:1358079/1‑68 (MQ=255)
                tAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCa                                                   >  1:2370078/1‑69 (MQ=255)
                tAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCa                                                   >  1:3053779/1‑69 (MQ=255)
                tAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCa                                                   >  1:2393654/1‑69 (MQ=255)
                tAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCa                                                   >  1:663876/1‑69 (MQ=255)
                tAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCa                                                   >  1:2613031/1‑69 (MQ=255)
                tAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCa                                                   >  1:330973/1‑69 (MQ=255)
                tAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCa                                                   >  1:3004350/1‑69 (MQ=255)
                tAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCa                                                   >  1:3160445/1‑69 (MQ=255)
                tAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCa                                                   >  1:3031922/1‑69 (MQ=255)
                tAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCa                                                   >  1:1054494/1‑69 (MQ=255)
                tAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCa                                                   >  1:2067173/1‑69 (MQ=255)
                tAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCa                                                   >  1:1632631/1‑69 (MQ=255)
                tAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCa                                                   >  1:19716/1‑69 (MQ=255)
                       cTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCACtgttgt                                            <  1:3216085/69‑1 (MQ=255)
                       cTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCACtgttgt                                            <  1:2258492/69‑1 (MQ=255)
                       cTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCACtgttgt                                            <  1:3378050/69‑1 (MQ=255)
                       cTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCACtgttgt                                            <  1:2514273/69‑1 (MQ=255)
                       cTGAAAGCGTAGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCACtgttgt                                            <  1:1395274/69‑1 (MQ=255)
                            aGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCttt                                                      >  1:1930645/1‑54 (MQ=255)
                            aGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCACTGTTGTAGCgg                                       >  1:440693/1‑69 (MQ=255)
                            aGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCACTGTTGTAGCgg                                       >  1:1732813/1‑69 (MQ=255)
                                     aGATCGCCTAAATCAGCATCGACAATGCCGGAGaagaag                                                            >  1:2005850/1‑39 (MQ=255)
                                     aGATCGCCTAAATCAGCATCGACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGcc                                >  1:2853565/1‑67 (MQ=255)
                                     aGATCGCCTAAATCAGCATCGACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGcc                                >  1:83737/1‑67 (MQ=255)
                                     aGATCGCCTAAATCAGCATCGACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGcc                                >  1:3017761/1‑67 (MQ=255)
                                     aGATCGCCTAAATCAGCATCGACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGcc                                >  1:1804617/1‑67 (MQ=255)
                                     aGATCGCCTAAATCAGCATCGACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGcc                                >  1:2452571/1‑67 (MQ=255)
                                     aGATCGCCTAAATCAGCATCGACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGcc                                >  1:1361172/1‑67 (MQ=255)
                                     aGATCGCCTAAATCAGCATCGACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGcc                                >  1:2683093/1‑67 (MQ=255)
                                     aGATCGCCTAAATCAGCATCGACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGcc                                >  1:2645461/1‑67 (MQ=255)
                                     aGATCGCCTAAATCAGCATCGACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGcc                                >  1:2625917/1‑67 (MQ=255)
                                     aGATCGCCTAAATCAGCATCGACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGcc                                >  1:1831858/1‑67 (MQ=255)
                                     aGATCGCCTAAATCAGCATCGACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGc                                 >  1:2749172/1‑66 (MQ=255)
                                     aGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGcc                                >  1:3246177/1‑67 (MQ=255)
                                     aGATCGCCTAAATCAGCATCGACAATGCAGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGcc                                >  1:2041321/1‑67 (MQ=255)
                                         cGCCTAAATCAGCATCGAGAATGCCCGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAGt                          >  1:1204707/1‑69 (MQ=255)
                                         cGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAg                           >  1:1417215/1‑68 (MQ=255)
                                         cGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAGt                          >  1:3010033/1‑69 (MQ=255)
                                         cGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAGt                          >  1:1904417/1‑69 (MQ=255)
                                            cTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAGTCGt                       <  1:2960048/69‑1 (MQ=255)
                                            cTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAGTCGt                       <  1:2221292/69‑1 (MQ=255)
                                                     cATCGACAATGCCCGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAGTCGTTATTCTGGt              <  1:747049/69‑1 (MQ=255)
                                                     cATCGACAATGCCCGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAGTCGTTATTCTGGt              <  1:444441/69‑1 (MQ=255)
                                                         taCAATGCCGGAGAAGAAGGTCTTTTAACTGTTGTAGCGGTCCAGCCa                               >  1:3053840/2‑48 (MQ=38)
                                                         gACTATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:1636070/1‑49 (MQ=39)
                                                         gACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:1163294/1‑49 (MQ=255)
                                                         gACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:105176/1‑49 (MQ=255)
                                                         gACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:733693/1‑49 (MQ=255)
                                                         gACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:355282/1‑49 (MQ=255)
                                                         gACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:1526222/1‑49 (MQ=255)
                                                         gACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:3047244/1‑49 (MQ=255)
                                                         gACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:2779213/1‑49 (MQ=255)
                                                         gACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:1813273/1‑49 (MQ=255)
                                                         gACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:2093907/1‑49 (MQ=255)
                                                         gACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:2300653/1‑49 (MQ=255)
                                                         gACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:2473706/1‑49 (MQ=255)
                                                         gACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:2343505/1‑49 (MQ=255)
                                                         gACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:1520715/1‑49 (MQ=255)
                                                         gACAATGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:1046118/1‑49 (MQ=255)
                                                         gACAATGCCGGACAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCAt                              >  1:2024373/1‑49 (MQ=255)
                                                            aaTGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAGTCGTTATTCTGGTAgccgcc       <  1:451213/69‑1 (MQ=255)
                                                            aaTGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAGTCGTTATTCTGGTAgccgcc       <  1:356382/69‑1 (MQ=255)
                                                            aaTGCCGGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAGTCGTTATTCTGGTAgccgcc       <  1:2104338/69‑1 (MQ=255)
                                                            aaTGCCCGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAGTCGTTATTCTGGTAgccgcc       >  1:1133968/1‑69 (MQ=255)
                                                               gCCCGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAGTCGTTAtt                   >  1:2634454/1‑54 (MQ=255)
                                                               gCCCGAGAAGAAGGTCTTTTCACTGCTGTAGCGGTCCAGCCATGAGTCGTTAtt                   >  1:2360853/1‑54 (MQ=255)
                                                                  cGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAGTCGTTATTCTGGTAGCCGCCGACaa  <  1:688652/68‑1 (MQ=255)
                                                                  |                                                                   
AATGCGCTGATATTCGTAACCGGCTGAAAGCGTCGTAAGATCGCCTAAATCAGCATCGACAATGCCCGAGAAGAAGGTCTTTTCACTGTTGTAGCGGTCCAGCCATGAGTCGTTATTCTGGTAGCCGCCGACAA  >  minE/740003‑740136

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: