Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 739,219 C→T 29.5% A530T (GCT→ACT)  fhuE ← ferric‑rhodotorulic acid outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE739,2190CT29.5% 92.2 / 39.2 92A530T (GCT→ACT) fhuEferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (20/44);  new base T (6/21);  total (26/66)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGC  >  minE/739157‑739279
                                                              |                                                            
cAGACGGCTATTCATCCAGTCCGATTTCAGACGCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGat                                                         <  1:2237713/68‑1 (MQ=255)
cAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACGGGTGATTGGAGCCAGat                                                         <  1:3228232/68‑1 (MQ=255)
cAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGata                                                        <  1:2796813/69‑1 (MQ=255)
cAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGata                                                        <  1:900883/69‑1 (MQ=255)
cAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGata                                                        <  1:3246615/69‑1 (MQ=255)
cAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGat                                                         <  1:2812247/68‑1 (MQ=255)
cAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGat                                                         <  1:3015586/68‑1 (MQ=255)
cAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGat                                                         <  1:3179915/68‑1 (MQ=255)
 aGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCTTAGTTGTTACCGGTGATTGTAGCCAGat                                                         <  1:302404/67‑1 (MQ=255)
         aTTCATCCAGTCCGATTTCAGACCCAGGTCGTAGTTGTTACCGGTGATTGGAGTCAGATATTTGCCTGa                                               >  1:842695/1‑69 (MQ=255)
         aTTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGTCAGATATTTGCCTg                                                >  1:936668/1‑68 (MQ=255)
         aTTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGTCAGATATTTGCCTGa                                               >  1:2244985/1‑69 (MQ=255)
         aTTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGTCAGATATTTGCCTGa                                               >  1:847883/1‑69 (MQ=255)
         aTTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGTCAGATATTTGCCTGa                                               >  1:3007428/1‑69 (MQ=255)
         aTTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGTCAGATATTTGCCTGa                                               >  1:1651534/1‑69 (MQ=255)
                cAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTc                                         >  1:3335276/1‑68 (MQ=255)
                cAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTc                                         >  1:979564/1‑68 (MQ=255)
                cAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTc                                         >  1:1853433/1‑68 (MQ=255)
                cAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTc                                         >  1:1714325/1‑68 (MQ=255)
                cAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTc                                         >  1:1487318/1‑68 (MQ=255)
                cAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTc                                         >  1:3362385/1‑68 (MQ=255)
                cAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTc                                         >  1:1328277/1‑68 (MQ=255)
                               cccAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGc                         <  1:2898520/69‑1 (MQ=255)
                               cccAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGc                         <  1:2905731/69‑1 (MQ=255)
                               cccAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGc                         <  1:1050185/69‑1 (MQ=255)
                               cccAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGc                         <  1:3411348/69‑1 (MQ=255)
                               cccAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGc                         <  1:2794004/69‑1 (MQ=255)
                               cccAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGc                         <  1:3402095/69‑1 (MQ=255)
                               cccAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGc                         <  1:2790096/69‑1 (MQ=255)
                               cccAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGc                         <  1:1316827/69‑1 (MQ=255)
                                    cTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGaa                    <  1:21686/69‑1 (MQ=255)
                                      cGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGGCTGa                                               <  1:3316653/40‑1 (MQ=255)
                                      cGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGa                                               <  1:1570544/40‑1 (MQ=255)
                                      cGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGa                                               <  1:3004109/40‑1 (MQ=255)
                                      cGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGa                                               <  1:709036/40‑1 (MQ=255)
                                         aGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATat               <  1:1776899/69‑2 (MQ=255)
                                         aGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATag               <  1:2254841/69‑1 (MQ=255)
                                         aGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATag               <  1:1120393/69‑1 (MQ=255)
                                         aGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATag               <  1:2611423/69‑1 (MQ=255)
                                         aGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATag               <  1:2654630/69‑1 (MQ=255)
                                         aGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATag               <  1:711988/69‑1 (MQ=255)
                                         aGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATag               <  1:1366789/69‑1 (MQ=255)
                                         aGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATag               <  1:348450/69‑1 (MQ=255)
                                         aGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATa                <  1:1032852/68‑1 (MQ=255)
                                          gttgttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATag               <  1:2746927/68‑1 (MQ=255)
                                            tgttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg            <  1:3000209/69‑1 (MQ=255)
                                            tgttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg            <  1:13221/69‑1 (MQ=255)
                                            tgttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg            <  1:1375689/69‑1 (MQ=255)
                                            tgttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg            <  1:3430374/69‑1 (MQ=255)
                                            tgttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg            <  1:2817257/69‑1 (MQ=255)
                                            tgttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg            <  1:618457/69‑1 (MQ=255)
                                            tgttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg            <  1:846970/69‑1 (MQ=255)
                                            tgttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg            <  1:2489577/69‑1 (MQ=255)
                                             gttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg            >  1:3393086/1‑68 (MQ=255)
                                             gttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg            >  1:1232695/1‑68 (MQ=255)
                                             gttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg            >  1:2095748/1‑68 (MQ=255)
                                             gttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg            >  1:252781/1‑68 (MQ=255)
                                             gttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg            >  1:2556158/1‑68 (MQ=255)
                                             gttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg            >  1:1793371/1‑68 (MQ=255)
                                             gttACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg            >  1:2560110/1‑68 (MQ=255)
                                                aCCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGtat         <  1:1531837/68‑1 (MQ=255)
                                                aCCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGtat         <  1:3323238/68‑1 (MQ=255)
                                                   ggTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTGTAGCt     <  1:2908996/69‑1 (MQ=255)
                                                    gTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGc                         <  1:1079661/48‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:2716324/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:2862964/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:3036523/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:1833541/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:199973/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:128798/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:844346/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:1231105/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:230880/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:326492/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:3116025/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:1441403/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:1452347/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:324977/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:3271228/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:3331189/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:15938/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:3362653/69‑1 (MQ=255)
                                                      gATTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:3371959/69‑1 (MQ=255)
                                                      gATTGGAGGCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGGATAGCTGGc  <  1:2847131/69‑1 (MQ=255)
                                                       aTTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:1616778/68‑1 (MQ=255)
                                                       aTTGGAGTCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGc  <  1:3425790/68‑1 (MQ=255)
                                                       aTTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTgg   >  1:420801/1‑67 (MQ=255)
                                                       aTTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTgg   >  1:1436130/1‑67 (MQ=255)
                                                       aTTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTgg   >  1:863512/1‑67 (MQ=255)
                                                       aTTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTgg   >  1:918262/1‑67 (MQ=255)
                                                       aTTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTgg   >  1:3150560/1‑67 (MQ=255)
                                                       aTTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTgg   >  1:951958/1‑67 (MQ=255)
                                                              |                                                            
CAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGC  >  minE/739157‑739279

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: