Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,740,334 C→T 6.7% V52I (GTC→ATC)  rseC ← RseC protein involved in reduction of the SoxR iron‑sulfur cluster

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,740,3340CT 6.7% 240.6 / 4.7 74V52I (GTC→ATC) rseCRseC protein involved in reduction of the SoxR iron‑sulfur cluster
Reads supporting (aligned to +/- strand):  ref base C (40/29);  new base T (3/2);  total (43/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTACCGCAACCGGCGCGTGAAGCGC  >  minE/1740264‑1740402
                                                                      |                                                                    
aGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGAc                                                                      <  1:3831668/71‑1 (MQ=255)
aGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGAc                                                                      <  1:2947659/71‑1 (MQ=255)
 gCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGAc                                                                      >  1:3103833/1‑70 (MQ=255)
  cAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACaa                                                                    >  1:981171/1‑71 (MQ=255)
     gCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGAc                                                                      <  1:2399382/66‑1 (MQ=255)
       tGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATgg                                                                 <  1:498157/69‑1 (MQ=255)
        gCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATg                                                                  >  1:1528173/1‑67 (MQ=255)
            tCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGt                                                          >  1:3830859/1‑71 (MQ=255)
            tCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGt                                                          >  1:3795279/1‑71 (MQ=255)
            tCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGt                                                          >  1:2944983/1‑71 (MQ=255)
            tCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGt                                                          >  1:213964/1‑71 (MQ=255)
            tCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGt                                                          >  1:209039/1‑71 (MQ=255)
                 gATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACaa                                                                    >  1:266225/1‑56 (MQ=255)
                   tCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGAc                                                                      <  1:1274541/52‑1 (MQ=255)
                    ccccAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACaa                                                                    >  1:1578999/1‑53 (MQ=255)
                     cccAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATgg                                                                 <  1:1138201/55‑1 (MQ=255)
                       cAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTa                                                               >  1:1477225/1‑55 (MQ=255)
                       cAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTa                                                               >  1:1585905/1‑55 (MQ=255)
                       cAATTCCACTTTTTGCCCCGGCACTAACGGCTCAACACAGGGTACGACAATGGTa                                                               >  1:3127404/1‑55 (MQ=255)
                        aaTTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTg                                                         <  1:3326666/60‑1 (MQ=255)
                          ttCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATgg                                                                 <  1:554199/50‑1 (MQ=255)
                            ccACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGt                                          >  1:1553139/1‑71 (MQ=255)
                            ccACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGt                                          >  1:3232346/1‑71 (MQ=255)
                            ccACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGt                                          >  1:3592484/1‑71 (MQ=255)
                            ccACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGt                                          >  1:47294/1‑71 (MQ=255)
                            ccACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGt                                          >  1:71040/1‑71 (MQ=255)
                            ccACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGt                                          >  1:2125995/1‑71 (MQ=255)
                              aCTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATgg                                                                 <  1:2077147/46‑1 (MQ=255)
                                tttttGCCCCGGCACTAACGGCTCATCACAGGGTACGAc                                                                      <  1:2642336/39‑1 (MQ=255)
                                tttttGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCgg                                                 >  1:2081139/1‑60 (MQ=255)
                                tttttGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGtt                                         <  1:3650560/68‑1 (MQ=255)
                                 ttttGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGc                                                   >  1:2251018/1‑57 (MQ=255)
                                 ttttGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCg                                                  >  1:3353274/1‑58 (MQ=255)
                                  tttGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGt                                                       >  1:4061469/1‑52 (MQ=255)
                                  tttGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGt                                                       >  1:657571/1‑52 (MQ=255)
                                  tttGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGtttatt                                     >  1:2960347/1‑70 (MQ=255)
                                  tttGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGtttatt                                     >  1:2961624/1‑70 (MQ=255)
                                  tttGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGtt                                         <  1:2835108/66‑1 (MQ=255)
                                   ttGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGtttattt                                    <  1:2104942/70‑1 (MQ=255)
                                   ttGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGtttattt                                    <  1:3347546/70‑1 (MQ=255)
                                   ttGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGtttattt                                    >  1:3279601/1‑70 (MQ=255)
                                     gCCCCGGCACTAACGGCTCATCACAGGGTACGACAATgg                                                                 <  1:574302/39‑1 (MQ=255)
                                     gCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTaa                                  <  1:2081510/70‑1 (MQ=255)
                                     gCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTaa                                  <  1:3554213/70‑1 (MQ=255)
                                        ccGGCACTAACGGCTCATCACAGGGTACGATAATGGTATGCGTgg                                                        <  1:1881157/45‑1 (MQ=255)
                                          ggCACTAACGGCTCATCACAGGGTACGATAATGGTATGCGTGGtt                                                      >  1:3403174/1‑45 (MQ=255)
                                            cACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGc                          <  1:137132/71‑1 (MQ=255)
                                              cTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGttt                                        <  1:231018/55‑1 (MQ=255)
                                              cTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGttt                                        <  1:3870200/55‑1 (MQ=255)
                                               tAACGGCTCATCACAGGGTACGATAATGGTATGCGTGGTTTGCGGGCCAAg                                           <  1:1500654/51‑1 (MQ=255)
                                               tAACGGCTCATCACAGGGTACGACAATGGTATGCg                                                           <  1:4110018/35‑1 (MQ=255)
                                               tAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGttt                                                     <  1:4056833/41‑1 (MQ=255)
                                               tAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAacac                               <  1:3397037/63‑1 (MQ=255)
                                                aaCGGCTCATCACAGGGTACGATAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTa                        >  1:3563116/1‑69 (MQ=255)
                                                aaCGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTAcc                      >  1:121553/1‑71 (MQ=255)
                                                aaCGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTAcc                      >  1:2462439/1‑71 (MQ=255)
                                                aaCGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTAcc                      >  1:126504/1‑71 (MQ=255)
                                                 aCGGCTCATCACAGGGTACGATAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTACCg                     >  1:2930776/1‑71 (MQ=255)
                                                 aCGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTAc                       <  1:55288/69‑1 (MQ=255)
                                                      tcatcaCAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGtttattt                                    <  1:3110359/51‑1 (MQ=255)
                                                       catcaCAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTACCGCAACCg               >  1:2755793/1‑71 (MQ=255)
                                                         tcaCAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTaa                                  <  1:864344/50‑1 (MQ=255)
                                                          cacaGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTACCGCAACCGgcg            >  1:1547593/1‑71 (MQ=255)
                                                          cacaGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTACCGCAACCGgcg            >  1:4112191/1‑71 (MQ=255)
                                                           acaGGGTACGACAATGGTATGCGTGGTTTGCGGGCCaa                                            >  1:3208204/1‑38 (MQ=255)
                                                           acaGGGTACGACAATGGTATGCGTGGTTTGCGGGCCaa                                            >  1:447171/1‑38 (MQ=255)
                                                           acaGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTaa                                  >  1:3172695/1‑48 (MQ=255)
                                                           acaGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTaa                                  >  1:1734028/1‑48 (MQ=255)
                                                           acaGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAaca                                >  1:631127/1‑50 (MQ=255)
                                                           acaGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTACCTCAACCGgcgc           >  1:2380141/1‑71 (MQ=255)
                                                           acaGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTACCGCAACCGgcgc           >  1:1696068/1‑71 (MQ=255)
                                                           acaGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTACCGCAACCGgcgc           >  1:3465804/1‑71 (MQ=255)
                                                           acaGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTACCGCAACCGgcgc           <  1:2082251/71‑1 (MQ=255)
                                                            caGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAg                                           <  1:3543068/38‑1 (MQ=255)
                                                            caGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAg                                           <  1:2907579/38‑1 (MQ=255)
                                                            caGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTACCGCAACCGGcgcg          <  1:696718/71‑1 (MQ=255)
                                                                 tACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTACCGCAACCGGCGCGTGAAg     <  1:1958434/71‑1 (MQ=255)
                                                                    gACAATGGTATGCGTGGTTTGCGGGCCAAGtttattt                                    <  1:49831/37‑1 (MQ=255)
                                                                    gACAATGGTATGCGTGGTTTGCGGGCCAAGtttattt                                    <  1:3785849/37‑1 (MQ=255)
                                                                    gACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTACCGCAACCGGCGCGTGAAgcgc  <  1:300467/71‑1 (MQ=255)
                                                                      cAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAaca                                >  1:3410694/1‑39 (MQ=255)
                                                                      |                                                                    
AGCAGGCTGCCTTCGGCGATCCCCAATTCCACTTTTTGCCCCGGCACTAACGGCTCATCACAGGGTACGACAATGGTATGCGTGGTTTGCGGGCCAAGTTTATTTAACACGCGGCTACCGCAACCGGCGCGTGAAGCGC  >  minE/1740264‑1740402

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: