Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,082,141 C→T 14.1% D10N (GAT→AAT)  ydhY ← predicted 4Fe‑4S ferridoxin‑type protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,082,1410CT14.1% 233.5 / 25.7 85D10N (GAT→AAT) ydhYpredicted 4Fe‑4S ferridoxin‑type protein
Reads supporting (aligned to +/- strand):  ref base C (41/32);  new base T (5/7);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

GGGCAATTGTTAGACCAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAACACAGTGCTGACAATTCAG  >  minE/1082072‑1082210
                                                                     |                                                                     
gggCAATTGTTAGACCAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCt                                                                      <  1:3804766/71‑1 (MQ=255)
gggCAATTGTTAGACCAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCt                                                                      <  1:2851478/71‑1 (MQ=255)
  gCAATTGTTAGACCAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTaa                                                                    <  1:1994131/71‑1 (MQ=255)
  gCAATTGTTAGACCAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTaa                                                                    <  1:3820272/71‑1 (MQ=255)
       tGTTAGACCAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAAt                                                                   >  1:2804225/1‑67 (MQ=255)
        gTTAGACCAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTg                                                               >  1:349029/1‑70 (MQ=255)
        gTTAGACCAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTg                                                               >  1:3599477/1‑70 (MQ=255)
         ttAGACCAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTgg                                                              <  1:3353392/70‑1 (MQ=255)
          tAGACCAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGAc                                                            >  1:1861018/1‑71 (MQ=255)
            gACCAGCAAGGCCTTTTCCTGATATACTCAAGAATTCTAATCTCGTTAAGCCAATATCTAATAGTGGACg                                                           <  1:3146443/70‑1 (MQ=255)
            gACCAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACg                                                           <  1:4018236/70‑1 (MQ=255)
            gACCAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACg                                                           <  1:355057/70‑1 (MQ=255)
            gACCAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACg                                                           <  1:2544462/70‑1 (MQ=255)
            gACCAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACg                                                           <  1:1881879/70‑1 (MQ=255)
              ccAGCATGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATc                                                        <  1:3352619/71‑1 (MQ=255)
               cAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATc                                                        >  1:607753/1‑70 (MQ=255)
               cAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCa                                                       <  1:1311815/71‑1 (MQ=255)
               cAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCa                                                       <  1:3845706/71‑1 (MQ=255)
                aGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAAt                                                                   <  1:4007652/58‑1 (MQ=255)
                aGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAAt                                                                   <  1:3997754/58‑1 (MQ=255)
                 gCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAAc                                                     >  1:3246963/1‑71 (MQ=255)
                       ccTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATTTAATa                                                                  <  1:3217444/52‑1 (MQ=255)
                       ccTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATTTAATa                                                                  <  1:1580808/52‑1 (MQ=255)
                       ccTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATTTAATa                                                                  <  1:3860631/52‑1 (MQ=255)
                       ccTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATTTAATa                                                                  <  1:2901919/52‑1 (MQ=255)
                       ccTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATTTAATa                                                                  <  1:696025/52‑1 (MQ=255)
                       ccTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATTTAATa                                                                  <  1:814016/52‑1 (MQ=255)
                       ccTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATTTAATa                                                                  <  1:1381397/52‑1 (MQ=255)
                         ttttCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCt                                                                      <  1:3365756/46‑1 (MQ=255)
                         ttttCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCt                                                                      <  1:3576402/46‑1 (MQ=255)
                         ttttCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCt                                                                      <  1:3632499/46‑1 (MQ=255)
                         ttttCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCt                                                                      <  1:2496270/46‑1 (MQ=255)
                         ttttCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCt                                                                      <  1:1640796/46‑1 (MQ=255)
                         ttttCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCt                                                                      <  1:4108497/46‑1 (MQ=255)
                         ttttCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCa                                             >  1:4079654/1‑71 (MQ=255)
                         ttttCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCa                                             >  1:1626547/1‑71 (MQ=255)
                         ttttCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCa                                             >  1:2580571/1‑71 (MQ=255)
                           ttCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCatt                                           >  1:3205757/1‑71 (MQ=255)
                           ttCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCatt                                           >  1:551502/1‑71 (MQ=255)
                             ccTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGAc                                                            >  1:1767071/1‑52 (MQ=255)
                                gATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCatttattt                                      >  1:239378/1‑71 (MQ=255)
                                gATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCatttattt                                      >  1:670513/1‑71 (MQ=255)
                                  tataCGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTc                                     <  1:2311962/70‑1 (MQ=255)
                                  tataCGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCg                                    <  1:3626070/71‑1 (MQ=255)
                                    taCGCAAGAATTCTAATCGCGTTAAGCCAATATCTaa                                                                    <  1:3062011/37‑1 (MQ=255)
                                    taCGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCa                                  >  1:2142023/1‑71 (MQ=255)
                                      cGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATc                                >  1:3048768/1‑71 (MQ=255)
                                      cGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATc                                >  1:3884132/1‑71 (MQ=255)
                                      cGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATc                                <  1:677083/71‑1 (MQ=255)
                                        cAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCct                              >  1:706982/1‑71 (MQ=255)
                                        cAAGAATTCTAATCGCGTTAACCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCct                              >  1:1110470/1‑71 (MQ=255)
                                          aGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCctct                            >  1:1948041/1‑71 (MQ=255)
                                                cTAATCGCGTTAAGCCAATATTTAATAGTGGACGATCAACCgg                                                  >  1:6541/1‑43 (MQ=255)
                                                cTAATCGCGTTAAGCCAATATTTAATAGTGGACGATCAACCgg                                                  >  1:4009360/1‑43 (MQ=255)
                                                cTAATCGCGTTAAGCCAATATTTAATAGTGGACGATCAACCgg                                                  >  1:865532/1‑43 (MQ=255)
                                                cTAATCGCGTTAAGCCAATATTTAATAGTGGACGATCAACCgg                                                  >  1:3494149/1‑43 (MQ=255)
                                                cTAATCGCGTTAAGCCAATATTTAATAGTGGACGATCAACCgg                                                  >  1:2926342/1‑43 (MQ=255)
                                                cTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTc                                     >  1:2130272/1‑56 (MQ=255)
                                                cTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTc                                     >  1:630859/1‑56 (MQ=255)
                                                 tAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATc                                <  1:1607966/60‑1 (MQ=255)
                                                 tAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATc                                <  1:2482723/60‑1 (MQ=255)
                                                   aTCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAaca                   >  1:1135006/1‑71 (MQ=255)
                                                   aTCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAaca                   <  1:1550976/71‑1 (MQ=255)
                                                   aTCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAaca                   >  1:542832/1‑71 (MQ=255)
                                                    tCGCGTTAAGCCAATATCTAATAGTGTACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAacac                  <  1:3857598/71‑1 (MQ=255)
                                                    tCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCa                       <  1:213364/66‑1 (MQ=255)
                                                      gcgTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCctc                             >  1:336215/1‑58 (MQ=255)
                                                      gcgTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCtt                           >  1:4000040/1‑60 (MQ=255)
                                                      gcgTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCtt                           >  1:3802596/1‑60 (MQ=255)
                                                         ttAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCa                                  <  1:2179452/50‑1 (MQ=255)
                                                         ttAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCa                                  <  1:3632213/50‑1 (MQ=255)
                                                          tAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAACACAGTGCt            >  1:1810585/1‑71 (MQ=255)
                                                           aaGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCaaa                     >  1:2146887/1‑61 (MQ=255)
                                                           aaGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCaaa                     >  1:2202369/1‑61 (MQ=255)
                                                           aaGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCaaa                     >  1:651830/1‑61 (MQ=255)
                                                                    tCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCaaa                     >  1:38898/1‑52 (MQ=255)
                                                                    tCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAACACAGTGCTGAc         >  1:2137016/1‑64 (MQ=255)
                                                                    tCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAACACAGTGCTGACAATTcag  >  1:4045619/1‑71 (MQ=255)
                                                                    tCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAACACAGTGCTGACAATTcag  >  1:2498443/1‑71 (MQ=255)
                                                                    tCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAACACAGTGCTGACAATTcag  >  1:649402/1‑71 (MQ=255)
                                                                    tCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAACACAGTGCTGACAATTcag  >  1:3033111/1‑71 (MQ=255)
                                                                    tCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAACACAGTGCTGACAATTcag  >  1:1574897/1‑71 (MQ=255)
                                                                    tCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAACACAGTGCTGACAATTcag  >  1:1526855/1‑71 (MQ=255)
                                                                    tCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAACACAGTGCTGACAATTca   >  1:451767/1‑70 (MQ=255)
                                                                    tCTAATAGTGGACGATCAACCGGGTTAatttattt                                      >  1:2700722/1‑35 (MQ=255)
                                                                     cTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAacac                  <  1:3858061/54‑1 (MQ=255)
                                                                     |                                                                     
GGGCAATTGTTAGACCAGCAAGGCCTTTTCCTGATATACGCAAGAATTCTAATCGCGTTAAGCCAATATCTAATAGTGGACGATCAACCGGGTTCATTTATTTCGCATCCTCTTCGCAAACACAGTGCTGACAATTCAG  >  minE/1082072‑1082210

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: